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ABSTRACT 



We present a comprehensive photochemistry model for exploration of the 
chemical composition of terrestrial exoplanet atmospheres. The photochemistry 
model is designed from the ground up to have the capacity to treat all types of 
terrestrial planet atmospheres, ranging from oxidizing through reducing, which 
makes the code suitable for applications for the wide range of anticipated ter- 
restrial exoplanet compositions. The one-dimensional chemical transport model 
treats up to 800 chemical reactions, photochemical processes, dry and wet deposi- 
tion, surface emission and thermal escape of O, H, C, N and S bearing species, as 
well as formation and deposition of elemental sulfur and sulfuric acid aerosols. We 
validate the model by computing the atmospheric composition of current Earth 
and Mars and find agreement with observations of major trace gases in Earth's 
and Mars' atmospheres. We simulate several plausible atmospheric scenarios of 
terrestrial exoplanets, and choose three benchmark cases for atmospheres from 
reducing to oxidizing. The most interesting finding is that atomic hydrogen is 
always a more abundant reactive radical than the hydroxyl radical in anoxic at- 
mospheres. Whether atomic hydrogen is the most important removal path for a 
molecule of interest also depends on the relevant reaction rates. We also find that 
volcanic carbon compounds (i.e., CH 4 and C0 2 ) are chemically long-lived and 
tend to be well mixed in both reducing and oxidizing atmospheres, and their dry 
deposition velocities to the surface control the atmospheric oxidation states. Fur- 
thermore, we revisit whether photochemically produced oxygen can cause false 
positives for detecting oxygenic photosynthesis, and find that in 1-bar C0 2 -rich 
atmospheres oxygen and ozone may build up to levels that have been previously 
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considered unique signatures of life, if there is no surface emission of reducing 
gases. The atmospheric scenarios present in this paper can serve as the bench- 
mark atmospheres for quickly assessing the lifetime of trace gases in reducing, 
weakly oxidizing, and highly oxidizing atmospheres on terrestrial exoplanets for 
the exploration of possible biosignature gases. 

Subject headings: radiative transfer — atmospheric effects — planetary systems 
— techniques: spectroscopic — astrobiology 

1. Introduction 

The search and characterization of super Earths is flourishing. A large number of super 
Earths are being detected (e.g., Rivera et al. 2005; Udry et al. 2007; Forveille et al. 2009; 
Leger et al. 2009; Charbonneau et al. 2009; Mayor et al. 2009; Vogt et al. 2010; Rivera 
et al. 2010; Dawson & Fabrycky 2010; Holman et al. 2010; Howard et al. 2011; Bonfils et 
al. 2011; Winn et al. 2011; Demory et al. 2011; Batalha et al. 2011; Lissauer et al., 2011; 
Borucki et al. 2011; Fressin et al 2012; Cochran et al. 2012). Attempts to observe super 
Earth atmospheres are growing (e.g., Batahla et al. 2011 for Kepler 10 b; Gillon et al. 2012 
and Demory et al. 2012 for 55 Cnc e), and one super Earth/mini Neptune GJ 1214b is being 
observed in as much detail as possible (e.g., Bean et al. 2010; Croll et al. 2011; Desert et al. 
2011; Berta et al. 2012). Although we distinguish between super Earths and mini Neptunes 
theoretically (super Earths are rocky planets with thin atmospheres whereas mini Neptunes 
are planets with massive gas envelopes), it is so far difficult to discriminate between the 
two categories observationally. Nonetheless, we focus on terrestrial planet atmospheres that 
include super Earths, since those planets hold the most interest on the path to finding and 
characterizing planets that might harbor life. In addition, the observational push to super 
Earth characterization has the potential to provide a handful of super Earth atmospheres 
to study in the coming years. In the more distant future, the community still holds hope 
that a direct-imaging space-based mission under the Terrestrial Planet Finder concept will 
enable planets to be observed in reflected light. 

Super Earth atmospheres require understanding of photochemistry for their study. This 
is because the amounts of trace gases in the atmospheres depend critically on the component 
gas sources (surface emission) and sinks (chemical reactions initiated by UV photolysis, as 
well as surface deposition). This is in contrast to giant exoplanets, where there is no surface 
for emission, and where UV photolysis leads largely to perturbations on atmospheric com- 
position, not to the definition of atmospheric composition. For super Earths and smaller 
terrestrial exoplanets, the key processes to consider in determining atmospheric composition 
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are photolysis, kinetics of the reactions between atmospheric components, vertical diffusion 
of molecules, atmospheric escape, dry and wet deposition, and condensation and sedimenta- 
tion of condensable species. An atmospheric spectrum may have strong features from spec- 
troscopically active trace gases whose lifetime depends on photochemistry; some of these 
trace gases may be hallmarks for specific atmospheric scenarios. Observational signatures of 
non-equilibrium chemistry have been observed for several transiting hot J up iters and "hot 
Neptunes" (e.g., Swain et al. 2008; Stevenson et al. 2010; Madhusudhan & Seager 2011). 

Theoretical studies of photochemistry on exoplanets have also been ongoing. One- 
dimensional photochemistry models have long been used to understand the atmospheres of 
planets and moons in the Solar System (see Yung & Demore 1999 and references therein), and 
such models and their networks of chemical reactions have been the foundation for exoplanet 
models. Photochemistry models have been very successful in understanding the atmosphere 
of Earth (e.g., Kasting et al. 1985; Zahnle 1986; Seinfeld & Pandis 2006 and references 
therein), Mars (e.g. Yung & Demore 1999; Zahnle et al. 2008), Venus (e.g. Krasnopol- 
sky & Pollack 1994; Zhang et al. 2012), and Titan (e.g. Atreya et al. 2006). After the 
discoveries of transiting hot Jupiters, photochemistry models were adapted and applied to 
exoplanets. Liang et al. (2003) suggest that, driven by the photolysis of water, hot Jupiters 
have much higher H concentrations than Jupiter in their upper atmospheres. Zahnle et al. 
(2009a) and Line et al. (2010) investigated the hydrocarbon photochemistry in hot Jupiters' 
H 2 -dominated atmospheres with a complete chemical network (as compared to the work of 
Liang et al. (2003) that was based on a simplified hydrocarbon chemistry scheme). Zahnle et 
al. (2009b) studied sulfur photochemistry in hot Jupiters' atmospheres, and suggested that 
H 2 S photolysis becomes important at altitudes above the ~100 Pa pressure level. Moses et 
al. (2011) studied atmospheric chemistry of HD 209458b and HD 189733b, the most well- 
characterized hot Jupiters, with a complete chemical network involving C, H, O, N species 
and derived observational signatures of disequilibrium processes in their atmospheres. Kop- 
parapu et al. (2012) presented a photochemistry model of a possibly carbon-rich hot Jupiter, 
WASP-12b, and suggested that C 2 H 2 and HCN, produced by photochemistry, may contribute 
significantly to the opacity of its atmosphere. More recently, Millier-Ricci Kempton et al. 
(2012) presented the first analysis of photochemistry in a H 2 -dominated atmosphere of mini 
Neptune GJ 1214b. 

A key challenge in developing a photochemistry model for terrestrial exoplanet atmo- 
spheres is handling the very wide range of plausible atmospheric compositions. The broad 
possibilities for super Earth atmospheric composition arise from the ideas for atmospheric 
origin. Super Earth atmospheres may originate from the capture of nebular gases, degassing 
during accretion, and/or degassing from subsequent tectonic activity (Elkins-Tanton & Sea- 
ger 2008; Schaefer & Fegley 2010). In principle, the atmospheres of super-Earths and ter- 
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restrial exoplanets may be reducing (H 2 -rich), oxidized (N 2 -rich, or C0 2 -rich), or oxidizing 
(0 2 -rich). It is also possible that some super-Earth atmospheres are water rich (Rogers & 
Seager 2010; Bean et al. 2010; Miller- Ricci Kempton et al. 2012). To date, photochemistry 
models are usually specific to certain types of atmospheres, because very different photo- 
chemical reactions dominate the chemistry of those atmospheres (for example in oxidized vs. 
reducing atmospheres). In particular, photochemistry of hydrogen-dominated atmospheres 
on terrestrial planets is yet to be explored; but hydrogen-dominated atmospheres may still 
create habitable surface temperature through collision-induced absorption (Pierrehumbert 
k Gaidos 2011; Wordsworth 2012). 

A second significant challenge for a terrestrial exoplanet photochemistry model is the 
range of free parameters that describe key physical processes, namely vertical diffusion, fluxes 
of surface emission, and rates of dry and wet deposition. While the fundamental equations 
of continuity that govern the atmospheric composition are the same, these coefficients must 
be kept as free parameters in order to explore the range of atmospheric mixing and the 
range of plausible surface conditions, which yield the wide variety of possible exoplanets' 
atmospheres. The motivation of this work is to provide a photochemistry model that can be 
applied to study atmospheres of a wide range of redox states and that can handle the range 
of free parameters for physical processes in terrestrial exoplanet atmospheres. 

In this paper, we present a photochemistry model developed from the ground up from 
basic chemical and physical principles and using both established and improved computer 
algorithms, which have the capacity of modeling very different terrestrial exoplanet atmo- 
spheres in a consistent way. The design of the photochemistry model aims at providing the 
maximum flexibility in treating atmospheres of very different compositions, especially oxi- 
dation states. Besides stellar parameters and planetary parameters (mass, radius, etc.), the 
photochemistry model accepts a series of input parameters including temperature structure 
of the atmosphere, dominant species and mean molecular mass of the atmosphere, coeffi- 
cients of vertical diffusion, fluxes of surface emission and fluxes of dry and wet depositions. 
These input parameters are treated as free parameters so that the photochemistry model 
can simulate a wide range of atmospheric compositions. 

In § 2 we describe our photochemistry model in detail. In § 3 we validate our photo- 
chemistry model by simulating the atmospheres of current Earth and Mars. In § 4 we present 
three benchmark cases for exoplanets, representing reducing, weakly oxidizing, and highly 
oxidizing atmospheres. We discuss key parameters that determine atmospheric compositions 
in § 4 and conclude in § 5. 
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2. A Photochemistry Model 

The purpose of the photochemical model is to compute the steady-state chemical com- 
position of an exoplanetary atmosphere. The system can be described by a set of time- 
dependent continuity equations, one equation for each species at each altitude. Each equa- 
tion describes: chemical production; chemical loss; diffusion (contributing to production 
or loss); sedimentation (for aerosols only); and specified gains and losses on the lower and 
upper boundaries of the atmosphere. The lower boundary has assigned geological or bio- 
logical source fluxes and assigned deposition rates of species and the upper boundary has 
diffusion-limited atmospheric escape, relevant for light species. 

With the assumption of a one-dimensional plane-parallel atmosphere, the unknowns to 
be computed in the model are the number densities of each species at each altitude. Starting 
from an initial state, the system is evolved to the steady state in which production and loss are 
balanced for each species at each altitude. Because the removal timescales of different species 
are very different, the inverse Euler method is employed for the numerical time stepping. 
We allow the time step to be adjusted according to how fast the atmospheric components 
change (i.e., the variation timescale), and determine if the solution converges to the steady- 
state solution by evaluating the variation timescale for each species at each altitude, as well 
as the global production and loss balance for each species. Once the steady-state solution is 
found, we use a separate code to compute the exoplanet atmosphere spectrum for thermal 
emission, reflected stellar radiation, and transmission of stellar radiation. 

We developed the photochemistry model from the ground up, and tested and validated 
the model by reproducing the atmospheric composition of Earth and Mars, in terms of 
globally- averaged ID vertical mixing ratio profiles. For one of the tests, we started with a 
temperature-pressure profile of Earth's atmosphere, and with a subset of Earth's composition 
at present-day values, specifically N 2 , 2 , and C0 2 with assigned mixing ratios. We included 
surface emission of major trace gases, assumed to this exercise to be emitted according to 
typical globally averaged emission rates, for CO, CH 4 , NH 3 , N 2 0, NO, N0 2 , S0 2 , OCS, 
H 2 S, and H 2 S0 4 . The code correctly predicts the amounts of these trace gases in Earth's 
troposphere, and the amounts of ozone in the atmosphere as well as the location of the ozone 
layer. The code also produces vertical profiles of active hydrogen species (OH and H0 2 ) and 
active nitrogen species (NO, N0 2 , and HN0 3 ) in the stratosphere that are consistent with 
observations. As another test, we reproduce the chemical composition of the current Mars 
atmosphere, in agreement with measured composition (e.g., Krasnopolsky 2006) and previous 
1-D photochemistry model results (Nair et al. 1994; Yung & DeMore 1999; Zahnle et al. 
2008). In particular, the code correctly illustrate the effect of HO x catalytic chemistry that 
stabilizes Mars' C0 2 atmosphere. The results of these tests are described in more detail 
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below (Section 3). 

The photochemistry model is designed for exoplanet applications, and has features that 
are not needed for studies of individual Solar System objects but which are useful for study 
of potential exoplanet atmospheres. The most important feature that yields the capacity of 
treating reducing and oxidized atmospheres with the same code is the flexibility of choosing 
chemical species to be considered in photochemical equilibrium. The species that have life- 
times less than the numerical time step can be treated as being in photochemical equilibrium. 
Our photochemistry model features a "burn-in" period in which all species' number densities 
are rigorously computed and then those species that satisfy the lifetime condition are treated 
in equilibrium so that the computation can be optimized for speed. We have coded the Ja- 
cobian matrix for the implicit Euler method analytically, which improves the numerical rigor 
of the code. In this way our photochemistry model can quickly find the steady-state solution 
of atmospheric composition starting from any sensible initial conditions for both reducing 
and oxidized atmospheres. We emphasize that the ability to compute atmospheric compo- 
sition of different oxidation states and the elimination of the need of providing fine-tuned 
initial conditions are essential for exoplanet exploration, because the oxidation states are 
unknown for terrestrial exoplanets and there are no benchmark cases to provide common 
initial conditions. 

Another important feature of our photochemistry model is that it treats a wide range 
of free parameters for terrestrial exoplanets. We do not hardwire any model parameters 
that meant to be free for exoplanet exploration, including those parameters for planets, such 
as stellar spectrum and surface gravity, those parameters for atmospheres, such as mean 
molecular mass, refractive indices, temperature profiles (which can also be self-consistently 
computed), and eddy diffusion coefficients, and those parameters that are specific to partic- 
ular molecules, such as solubility, rainout rate, saturation vapor pressure, surface emission 
rate, and dry deposition velocity. The code thus has a clear structure that makes extensive 
parameter exploration possible. 

The third important feature of our photochemistry model is that we offer the flexibility 
of choosing a subset of chemical species and chemical reactions for the computation. It is in 
particular important for exoplanet exploration to be able to isolate a chemical system from 
the complex network of atmospheric chemistry; and our photochemistry model responds to 
this need. We design our code to be able to include or exclude the effect of any species 
or reactions in the complex chemical network, in order to understand the fundamentals of 
atmospheric chemistry on terrestrial exoplanets. Also, we pay attention to the temperature 
range in which chemical kinetic rates and photochemical cross sections are valid. The tem- 
perature range of exoplanet atmospheres can be significantly larger than the Solar System 
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planet atmospheres. We label the chemical reactions whose rates are measured only at low 
temperatures (T < 400 K) and the chemical reactions whose rates are measured only at 
high temperatures (T > 1000 K), and use them only at appropriate temperature ranges. 
In addition, we consider the temperature dependence of photochemical cross sections in the 
model. Finally, we include a basic aerosol formation scheme in the model, the computation 
of which only relies on the saturation pressure of relevant substances and a timescale of 
condensation and sublimation. The aerosol scheme is simple, but could easily be expanded 
to treat non-conventional aerosols in the atmosphere. 

We now present the detailed formulation of the photochemistry model. 



2.1. Fundamental Equations 

The coupled 1-dimensional continuity-transport equation that governs the chemical com- 
position is 

d S L = P-nL- d * (1) 

dt dz' K J 

where n is the number density of a certain species (cm -3 ), z is the altitude, P is the pro- 
duction rate of the species (cm -3 s _1 ), L is the loss rate of the species (s _1 ), and $ is the 
upward vertical transport flux of the species (cm -2 s -1 ). The flux can be approximated by 
the eddy diffusion together with the molecular diffusion, viz., 

where K is the eddy diffusion coefficient (cm 2 s _1 ), D is the molecular diffusion coefficient 
(cm 2 s _1 ), N is the total number density of the atmosphere, / = n/N is the mixing ratio 
of the species, Hq is the mean scale height, H is the molecular scale height, T is the tem- 
perature (K), and is the thermal diffusion factor (Banks & Kockarts 1973; Levine 1985; 
Krasnopolsky 1993; Yung & DeMore 1999; Bauer & Lammer 2004; Zahnle et al. 2006). 
The first term of equation fl2]) represents eddy diffusion, and the last two terms represent 
molecular diffusion. The vertical transport flux is written in terms of the derivatives of the 
mixing ratio rather than in terms of the number density because this form is simpler and 
more straightforward to be implemented in numerical schemes. We show that equation (j2J) 
is equivalent to the standard vertical diffusion equation in Appendix [A] 

The molecular diffusion coefficient D and the thermal diffusion factor ax can be deter- 
mined from the gas kinetic theory, but the eddy diffusion coefficient is a more speculative 
parameter and must be estimated empirically. For molecular diffusion, we use the following 
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expressions in cm 2 s 1 for H and H 2 in a N 2 -based atmosphere 

4.87 x 10 17 x T - 698 



D(R, N 



2) 



N 

2.80 x 10 17 x T - 740 



^(H 2 ,N 2 ) = ; (3) 



and in a C0 2 -based atmosphere 

L>(H,CO 



3.87 x 10 17 x T - 750 



N 

nm _,_ . 2.15 x 10 17 x T - 750 (A . 
D(H 2 ,C0 2 ) = ; (4) 



N 



and in a H 2 -based atmosphere 



. 8.16 x 10 17 x T - 728 
D(R, H 2 ) = . (5) 

In equations ([MS]), T has a unit of K and N has a unit of cm -3 . The functional form 
of equations ([MSD is derived from the gas kinetic theory and the coefficients are obtained 
by fitting to experimental data (Marreo & Mason 1972; Banks & Kockarts 1973). The 
expressions for molecular diffusion coefficients are valid for a wide range of temperatures, 
up to 2000 K, except for -D(H 2 , C0 2 ) that is only valid at low temperatures, i.e., T < 550 
K (Marreo & Mason 1972). -D(H, C0 2 ) is assumed to be 1.8 times larger than D(H 2 ,C0 2 ) 
(Zahnle et al. 2008). The thermal diffusion factor is taken as a constant for H and H 2 as 
qt = —0.38 (Banks & Kockarts 1973). The negative sign of ax corresponds to the fact that 
the lightest molecules tend to diffuse to the warmest region. 

The eddy diffusion coefficient is the major uncertainty in the 1-D photochemistry model. 
For Earth's atmosphere, the eddy diffusion coefficient can be derived from the number den- 
sity profile of trace gases (e.g. Massie & Hunten 1981). The eddy diffusion coefficient of 
Earth's atmosphere is characterized by the convective troposphere and the non-convective 
stratosphere, which is not necessarily applicable for the exoplanets. The eddy diffusion 
coefficient may be typically parameterized as 

K = K T (z < z T ), 

(NM\ 1/2 



K = minlK H , Kt ( j\r J ) ( z > z t), (6) 

where Kt, Kh and zt are independent parameters satisfying Kh > Kt- This formula is 
adapted from Yung & DeMore (1999) with the requirement of continuity in K. Our code 
can either import eddy diffusion coefficient from a file or specify eddy diffusion coefficient 
according to equation (jSJ). 
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The goal of the photochemistry model is to obtain a steady-state solution of each species, 
or a set of n(z), which makes the left hand side of the equation ([1]) vanish and satisfies the 
boundary conditions. Assuming N x species in the model and N equally stratified layers of 
the atmosphere, we transform the continuity equation ([T]) into a discrete form as 



at 



$i+l/2 — $j-l/2 

Az ' 



(7) 



where the subscript i denotes physical quantities in the zth layer, and the subscripts i + 1/2 
and % — 1/2 mean that the flux is defined at the upper and lower boundary of each layer. It is 
physically correct to define the flux term at the boundary of each layer, which is also crucial 
for numerically preserving hydrostatic equilibrium for an atmospheric transport scheme that 
uses number density as independent variables. According to equation (J21), 
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where m a is the mean molecular mass of the atmosphere, m is the molecular mass of the 
species, g is the gravitational acceleration, and &b is the Bolzmann constant. We have 
approximated / i+1 / 2 by (fi+i + /»)/2 in equation (J2J). 

The combination of equation ([7]) and (jHJ) gives the following 2nd ordered centered discrete 
differential equation to be solved numerically, 
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where 



a: 



c i+l/2 
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In principle, dn/dt = is equivalent to a set of N x Ni nonlinear algebraic equations. The 
set of nonlinear algebraic equations may be solved numerically by Newton-Raphson methods 
(Press et al. 1992). However, in practice, we find that Newton-Raphson methods require 
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an initial guess to be in the vicinity of the solution. Instead, we treat the problem as a 
time-stepping problem by evolving the system according to equation (jUJ) to the steady state. 

Implicit numerical methods are implemented to solve equation OH]). Due to the orders- 
of-magnitude differences in chemical loss timescales, the system is numerically stiff. We use 
the inverse-Euler method for the time stepping, as in most previous photochemistry models 
(e.g., Kasting et al. 1985; Nair et al. 1994; Zahnle et al. 2008). For each time step, we need 
to invert a matrix of dimension N x Ni. As seen in equation ([9]), the variation in any layer only 
depends on the number density in that layer and adjacent layers. As a result, the matrix 
is by nature block tridiagonal with a block dimension of N x , which is solved efficiently by 
the Thomas algorithm (Conte & DeBoor 1972). The time step is self-adjusted in a way that 
the code updates the time step after each iteration according to the variation timescale of 
the whole chemical-transport system, i.e., the minimum variation timescale of each species 
at each altitude. As the system converges, larger and larger time steps are chosen. 

The criteria of convergence to the steady-state solution is that the variation timescale 
of each species at each altitude is larger than the diffusion timescale of whole atmosphere, 
and the fluxes of gain and loss balance out for all species. The gain fluxes include surface 
emission, chemical production and condensation; whereas the loss fluxes include chemical 
loss, dry deposition, wet deposition, atmospheric escape, and condensation. We require that 
the ratio between the net global flux and the column-integrated number density for all long- 
lived species is small compared to the diffusion timescale. This condition is very important 
in the determination of convergence. If any possible long-term trends of major species can 
be detected, we run the code for an extended time period to test the convergence rigorously. 

The current photochemistry model can compute concentrations of 111 molecules or 
aerosols made of C, H, O, N, S elements. These species are C0 2 , H 2 , 2 , N 2 , H 2 0, O, 
0( 1 D), 3 , H, OH, H0 2 , H 2 2 , N, NH, NH 2 , NH 3 , N 2 H 2 , N 2 H 3 , N 2 H 4 , N 2 0, NO, N0 2 , 
N0 3 , N 2 5 , HNO, HN0 2 , HN0 3 , HN0 4 , C, CO, CH 4 , CH, CH 2 , GE\, CH 3 , CH 2 0, CHO, 
CH 3 0, CH 3 2 , CH0 2 , CH 2 2 , CH 4 0, CH 4 2 , C 2 , C 2 H, C 2 H 2 , C 2 H 3 , C 2 H 4 , C 2 H 5 , C 2 H 6 , 
C 2 HO, C 2 H 2 0, C 2 H 3 0, C 2 H 4 0, C 2 H 5 0, HON, CN, CNO, HCNO, S, S 2 , S 3 , S 4 , S 8 , SO, S0 2 , 
SO*, SO^, S0 3 , H 2 S, HS, HSO, HS0 2 , HS0 3 , H 2 S0 4 , OCS, CS, CH 3 S, CH 4 S, CH 3 N0 2 , 
CH 3 ON0 2 , CH 5 N, C 2 H 2 N, C 2 H 5 N, C 3 H 2 , C 3 H 3 , CH 3 C 2 H, CH 2 CCH 2 , C 3 H 5 , C 3 H 6 , C 3 H 7 , 
^Hg, C 4 H, C 4 H 2 , C 4 H 3 , C 4 H 4 , C 4 H 5 , 1-C 4 H 6 , 1,2-C 4 H 6 , 1,3-C 4 H 6 , C 4 H 8 , C 4 H 9 , C 4 H 10 , 
C 6 H, C 6 H 2 , C 6 H 3 , C 6 H 6 , C 8 H 2 , H 2 S0 4 aerosols, S 8 aerosols and organic hazes. 

An important feature of the photochemistry code is the flexibility to choose a subset of 
molecules and aerosols to study a particular problem. Since the computation time for each 
time step scales with N^, decreasing the number of chemical species in the main loop greatly 
reduces the computation time. As is common practice in the effort to reduce the stiffness of 
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the system and improve the numerical stability, "fast" species with relatively short chemical 
loss timescales are computed directly from the diagnostic equation, namely 

n=-, (10) 

which implies that photochemical equilibrium can be achieved within each time step. The 
choice of fast species varies with the atmospheric composition and should be considered 
on a case-by-case basis. The choices of "fast" species were usually hard-wired in previous 
photochemistry codes, which fundamentally limited their application to certain specific types 
of atmospheres (e.g., to Earth and to Mars). Our photochemistry code offers the flexibility 
to adjust the choice of species in photochemical equilibrium, which yields the capacity of 
treating atmospheres having very different oxidation states with the same code. We evaluate 
P and L in equation fTiTJj) at each time step using the number density of all species from 
previous time step and determine the number density of fast species according to equation 
( fit)]) after the current time step linearly. We neglect nonlinearity due to multiple species 
in photochemical equilibrium; as a result, we do not allow strongly inter-dependent species 
to be concurrently considered as in photochemical equilibrium, for example for S 2 , S 3 , and 
S 4 . Our approach yields stable convergence to the steady-state solution from virtually any 
initial test solutions. We always verify the mass balance of the steady-state solution. 



2.2. Chemical Kinetics 

We compiled a comprehensive list of chemical and photochemical reactions from the 
literature. The production and loss rates in equation ([7]) are provided from all chemical 
and photochemical reactions that produce or consume the relevant molecule. In the generic 
model, we included 645 bimolecular reactions, 85 termolecular reactions, and 93 thermal 
dissociation reactions. We included the thermal dissociation reactions to make the photo- 
chemistry model potentially adaptable to simulate hot planets. We mainly used the updated 
reaction rates from the online NIST database (http:/ /kinetics. nist.gov) and refer to the JPL 



publication (Sander et al. 2011) for the recommended rate when multiple measurements are 
presented in the NIST database. For termolecular reactions, we used the complete formula 
suggested by the JPL publication (Sander et al. 2011) when the low-pressure-limiting rate 
and the high-pressure-limiting rate are available. For all reaction rates we use the Arrhenius 
formula to account for the dependence on temperature when the activation energy data are 
available, otherwise we adopted the value of experimental measurements, which are usually 
under ~298 K. We performed careful comparisons of our reaction lists and corresponding 
reaction rates with those used by Nair et al. (1994), Pavlov & Kasting (2002), and Zahnle et 
al. (2008). Most of our reaction rates are the same as those used by other codes. For a dozen 
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reactions, we find updated reaction rates listed in either the JPL publication or the NIST 
database updated after the publication of the cited photochemistry models. Also, for the 
general purpose of our photochemistry model, we include the reactions that are very slow at 
low temperatures (i.e., 200 - 400 K) but may be become important at higher temperatures 
> 1000 K. In addition, a number of reactions that lack laboratory-measured rates may be 
important for low-temperature (i.e., 200 - 400 K) applications. We have adopted the rates 
for unmeasured sulfur reactions from Turco et al. (1982), Kasting (1990), and Moses et al. 
(2002), and the rates for reactions of C>2 hydrocarbons from Yung & DeMore (1999). All 
chemical reactions and their reaction rates are tabulated in Table [TUl 

Chemical reaction rates may sensitively depend on temperature. The reported reac- 
tion rates are usually valid in certain specific temperature ranges. For example, reaction 
rates recommended by the JPL publication (Sander et al. 2011) are the reactions relevant 
to Earth's atmosphere, in general valid within the small range of 200 - 300 K. Hence it 
may be problematic to extend the rate expressions to high-temperature cases, such as hot 
planets. At the high-temperature end, the valid temperature range for some of the reaction 
rates measured in the combustion chemistry is a few thousand K. These reactions may be 
important for modeling hot Jupiters (e.g., Moses et al. 2011), but they should not be in- 
cluded in the computation of low-temperature cases, e.g., for Earth-like planets. We annotate 
the low-temperature reactions and the high-temperature reactions in Table [10J Moreover, 
most chemical reactions involving free radicals have very small activation energy, and thus 
weak temperature dependencies. Usually the reactions rate of free radicals are measured at 
room temperature and can be used in photochemistry models under different temperatures 
(Sander et al. 2011). Common atmospheric free radicals are: 0( 1 D), OH, H, N, NH 2 , N 2 H 3 , 
C, CH, CH 2 , CH 2 , CH 3 , CHO, CH0 2 , CH 3 0, CH 3 2 , C 2 , C 2 H, C 2 H 3 , C 2 H 5 , C 2 HO, C 2 H 3 0, 
C 2 H 5 0, CN, CNO, S, S 2 , S0 3 , HS, HSO, HS0 2 , HS0 3 , CH 3 S and CS. For reactions with 
these molecules, if no activation energy is reported in the literature, we assume that the 
activation energy is negligible and their kinetic rates do not depend on temperature. 

2.3. Interaction with Radiation 

Molecules in the atmosphere absorb ultraviolet (UV) and visible light from the host 
star. If the absorbed photon carries enough energy, the molecule may be photodissociated 
to form free radicals. The photodissociation rate is proportional to the number density of 
photons with UV and visible wavelengths at each altitude. For direct stellar radiation, the 
optical depth r includes the contribution from molecular absorption r a , Rayleigh scattering 
r r , and aerosol particle extinction r m . For the multiple-scattered (diffusive) radiation, we 
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use the 5-Eddington 2-stream method implemented based on the formulation of Toon et al. 
(1989). The actinic flux of the diffusive radiation is F^s = 2(F + + F~) where F + and F~ 
are the diffusive flux in the upward and downward direction. The photolysis flux at a certain 
altitude includes both the direct radiation and the diffusive radiation, i.e., 



where F is the radiation flux at the top of the atmosphere where r = 0, and Ho is the angle 
of the path of sunlight. By default we assume the zenith angle of the star to be 57.3° (see 
Appendix IB1 for justification). The photodissociation rate is then 



where a a is the absorption cross section, g(A) is the quantum yield that is defined as the 
ratio between the yield of certain photodissociation products and the number of photons 
absorbed, and L(X , z) is the actinic flux with units of quanta cm -2 s _1 nm _] 0. The 1/2 
factor is included to account for diurnal variation (e.g. Zahnle et al. 2008), which is not 
used to model the dayside of a tidally lock planet. With the general-purpose model, we treat 
70 photodissociation reactions. Different branches resulting from the photodissociation of 
one molecule are treated as different photodissociation reactions. The photodissociation 
reactions, the sources of data for cross sections and quantum yields, and the rates on the 
top of Earth's atmosphere are tabulated in Table [71 

For the UV and visible cross sections and the quantum yields, we use the recommended 
values from the JPL publication (Sander et al. 2011) when available. We also use the cross 
sections from the MPI-Mainz-UV-VIS Spectral Atlas of Gaseous Molecules^] when the JPL 
recommended values are not available or incomplete. There are a number of molecules of 
atmospheric importance that lack UV and visible cross sections and quantum yields, and 
we have estimated photolysis rates for them in the following way. The photodissociation 
timescale of S 2 has been measured to be ~ 250 s at 1 AU from the Sun (DeAlmeida & 
Singh 1986), from which we estimate the photolysis rate of S 2 . For other molecules that 
have no UV or visible-wavelength cross sections available in the literature, we assume their 



lr The actinic flux is the quantity of light available to molecules at a particular point in the atmosphere 
and which, on absorption, drives photochemical processes in the atmosphere. It is calculated by integrating 
the spectral radiance over all directions of incidence of the light. The actinic flux is distinguished from the 
spectral irradiance in the way that it does not refer to any specific orientation because molecules are oriented 
randomly in the atmosphere. 

2 www. atmosphere . mpg . de / spectral- atlas- mainz 
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photolysis rate to be the same as another molecule that has similar structure. For example, 
the photolysis rate of HNO is assumed to be the same as HN0 2 (Zahnle et al. 2008); the 
photolysis rate of HSO is assumed to be the same as H0 2 (Pavlov & Kasting 2002); and the 
photolysis rate of N 2 H 2 is assumed to be the same as N 2 H 4 . 

Temperature dependencies of photolysis cross sections and quantum yields are consid- 
ered. Notably, at 200 K, compared to room temperature, N 2 0, N 2 5 , HN0 3 , OCS, C0 2 
have smaller UV cross sections, leading to photolysis rates more than 10% lower; whereas 
N0 3 has a larger UV cross section, leading to a photolysis rate more than 10% higher. S0 2 
has complex band structures in its UV spectrum and the cross sections depend on temper- 
ature as well. We take into account any other temperature dependencies reported for a gas 
at temperatures outside 290 - 300K (see Table [7] for notes on temperature dependencies). In 
most cases, cross sections are also measured at lower temperature such as 200 K, primarily 
for Earth investigations (Sander et al. 2011). We use linear interpolation to simulate cross 
sections between 200 and 300 K, and do not extrapolate beyond this temperature range. It 
is worthwhile noting that UV cross sections of almost all gases at temperatures significantly 
higher than room temperature are unknown, and temperature dependencies of cross sections 
are not considered in previous high-temperature photochemistry models of hot exoplanets 
(e.g., Moses et al. 2011; Miller-Ricci Kempton et al. 2012). This may be a plausible simpli- 
fication for gases without significant band structure in UV, but for gases such as H 2 , C0 2 
and S0 2 , whose band structures are sensitive to temperature, extrapolation of cross sections 
to high temperature might induce significant errors. One needs to be cautious about the 
uncertainty of photolysis rate at temperatures much higher than room temperature. 

Rayleigh scattering from atmosphere molecules introduces additional optical depth, par- 
ticularly important for attenuation of UV radiation. The optical depth due to Rayleigh 
scattering is 



where m r is the real part of the refractive index of the molecule, C r is a corrective factor to 
account for the anisotropy of the molecule, and N s is the number density at the standard 
condition (1 atm, 273.15 K). The refractive index depends on the main constituent in the 
atmosphere. The refractive index of Earth's atmosphere is from Seinfeld & Pandis (2006); 
the refractive index of H 2 is from Dalgarno & Williams (1962); the refractive index of N 2 is 
from Cox (2000); the refractive index of C0 2 is from Old et al. (1971); and refractive indices 
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in which the Rayleigh scattering cross section is (e.g., Liou 2002) 




s 



(14) 



- 15 - 



of CO, and CH 4 are given in Sneep & Ubachs (2005). In principle the correction factor C r 
depends on the molecule and the wavelength; but C r is usually within a few percent with 
respect to the unity, except for C r ~ 1.14 for C0 2 (Sneep & Ubachs 2005). In the following 
we assume C r = 1.061, the value for Earth's atmosphere at ~ 200 nm (Liou 2002). 

For the spectrum of a solar-type star (G2V), we use the Air Mass Zero (AM0) reference 
spectrum produced by the American Society for Testing and Material^. The AM0 spectrum 
covers a wavelength range from 119.5 nm to 10 microns. For the extreme-UV spectrum, we 
use the average quiet-Sun emission provided by Curdt et al. (2004). 



Microphysical processes involved in the formation of atmospheric aerosols are nucleation, 
condensational growth, and coagulation (Toon & Farlow 1981; Seinfeld & Pandis 2006). It 
is beyond the purpose of this work to simulate these microphysical processes in detail. For 
photochemically produced aerosols, the competition between coagulation and sedimentation 
mainly determines the particle size distribution. A complete treatment of the atmospheric 
aerosols involves solving the steady-state size distribution function (e.g. Seinfeld & Pandis 



In our photochemical model, we simplify the problem by assuming the particle radius to 
be a free parameter. In Earth's atmosphere, aerosols formed in the atmosphere often have a 
lognormal size distribution around 0.1-1 /im (Seinfeld & Pandis 2006). We define a particle 
radius parameter, r p , to be the surface area average radius. This parameterization allows us 
to separate the complexity of aerosol formation from the photochemistry model, as well as 
to explore how the particle size affects the overall chemical composition. 

We compute the production and loss rate of a molecule in the condensed phase (i.e., 
aerosols) based on its condensation timescale. When the molecule becomes supersatu- 
rated at a certain altitude, condensation can happen and aerosols form. The condensa- 
tion/evaporation timescale is given by Hamill et al. (1977) and Toon & Farlow (1981) as 



where t c is the condensation/evaporation timescale, m is the mass of molecule, p p is the 
particle density, is the Boltzmann constant, T is the atmospheric temperature, n g is the 
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Treatment of Aerosols 



2006). 
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3 http: / /rredc. nrel.gov / solar /spectra/amO / 
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number density of the corresponding gas, and n v is the saturated vapor number density at the 
corresponding pressure. The formula ffl5|) is suitable for both condensation and evaporation 
processes. When n g > n v , the gas phase is saturated, so the condensation happens and 
t c > 0. When n g < n v , the gas phase is unsaturated, so the evaporation is possible and 
t c < 0. The production or loss rate of the molecule in the condensed phase is 

P = ^, L = -. (16) 

We include gravitational settling in the mass flux term of the continuity-transport equa- 
tion (equation [TJ) for aerosol particles in addition to eddy mixing. The additional gravita- 
tional downward flux of the aerosol particle is 

$ F = -v F n C} (17) 

where v-p is the settling velocity of the particle in the atmosphere. The settling velocity is 
reached when the gravitational force is balanced by the gas drag. For aerosols with diameter 
of order of 1 /an, the settling velocity is reached within 10~ 5 s in Earth's atmosphere (Seinfeld 
& Pandis 2006). Therefore we assume the falling velocity to be the settling velocity. The 
settling velocity can be derived from the Stokes' law (Seinfeld & Pandis, 2006) as 

2 rJppgCc 
9 fi 

where g is the gravitational acceleration, /i is the viscosity of the atmosphere, and C c is the 
slip correction factor related to the mean free path (A) of the atmosphere as 



C c = l + - 



1.257 + 0.4 exp 



l.lr p 



(19) 



The treatment of photochemically produced aerosols described above is applicable to 
any molecules that could reach saturation as a result of photochemical production. In Earth's 
atmosphere, the photochemical aerosols include sulfate aerosols (H 2 S0 4 ), sulfur aerosols (S 8 ), 
organic hazes, nitric acid aerosols (HN0 3 ), and hydrochloric acid aerosols (HC1). For now 
we have implemented sulfate aerosols (H 2 S0 4 ) and sulfur aerosols (S 8 ). These aerosols and 
water vapor are the common condensable materials at habitable temperatures that commonly 
exist in planetary atmospheres (e.g., Kasting et al. 1989; Pavlov et al. 2000; Seinfeld & 
Pandis 2006). The required data for including an aerosol species in the photochemistry 
model is the saturation vapor pressure. Saturation vapor pressure of H 2 S0 4 is taken as 
recommended by Seinfeld & Pandis (2006) for atmospheric modeling, with a validity range 
of 150 - 360 K. Saturation pressure of S 8 is taken as the total sulfur saturation pressure 
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against liquid sulfur at T > 392 K and solid (monoclinic) sulfur at T < 392 K tabulated and 
expressed by Lyons (2008). In addition to aerosols, we use equation f Ti5]fT6]) to compute the 
process of condensation of water vapor in the atmosphere. We do not consider evaporation 
of condensed water in the atmosphere because water droplets may grow by aggregation and 
rapidly precipitate out (Seinfeld & Pandis 2006). Saturation pressure of water is taken as 
that against ice at temperature lower than 273.16 K (Murphy & Koop 2005), and that 
against liquid water at temperature higher than 273.16 K (Seinfeld & Pandis 2006). 

We treat the optical effect of aerosols by considering both scattering and absorption. 
Assuming a homogeneous sphere, the Mie theory (see Van de Hulst, 1981 for a detail de- 
scription) computes extinction cross section (cr ext ), single scattering albedo (w B ) and asym- 
metric factor ((7asym); based on the following parameters: the refractive index of the material 
(m r + rriii), r p , and the wavelength. In this paper, we use the refractive index of S 8 aerosols 
from Tian et al. (2010) for the UV and visible wavelengths and from Sasson et al. (1985) 
for infrared (IR) wavelengths. We use the refractive index of H 2 S0 4 aerosols (assumed to 
be the same as 75% sulfuric acid solution) from Palmer & William (1975) for UV to IR 
wavelengths, and Jones (1976) for far IR wavelengths. 



2.5. Boundary Conditions 

The atmospheric chemical composition of a terrestrial exoplanet is ultimately deter- 
mined by boundary conditions. The upper boundary conditions describe the atmospheric 
escape. The lower boundary conditions describe the surface emission, the deposition of 
molecules and aerosols to the surface, or the presence of a large surface reservoir of certain 
molecule (e.g. H 2 0). The boundary conditions need to be properly provided to capture the 
physics of an exoplanet atmosphere. 

The upper boundary conditions describe the atmospheric escape of an terrestrial exo- 
planet. The escape rates of exoplanet atmospheres are fairly uncertain depending on stellar 
soft X-ray and UV luminosity, exosphere chemistry, existence of magnitude fields, etc (e.g., 
Yelle et al. 2008; Tian 2009). We thus provide the following options of specifying escape 
rates in the photochemistry code: 

Type 1: &N t +i/2 = 0, or no escape; 

Type 2: <3>at ;+ i/ 2 = n^Vnm, where Vu m is the diffusion-limited escape velocity; 
Type 3: &Ni+i/2 is a assigned nonzero value. 
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Here we use the same notation for flux as in equation ([7HHJ), such that the upper boundary 
condition replaces the flux at the upper boundary of the layer Ni in equation (jSJ) . For the 
type 2 upper boundary condition (atmospheric escape), the diffusion-limited velocity (VJ im ) 
is 



where D Nl+ u 2 is the molecular diffusion coefficient evaluated at the top of the atmosphere. 
The diffusion-limited flux is the highest escape flux of an atmosphere in hydrostatic equi- 
librium (Hunten 1974). For Mars, the Jeans escape of H 2 reaches the diffusion-limited flux 
when the exobase temperature is above 400 K (Zahnle et al. 2008). We use the Type 2 upper 
boundary condition for the escape of H and H 2 in our model, and generally we assume no 
escape for all other species, i.e., the Type 1 upper boundary condition. The Type 3 upper 
boundary condition may be used when processes above the neutral atmosphere are impor- 
tant. For example, an influx of atomic N can represent the photodissociation of N 2 in the 
upper atmosphere. The Type 3 boundary condition may also be used when hydrodynamic 
escape has to be considered. 

The lower boundary conditions describe the interaction between the atmosphere and 
the surface, which includes surface emission and surface deposition. The three types of lower 
boundary conditions are: 

Type 1: n\ is assigned; 

Type 2: $i/ 2 = — niVoEP where Vdep is molecule-specific dry deposition velocity; 
Type 3: $i/ 2 is assigned. 

Again we use the same notation for flux as in equation flTJlE]), and replaces the flux at 
the lower boundary of the layer 1 in equation ([8]). 

The Type 1 lower boundary condition presents a large reservoir at the surface and for 
this paper we use this condition for water vapor to simulate the effect of a surface with 
oceans. This approach is equivalent to setting the relative humidity at the surface to be 
a constant. Note that specifying the Type 1 lower boundary condition means decreasing 
the number of free variables; n\ is fixed as the lower boundary condition and is no longer 
considered as a variable in the main computation loop. 

The Type 2 lower boundary condition specifies the dry deposition velocity, which is a 
key parameter that determines the chemical composition of the atmosphere and is a major 
unknown. The deposition velocity depends on both the dynamical properties of the lower 




(20) 



- 19 - 



atmosphere and the chemistry of the planet's surface. With the number density at the bottom 
layer (rij.) computed from the photochemistry model, the interaction between the bottom 
layer of atmosphere and the surface consists of two steps: first, the molecular transport 
across a thin stagnant layer of air adjacent to the surface, called the quasi-laminar sublayer 
; second, the uptake at the surface (Seinfeld & Pandis 2006). A parameterization of dry 
deposition velocity involving these two steps is described in Appendix 

Physically the dry deposition requires a sink at the surface; for a gas without effective 
surface sink the dry deposition velocity should be zero. In particular, the surface deposition 
of a number of gases, including CO, H 2 , CH 4 , and NH 3 , is primarily removed by microorgan- 
isms on Earth (e.g., Kharecha et al. 2005; Seinfeld & Pandis 2006; Kasting, 2012, private 
communication), and the canonical values for their dry deposition velocities are not applica- 
ble for presumably abiotic planets. Take carbon monoxide for an example. If a planet has no 
ocean, there is no known reaction that can consume CO at the surface, and therefore the dry 
deposition velocity of CO for an abiotic desiccated planet should be zero (like Mars). If the 
surface has an ocean, the dissolved CO may be naturally and slowly converted into acetates 
by OH in sea water, or biologically converted into acetates at a much faster rate. The CO 
deposition velocity has been estimated to be 1.2 x 10 4 cm s _1 for the most efficient biological 
removal and 10~ 8 ~ 10~ 9 cm s _1 on an abiotic ocean planet (Kasting 1990; Kharecha et 
al. 2005). For another example, it is common to use a fairly large dry deposition velocity 
for S0 2 (~ 1 cm s _1 ) to study the sulfur cycle in Earth's marine atmosphere (e.g., Toon 
et al. 1987). In contrast, the dry deposition velocity is assumed to be zero, or reduced by 
an arbitrary factor of up to 1000, to mimic a putative Mars ocean that is believed to be 
saturated with dissolved S0 2 and other sulfur species (e.g., Halevy et al. 2007; Tian et al. 
2010). In general, the dry deposition velocity depends on a broad context of planetary geo- 
chemistry, notably on the surface mineralogy, the acidity of ocean, and the surface pressure 
and temperature. Therefore it is critical for exoplanet exploration to understand the inter- 
action between atmospheric chemistry composition of terrestrial exoplanets and the surface 
deposition. 

In addition to the dry deposition at the surface, we also include the wet deposition 
throughout the atmosphere as a removal process for soluble species. We use the parameter- 
ization of Giorgi & Chameides (1985) as 

kR[Z) ~ tRX 55A V [L x 10-* + (H'RT(z))-iy { ] 

where fa is a reduction factor adjustable in the model, &h 2 o is the precipitation rate taken 
to be 2 x 10~ 6 s _1 , Ay is the Avogadro's number, L is the liquid water content taken to 
be 1 g m -3 in the convective layer near the surface, and H' is the effective Henry's Law 
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constant measuring the solubility of the molecule in the unit of mol dm -3 atm -1 . The 
effective Henry's law constant may differ from the standard Henry's law constant when 
taking into account dissociation in the aqueous phase (Giorgi & Chameides 1985; Seinfeld & 
Pandis 2006). In the model we use the effective Henry's Law constant published in Giorgi & 
Chameides (1985) as well as the standard Henry's Law constants from the NIST Chemistry 
Webbook (http://webbook.nist.gov/chemistry/). Since the parameterization of equation 
( 12T]) is primarily for modeling Earth's atmosphere, the specific reduction factor should 
be applied when it is reasonable to believe the hydrological cycle is reduced on an exoplanet 
(e.g., Tian et al. 2010). 

The Type 3 lower boundary condition has assigned flux that represents the surface 
emission. Note that the dry deposition (Type 2) and assigned flux (Type 3) lower boundary 
conditions can be used at the same time, but assigned number density (Type 1) boundary 
condition over-rules other lower boundary conditions. For example, S0 2 may be deposited 
to the surface at a rate proportional to the number density of the bottom layer, and also be 
emitted from the surface to the atmosphere at an assigned flux. 

We finish this section with a definition of the so-called "redox balance" , the requirement 
of which arises when the fixed mixing ratio lower boundary condition (Type-1) is used in 
the photochemistry model. Following Kasting & Brown (1998), Zahnle et al. (2006), and 
Segura et al. (2007), we define H 2 0, N 2 , C0 2 , S0 2 as redox neutral, and assign the redox 
number (7V) of any H-O-C-N-S molecule as the number of hydrogen in excess, i.e., 

ft(H ai O a2 N a3 C a4 S a5 ) = a x - 2a 2 + 4a 4 + 4a 5 . (22) 

Note that our definition of the redox number differs from Zahnle et al. (2006) and Segura 
et al. (2007) in the way that we count the number of H and they count the number of H 2 . 
With our definition, for example, the redox number of H 2 S is 6, and the redox number of 
H 2 S0 4 is -2. The redox balance says that the total redox influx to the atmosphere (i.e., 
surface emission) should be balanced by the total redox outflux from the atmosphere (i.e., 
atmospheric escape, dry deposition, and wet deposition), otherwise the atmosphere is being 
oxidized or reduced. The redox balance is equivalent to the conservation of total budget of 
hydrogen in the atmosphere, and is equivalent to the conservation of the total number of 
electrons in the atmosphere. One might consider the redox balance to be redundant. The 
redox balance is redundant if the lower and upper boundary conditions for all species are 
specified in fluxes (and not in mixing ratios). In many applications, however, it is useful 
to use the Type-1 lower boundary condition, which specifies mixing ratio. In that case, the 
mass conservation would not necessarily guarantee the redox balance. After all, an imbalance 
in the redox budget for a steady-state solution indicates either bugs in chemistry kinetics 
or unphysical boundary conditions. We explicitly check the redox balance for all our model 
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outputs. 

3. Model Validation 

We validate the photochemistry model by computing atmospheric compositions of cur- 
rent Earth and Mars. Due to the one- dimensional nature of the photochemistry model, we 
intend to compare with global averages of observations, if available. In most cases, however, 
the global average is only an order-of-magnitude representation of chemical composition of 
Earth's and Mars' atmospheres, because of strong diurnal and spatial variations. We should 
therefore expect our one-dimensional chemical-transport model to match quantitatively with 
the global average values to within a factor of two or three. For Mars, vertical profiles of 
most trace gas content have not yet been measured. We therefore compare our results with 
other photochemistry model outputs from the literature. 

3.1. Transport-Only Model for an Earth-like Atmosphere 

As a test of transport-related schemes of the photochemical model, we first consider a 
transport-only case. In theory, for species whose removal timescales are significantly larger 
than their transport timescales, such as C0 2 in Earth's atmosphere, their mixing ratios do 
not change with altitude, i.e., they are well-mixed. 

We turn off the chemical network deliberately and compute a model using Earth's tem- 
perature profile to test the eddy diffusion transport, molecular diffusion transport, condensa- 
tion, and rainout schemes of our chemical-transport model. A valid transport model should 
preserve hydrostatic equilibrium, predict well mixed mixing ratios for long-lived species, and 
predict a mixing ratio gradient for species that is rapidly removed in the atmosphere. We 
have verified that the transport scheme, although written in terms of number density rather 
than mixing ratio, preserves hydrostatic equilibrium. The fact that the code maintains hy- 
drostatic equilibrium indicates that the transport scheme is numerically correct. We verify 
that long-lived species are well mixed, such as C0 2 throughout the atmosphere, as shown in 
Figure [TJ We assign a mixing ratio for C0 2 at the surface and the code is able to predict a 
well-mixed vertical profile. We also verify that for H 2 molecular diffusion tends to increase 
the mixing ratio when the molecular diffusion coefficient is comparable with the eddy diffu- 
sion coefficient (i.e., the homopause; see Figured]). With a diffusion-limited escape, the effect 
of molecular diffusion and the effect of escape on the H 2 vertical profile near the homopause 
largely cancel out (Figured]), consistent with the definition of diffusion-limited escape. We 
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also show the behaviour of H 2 S, a poorly mixed species, to compare with well-mixed C02. 
In this simplified model the only way to remove H 2 S is photolysis, which requires photons 
with wavelengths shorter than 260 nm. Most of the photons in this wavelength range are 
effectively shielded by the major species 2 , so that removal rate of H 2 S by photolysis near 
the surface is smaller than that in the stratosphere (see Figured]). 

Water vapor content in a rocky exoplanet's atmosphere is controlled by the water reser- 
voir at the planetary surface, vertical transport, and condensation. In our model, the mix- 
ing ratio of water vapor at the surface is assigned according to appropriate temperature- 
dependent relative humidity. Water vapor can be transported up into the atmosphere by 
eddy diffusion, and as the temperature decreases the atmosphere may become supersaturated 
in water as the altitude increases. For an Earth-like planet atmosphere, it is appropriate to 
assume a water vapor mixing ratio of 0.01 at the surface (which corresponds to a relative 
humidity of about 60%). The condensation scheme becomes effective when water vapor is 
supersaturated, which keeps the water vapor profile along the saturation profile. Such an 
approach to computing the atmosphere water vapor content is commonly adopted by other 
previously described photochemistry models of terrestrial planets (e.g., Nair et al. 1994; 
Yung & Demore 1999; Zahnle et al. 2008). 

The one-dimensional transport model tends to saturate the tropopause and over-predict 
the amount of water vapor in the stratosphere. In fact, when using the US Standard Atmo- 
sphere 1976 as the temperature profile, the tropopause temperature is 217 K, much warmer 
than the required "cold trap" temperature (~ 200 K) in order to maintain a dry stratosphere 
(water vapor mixing ratio in a few ppm). Common 1-D transport-condensation schemes may 
over-predict the amount of water vapor above the cold trap by a few orders of magnitude, and 
therefore over-predict the strength of HO x cycle in the stratosphere. This is a well known 
problem in modeling Earth's atmosphere, mainly due to spatial and temporal variability 
of tropopause temperatures (K. Emmanuel, 2012, private communication). Tropospheric 
convection is most effective in transporting water vapor in tropics, where the tropopause 
temperature is the lowest and the number density of water vapor at the tropopause is the 
lowest compared with other latitudes. It is therefore likely that the tropopause is highly 
unsaturated as a global average. We use the temperature profile of the equatorial region in 
January (see Figured]), and limit the water vapor saturation ratio to within 20%. In this 
way we reproduce the water vapor profile of the US Standard Atmosphere 1976, which has 
a dry stratosphere (Figured]). The classic Manabe-Wetherald relative humidity profile for 
1-D photochemistry models of Earth's atmosphere also has 20% relative humidity at the 
tropopause (Manabe & Wetherald 1967). 

As a result of water vapor condensation and rain, soluble gases are removed from the 
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atmosphere by wet deposition (rainout). We expect to see a decreasing slope in mixing ratio 
of the species being rained out, which depends on the rainout rate, itself dependent on the 
solubility of the species. For example, S0 2 is much more soluble than H 2 S, so that S0 2 is 
more effectively rained out in the troposphere and cannot accumulate in the troposphere 
(see Figured]). 



3.2. Present-Day Earth 

We compute a photochemical model to simulate the present-day Earth, and compare 
our results with globally averaged measurements. The most important photochemical pro- 
cess in Earth's atmosphere is formation of the ozone layer, which leads to the temperature 
inversion in the stratosphere (e.g., Seinfeld & Pandis 2006). On the other hand, the tro- 
pospheric chemistry is mostly controlled by the hydrological cycle and surface emission of 
CH 4 , NO x (i.e., NO and N0 2 ) and S0 2 , and involves coupled processes of photochemistry 
and deposition. The reproduction of current Earth atmospheric composition is therefore a 
comprehensive validation of various aspects of the photochemistry model. 

We start with a nominal 80% N 2 and 20% 2 composition, having well-mixed C0 2 
of 350 ppm. We assume the water vapor mixing ratio at the surface is 0.01 (i.e., type- 
1 lower boundary condition), and let the water vapor be transported in the atmosphere 
and condense near the tropopause. The temperature profiles are assumed to be the 1976 
US Standard Atmosphere (subsequently referred to as "Model A") or a reference tropical 
atmosphere in January of 1986 from COSPAR0 International Reference Atmosphere (CIRA, 
subsequently referred to as "Model B"; ), and the eddy diffusion coefficient is adopted from 
Massie & Hunten (1981), who have derived it from vertical profiles of trace gases. We model 
the atmosphere from to 86 km altitude as 43 equally spaced layers. The incoming solar 
radiation is cut off for wavelengths shorter than 100 nm to account for the thermosphere 
absorption. The photolysis rates at the top layer of the atmosphere are tabulated in Table 
[7J Key trace species are assigned emission rates from the surface according to typical global 
values (Seinfeld & Pandis 2006). Emission rates and dry deposition velocities of these species 
are tabulated in Table [H 

Our models correctly simulate the photochemistry in Earth's stratosphere. First, the 
formation of the ozone layer is correctly predicted. The modeled vertical profiles and amounts 
of ozone are consistent with the globally averaged measurements, as shown in Figure [2j 
Second, the nitrogen oxide cycle and the HO x cycle are properly simulated. In Earth's 
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Table 1: Photochemistry model validation. For the surface emission and dry deposition 
velocities, measured surface mixing ratios are compared to our models using the 1976 US 
Standard Atmosphere (Model A) and using the reference tropical temperature profile (Model 
B). 



Species 


Surface Emission 3, 


Vdep 


Surface Mixing 


Ratio 






(molecule cm -2 s _1 ) 


(cm s _1 ) 


Measured 6 


Model A 


Model B 


CO 


3.7 x 10 11 


0.03 b 


40-200 ppb 


113 ppb 


101 ppb 


CH 4 


1.4 x 10 11 


0.03 b 


700-1745 ppb f 


1939 ppb 


1235 ppb 


NH 3 


7.7 x 10 9 


l d 


0.1-10 ppb 


0.24 ppb 


0.24 ppb 


N 2 


1.0 x 10 9 


b 


276-315 ppb f 


302 ppb 


290 ppb 


NO 


7.0 x 10 9 


0.016 b 


0.02-10 ppb g 


0.024 ppb 


0.025 ppb 


so 2 


9.0 x 10 9 


l c 


30-260 ppt 


237 ppt 


239 ppt 


OCS 


5.0 x 10 7 


0.01 d 


510 ppt 


188 ppt 


185 ppt 


H 2 S 


2.0 x 10 s 


0.015 c 


1-13 ppt 


3.92 ppt 


3.62 ppt 


H 2 S0 4 


7.0 x 10 s 


l d 


5-70 ppt h 


127 ppt h 


126.7 ppt h 



"Typical globally averaged emission rates are taken from Seinfeld & Pandis (2006). 

fc Typical dry deposition velocities are taken from the compilation of Hauglustaine et al. (1994). 

c Typical dry deposition velocities are taken from the compilation of Sehmel (1980). 

d Dry deposition velocities are assumed by considering the solubility and reactivity of gases (see Appendix 
e Mixing ratios at the surface are taken from Seinfeld & Pandis (2006). 
/Ran ges are from the preindustrial mixing ratio to the present-day value. 
9 Ranges are mixing ratio of NO x at the planetary boundary layer. 
''Mixing ratio of S0 4 ~ in both gaseous and aqueous phases. 

stratosphere, chemical species are rapidly converted to each other within the NO x group and 
the HO x group, and both cycles lead to catalytic destruction of ozone (e.g., Seinfeld & Pandis 
2006). We reproduce the vertical profiles of NO x and HO x in the stratosphere (see Figure |2J, 
qualitatively consistent with mid-latitude measurements. It appears that the models tend 
to over-predict the amounts of NO in the lower stratosphere. Third, for relatively long lived 
species in the stratospheres like N 2 and CH 4 , our models produce correct altitude gradient 
of mixing ratios. 

The key element of the tropospheric chemistry is the production of the hydroxyl radical 
OH, because OH is the major removal pathway for most species emitted from the surface. 
It has been established that OH in Earth's troposphere is produced by reactions between 
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0( 1 D) and H 2 0, and 0( 1 D) is produced by photolysis of 3 (Seinfeld & Pandis, 2006). In 
turn, the main source of tropospheric ozone is the NO x cycle, with non-neglible effects of 
HO x cycle and hydrocarbon chemistry as well. We find the surface ozone mixing ratio is 19 
ppb and 22 ppb for Model 1976 and Model 1986 respectively, and the surface OH number 
density is 2.0 x 10 6 cm -3 for Model 1976 and 2.2 x 10 6 cm -3 for Model 1986, very close to 
the commonly adopted OH number density of 1.0 x 10 6 cm -3 . We also confirm that the 
steady-state mixing ratios of major trace gases near the surface are consistent with ground 
measurements (see Tabled]). The main removal mechanism of CO, CH 4 , and NH 3 is through 
reactions of OH, and N 2 is long-lived in the troposphere and is transported up to the 
stratosphere. Considering the significant spatial and temporal variability of the amount of 
OH, we conclude that our photochemistry model correctly computes the chemistry in Earth's 
troposphere. 

Our photochemistry model correctly treats the sulfur chemistry in Earth's atmosphere 
(see Table[T]). All sulfur-bearing emission, if not deposited, is oxidized in multiple steps in the 
troposphere and eventually converted into sulfate. Sulfite (S 4+ ) and sulfate (S 6+ ) are soluble 
in water and effectively removed from the atmosphere by rainout. Our photochemistry 
models simulate these processes, and find steady-state mixing ratios of sulfur-bearing species 
(e.g., H 2 S and S0 2 ) consistently with the ground measurements. The models also predict 
the saturation of sulfate and the formation of sulfate aerosols as expected. The modeled 
mixing ratio of sulfate is slightly larger than observations, which might be related to the fact 
that we do not consider the hydration of H 2 S0 4 , nor the sulfate-facilitated cloud formation 
in our models. 

In summary, we validate our photochemistry model by reproducing the chemical com- 
position of current Earth's atmosphere (see Table [1] and Figure [2]). The photochemistry 
model successful in predicting the formation of the ozone layer, treating key chemical cycles 
in both the stratosphere and the troposphere, computing oxidation of hydrocarbon, ammo- 
nia and sulfur-bearing species in the troposphere, and transporting long-lived species from 
the troposphere to the stratosphere. These aspects involve all chemical kinetics, photolysis, 
and transport processes, which not only verify that our photochemistry model is suitable 
for applications in oxidizing atmospheres, but also allows our model to be applied to other 
atmospheric scenarios. 

3.3. Present-Day Mars 

We validate our photochemistry model by simulating the atmosphere of present-day 
Mars. The current atmosphere of Mars is a thin C0 2 -based atmosphere and its bulk chemical 
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composition is known from extensive ground-based observations (see Krasnopolsky 2006 and 
references therein). We validate our photochemistry code by reproducing the observed mixing 
ratios of CO and 3 , and comparing with previous results regarding Martian atmosphere 
photochemistry. Key parameters and results of the present- Mars atmosphere model are 
shown in Table [2] and Figure O We emphasize that it is important to use temperature- 
dependent UV cross sections for modeling the photolysis and UV penetration in the Martian 
environment (see Anbar et al. 1993 for an error analysis of temperature-dependence of C0 2 
in the Martian atmosphere). 

The most important aspect of Martian atmospheric photochemistry is the stabilization 
of the C0 2 atmosphere by H, OH, and H0 2 (commonly referred to as odd hydrogens; Nair et 
al. 1994; Krasnopolsky 2006; Zahnle et al. 2008). Odd hydrogens, as catalysts, facilitate the 
recombination of photochemically produced CO and 2 and maintain the C0 2 dominance 
of the Mars' atmosphere (e.g. Nair et al. 1994; Krasnopolsky 2006; Zahnle et al. 2008). A 
pure C0 2 atmosphere is unstable against photolysis, because C0 2 can only be restored with 
a slow three-body reaction CO + O + M — > C0 2 + M (e.g., Yung & DeMore 1999). Odd 
hydrogen species, including H, OH and H0 2 , are produced by photolysis of water vapor; and 
trace amounts of odd hydrogens can effectively stabilize the C0 2 dominant atmosphere (e.g. 
Nair et al. 1994; Krasnopolsky 2006; Zahnle et al. 2008). Nonetheless, one-dimensional 
photochemistry models tend to over- address the problem, predicting amounts of CO, 2 
and 3 several times smaller than the observed values (e.g. Nair et al. 1994; Krasnopolsky 
2006). By considering a slow dry deposition of H 2 2 and 3 to the surface, Zahnle et al. 
(2008) are able to predict the amount of CO, 2 and H 2 that match observations. We follow 
the assumptions of Zahnle et al. (2008), and confirm the finding of appropriate amounts 
of 2 , CO and H 2 (see Table |2] and Figure |3]). We therefore reproduce the photochemical 
stability of the current Martian atmosphere. 

The timescale of CH 4 removal can also be used for model validation. Rapid variation 
(over several years) of the amount of CH 4 in Mars' atmosphere has been reported (e.g., 
Mumma et al. 2009), but the modeling of coupled general circulation and gas-phase chem- 
istry find no known gas-phase chemistry path that allow such a rapid removal (e.g., Lefevre 
& Forget 2009). Based on our fiducial model of Mars' atmosphere, the loss timescale for 
CH 4 is computed to be about 240 years, within the same order of the magnitude of Lefevre 
& Forget (2009). 

In summary, we validate our photochemistry model by reproducing the atmospheric 
compositions of current Earth and Mars. We find that the model gives consistent results 
compared to observations and previous photochemistry models. All physical and chemical 
processes including photolysis, chemical reactions, transports, condensation, and deposition 
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Table 2: Photochemistry model validation. For the upper boundary flux and dry deposition 
velocities, measured mixing ratios of major trace gases on Mars are compared to our model. 



Species 


Upper Boundary Flux a 
(molecule cm -2 s _1 ) 


v b 

'DEP 

(cm s _1 ) 


Column Averaged Mixing Ratio 
Measured 6 Modeled 


o 2 








1200-2000 ppm 


1545 ppb 


CO 


-2.0 x 10 7 





800 ppm 


572 ppb 


H 2 


3.6 x 10 8 





17 ppm 


23 ppb 


H 2 2 





0.02 


~ 40 ppb f 


18 ppb 


3 





0.02 


~ 120 ppb f 


18 ppb 



a H 2 upper boundary flux is computed from diffusion-limited escape velocity after finding the steady state. 
Besides tabulated values, influx of N of 2.0 x 10 6 molecule cm" - 2 s"\ NO of 2.0 x 10 7 molecule cm s , 
and outflux of O of 2.0 x 10 7 molecule cm~ 2 s _1 , are considered as input to the model. 

h Dry deposition velocity is assumed to be 0.02 for reactive species, including H, O, 0( 1 D), 3 , OH, H0 2 , 
H 2 2 , CH0 2 , CH 2 2 , CH 3 2 , CH 4 2 , NH 3 , N0 3 , N 2 5 , HN0 3 , and HN0 4 (e.g., Zahele et al. 2008). Dry 
deposition velocity is assumed to be zero for all other species. 

e Mars data are from the compilation of Krasnopolsky (2006). 

Mixing ratios of 3 and H 2 2 have significant diurnal, seasonal and latitudinal variations, 
are rigorously tested in these examples. 

4. Exoplanet Benchmark Cases 

We now present three benchmark atmospheric scenarios of rocky exoplanets and sum- 
marize the effects of key photochemistry processes. The goal is to provide baseline models 
to assess the stability of molecules in different kinds of atmospheres in order to: identify the 
dominant stable molecules; calculate the lifetime of spectrally significant gases; and identify 
the amounts of the main reactive species that control molecule lifetimes. The benchmark 
cases are also intended to serve as the test cases for independent photochemistry models for 
rocky exoplanet atmospheres. 

The key to assess molecular stability is the oxidation state of the atmosphere, because 
the main reactive species are linked to the oxidation state of the atmosphere. In an oxidizing 
atmosphere, OH and O are created by photochemistry and are the main reactive radicals. 
In a reducing atmosphere, H, also created by photochemistry, is the main reactive species. 
Although we expect the atmospheric composition will be highly varied, based on the nearly 
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continuous range of masses and orbits of exoplanets, we believe that the primary dimension 
of chemical characterization for terrestrial exoplanet atmospheres is their oxidation states. 
For the benchmark cases, therefore, we have chosen three endmembers in terms of atmo- 
spheric oxidizing power. The scenarios are a reducing (90%H 2 -10%N 2 ) atmosphere, a weakly 
oxidizing N 2 atmosphere (> 99%N 2 ), and a highly oxidizing (90%CO 2 -10%N 2 ) atmosphere. 
We consider Earth-like volcanic gas composition that consists of C0 2 , H 2 , S0 2 , CH 4 , and 
H 2 S, with emission rates comparable to current Earth. We assume that the planet surface 
has a substantial fraction of its surface covered by a liquid water ocean so that water is 
transported from the surface and buffered by the balance of evaporation/condensation. Key 
assumptions of the parameters of the three atmospheric scenarios are summarized in Table 
[3] and rationals of important model parameters are given in § 14.31 

We now describe the key results of the three atmosphere scenarios, for the 1-bar atmo- 
sphere on an Earth-size and Earth-mass habitable rocky exoplanet around a Sun-like star. 
We also describe qualitatively the behaviors of the atmospheres in the habitable zone of a 
quiet M dwarf, based on our numerical explorations. Mixing ratios of emitted gases, pho- 
tochemical products, and reactive species in the three scenarios are tabulated in Table @] 
and shown in Figure HI Schematic illustrations of key non-equilibrium processes in the three 
scenarios are shown in Figure [51 We start with a qualitative overview of the key results, and 
then present the chemical properties of the three benchmark scenarios. 

4.1. General Results 

We here list our main findings on general chemistry properties of atmospheres on hab- 
itable terrestrial exoplanets. 

Atomic hydrogen (H) is a more abundant reactive radical than hydroxyl radical (OH) in 
anoxic atmospheres. Atomic hydrogen is mainly produced by water vapor photodissociation; 
and in anoxic atmospheres the main ways to remove atomic hydrogen is its recombination 
and reaction with CO. This is in contrast to oxygen-rich atmospheres (e.g., current Earth's 
atmosphere) in which H is quickly consumed by 2 . As a result, removal of a gas by H is 
likely to be an important removal path for trace gases in an anoxic atmosphere. Atomic 
oxygen is the most abundant reactive radical in C0 2 dominated atmospheres. Due to the 
photochemical origin of the reactive species including H, OH, and O, their abundances in the 
atmosphere around a quiet M dwarf are 2 orders of magnitude lower than their abundances 
around a Sun-like star. 

Dry deposition velocities of long-lived compounds, notably major volcanic carbon com- 
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Table 3: Basic parameters for the terrestrial exoplanet benchmark scenarios. The benchmark 
scenarios are H 2 , N 2 , C0 2 atmospheres on habitable terrestrial exoplanets with Earth-like 
volcanic emissions. Note that we do not consider any known biosignature gas emission or 
biotic contribution to the dry deposition velocities. 



Parameters 


Reducing 




Weakly Oxidizing 


Highly Oxidizing 


Main component 


90%H„. 10%N„ 
2' 2 




> 99%N 2 


90%CO o , 10%N„ 
2' 2 


Mean molecular mass 


4.6 




28 


42.4 


Planetary parameters 










Stellar type 


G2V 




G2V 


G2V 


Semi-major axis 


1.6 AU 




1.0 AU 


1.3 AU 


Mass 


M e 




M e 


M e 


Radis 






A® 


R(B 


Temperature profile 










Surface temperature 


288 K 




288 K 


288 K 


Surface pressure 


10 5 Pa 




10 5 Pa 


10 5 Pa 


Tropopause altitude 


120 km 




13.4 km 


8.7 km 


Temperature above tropopause 


160 K 




200 K 


175 K 


Maximum altitude 


440 km 




86 km 


51 km 


Eddy diffusion coefficient 










In the convective layer 


6.3 X 10 5 cm 2 s~ 


1 


1.0 X 10 5 cm 2 s _1 


6.8 X 10 4 cm 2 s" 1 


Minimum 


2.5 X 10 4 cm 2 s" 


1 


3.9 X 10 3 cm 2 s" 1 


2.7 X 10 3 cm 2 s" 1 


Altitude for the minimum 


107 km 




17.0 km 


11.6 km 


Near the top of atmosphere 


7.1 X 10 5 cm 2 s" 


1 


1.1 X 10 5 cm 2 s _1 


7.6 X 10 4 cm 2 s _1 


Water and rainout 










Liquid water ocean 


Yes 




Yes 


Yes 


Water vapor boundary condition 


/ (H 2 0) = 0.01 




/(H 2 OJ = 0.01 


/(H 2 0) = 0.01 


Rainout ratc a 


Earth-like 




Earth-like 


Earth-like 


Gas emission* 3 










co 2 


3 X 10 11 cm~ 2 




3 X 10 11 cm- 2 s" 1 


N/A 


H 2 


N/A 




3 X 10 10 cm -2 s _1 


3 X 10 10 cm -2 s _1 


so 2 


3 X 10° cm -2 




3 X 10° cm -2 


3 X 10 9 cm -2 


CH 4 


3 X 10 s cm -2 s~ 


1 


3 X 10 s cm- 2 


3 X 10 8 cm~ 2 


H 2 S 


3 X 10 S cm -2 




3 X 10 S cm -2 


3 X 10 8 cm~ 2 



Dry deposition velocity 11 








H 2 









CH 4 









C 2 H 6 


1.0 X 10" 


■5 


(Assumed) 


CO 


1.0 X 10" 


■8 


cm s _1 (Kharecha et al. 2 


CH 2 


0.1 cm s~ 


-1 


(Wagner et al. 2002) 


co 2 


1.0 X 10" 


■4 


cm s _1 (Archer 2010) 


o 2 









°3 


0.4 cm s" 


-1 


(Hauglustainc et al. 1994) 


H 2 2 


0.5 cm s~ 


-1 


(Hauglustainc et al. 1994) 


H 2 S 


0.015 cm 


s~ 


_1 (Sehmel 1980) 


so 2 


1.0 cm s~ 


-1 


(Sehmel 1980) 


S S (A) 


0.2 cm s~ 


-1 


(Sehmel 1980) 


H,S0 4 (A) 


0.2 cm s" 


-1 


(Sehmel 1980) 



"Rainout rates for H 2 , CO, CH 4 , C 2 H 6 , and 2 are generally assumed to be zero to simulate an ocean surface 
saturated with these gases on an abiotic exoplanet. 

& The volcanic gas emission rates from the planetary surface are assigned for each model scenario. H 2 
emission is not explicitly considered because the surface has a large water reservoir, i.e., an ocean. 

c We here list the dry deposition velocities (with references) for emitted gases and their major photochemical 
byproducts, and dry deposition velocities that are important for the mass and redox balance of the atmo- 
sphere. Dry deposition velocities are assumed to be identical for the three scenarios. C 2 H 6 dry deposition 
velocity is assumed to take into account the loss of carbon due to organic haze formation and deposition. 
C0 2 dry deposition velocity is estimated from a 10, 000- year lifetime that matches with the lifetime of silicate 
weathering on Earth (Archer 2010; S. Solomon, private commonication). 
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pounds including methane, carbon monoxide, and carbon dioxide, have significant effects on 
the atmospheric oxidation states. The specific choice of dry deposition velocities for emitted 
gases and their major photochemical byproducts in the atmosphere is critical to determine 
the atmospheric composition on terrestrial exoplanets. 

Volcanic carbon compounds (i.e., CH 4 and C0 2 ) are chemically long-lived and tend to be 
well mixed in terrestrial exoplanet atmospheres, whereas volcanic sulfur compounds (i.e., H 2 S 
and S0 2 ) are short-lived. CH 4 and C0 2 have chemical lifetime longer than 10,000 years in all 
three benchmark atmospheres ranging from reducing to oxidizing, implying that a relatively 
small volcanic input can result in a high steady-state mixing ratio. The chemical lifetime 
CO, another possible volcanic carbon compound, ranges from 0.1 to 700 years depending on 
the OH abundance in the atmosphere. 

We find abiotic 2 and 3 , photochemically produced from C0 2 photolysis, build up 
in the 1-bar C0 2 dominated atmosphere if volcanic emission rates of reducing gases (i.e., H 2 
and CH 4 ) is more than one order of magnitude lower than current Earth's volcanic rates. 
Abiotic 2 can be a false positive for detecting oxygenic photosynthesis, but the combination 
of 2 /0 3 and reducing gases remains rigorous biosignature. 



4.2. Chemistry of H 2 , N 2 , and C0 2 Dominated Atmospheres 

4-2.1. H 2 Dominated Atmospheres 

The main reactive agent in the H 2 dominated atmosphere is atomic hydrogen (H). The 
abundance of atomic hydrogen is five orders of magnitudes higher than that of hydroxyl 
radical (Figure H]). The source of both H and OH is water vapor photodissociation. In the 
H 2 dominated atmosphere, most of the OH molecules produced from water vapor photodis- 
sociation react with H 2 to reform H 2 and produce H through the reactions: 

H 2 + hv — H + OH, (CI) 

OH + H 2 — ► H + H 2 0. (C2) 

Therefore the abundance of atomic hydrogen is much higher than that of OH. We note that 
H production via water photodissociation is much more efficient than the direct photodisso- 
ciation of H 2 , which requires radiation in wavelengths less than 85 nm. As water vapor is the 
primary source of H and OH in anoxic atmospheres, the amounts of H and OH depends on 
the mixing ratio of water vapor above the cold trap, which is in turn sensitively controlled 
by the cold trap temperature. Water vapor mixing ratio spans 3 orders of magnitudes for 
the cold trap temperature ranging in 160 - 200 K; consequently, the number densities of H 
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Table 4: Atmospheric compositions of terrestrial exoplanet benchmark scenarios. For the 
surface emission in 1-bar H 2 , N 2 , and C0 2 dominated atmospheres, mixing ratios of emitted 
gases, photochemical products and reactive agents are computed by the photochemistry 
model. 



Scenario 


Column Averaged Mixing Ratio 








Emitted Gases 


Photochemical Products 3, 


Reactive Agents b 




Reducing 


C0 2 : 8.9 X 10" 5 


CO: 8.0 X 10~ 6 


H: 2.1 X 10 -9 




90%H 2 , 10%N 2 


S0 2 : 9.9 X lCT 12 


C 2 H 6 : 4.7 X 10" 10 


OH: 2.0 X 10~ 14 




CH 4 : 5.9 X 1CT 6 


S g : 3.5 X 10- 10 


O: 1.2 X 10" 11 






H 2 S- 9 1 X KT 10 


CH 2 0: 2.9 X 10- 10 
CH 4 0: 5.6 X 10 -11 


0( 1 D): 2.2 X 10" 


21 


Weakly oxidizing 


C0 2 : 1.3 X 10" 4 


CO: 1.7 X 10 -7 


H: 1.2 X 10 -9 




N 2 


H 2 : 4.5 X 10~ 4 


C 2 H 6 : 9.0 X 10~ 9 


OH: 9.3 X 10~ 14 






S0 2 : 8.9 X 10- 12 


CH 4 0: 1.4 X 10~ 9 


O: 6.5 X 10~ 10 






CH 4 : 3.1 X 1CT 5 


2 : 3.4 X lO" 10 


0( 1 D): 1.8 x 10" 


20 




H 2 S: 1.1 X 10- 14 


S g : 3.0 X 10- 10 
CH 2 0: 4.0 X 10- 11 
C 2 H 2 : 1.5 X lO" 11 






Highly oxidizing 


H 2 : 1.0 X 10~ 3 


CO: 7.7 X 10~ 3 


H: 6.0 X lO" 11 




90%CO 2 , 10%N 2 


S0 2 : 1.6 x 10- 10 


2 : 6.4 X 10~ 7 


OH: 7.8 X 10~ 15 






CH 4 : 3.7 X 10- 5 


C 2 H 6 : 6.1 X 10- 10 


O: 2.0 X 10~ 8 






H 2 S: 1.4 X 10~ 10 


H 2 2 : 3.7 X lO" 10 
H 2 S0 4 : 5.0 X 10- 11 
CH 2 0: 2.5 X 10- 12 


0( 1 D): 3.0 x 10" 


18 



a Species produced in the atmosphere from the emitted volcanic gases via photochemistry and subsequent 
series of chemical reactions, listed in the order of decreasing abundance. S 8 and H 2 S0 4 mixing ratios include 
both the gas phase and the condensed phase; whereas the condensed phase is found to contain more than 
85% of mass. 

'Common reactive agents in the atmosphere are H, OH, O, 0( 1 D). We list the abundance of these gases in 
all three scenarios for as useful reference for the future assessment of chemical lifetime of trace gases. 
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and OH in the atmosphere can easily vary by one order of magnitude, depending on the cold 
trap temperature. The removal of atomic hydrogen is mainly by recombination to H 2 , which 
can be more efficient with the presence of CO, via 



As a result of low OH abundance in the H 2 dominated atmosphere, both CH 4 and CO 
are long-lived and therefore well mixed. CH 4 emitted from the surface is slowly oxidized in 
the atmosphere into methanol, which gives methane a chemical lifetime of 8 x 10 4 years. CO 
is produced by C0 2 photodissociation in our benchmark model; it is emitted by volcanoes 
on Earth at a much lower rate than C0 2 , but it can presumably be the main carbon-bearing 
gas produced by volcanoes if the upper mantle is more reduced than the current Earth's 
(e.g., Holland 1984). We find that CO is long-lived in the H 2 atmosphere as well. 

The lack of efficient atmospheric sink of CH 4 and CO in the H 2 atmosphere implies that 
surface sink, if any, is the major sink for these two carbon compounds. CH 4 and CO have zero 
or very small dry deposition velocities on an abiotic planet, so they can build up to significant 
amounts in the H 2 atmosphere (Figure H]). If a nonzero deposition velocity is adopted for 
CH 4 and CO, their steady-state mixing ratios will be much lower than the benchmark model. 
For example, using a deposition velocity for CH 4 standard on Earth (~ 0.01 cm s _1 ) results 
in a mixing ratio of less than 1 ppb for CH 4 , compared to a mixing ratio of 6 ppm in 
the benchmark model. The dry deposition velocity is indeed the controlling factor for the 
steady-state abundance of the long-lived carbon compounds. In comparison, the emitted 
sulfur compounds (H 2 S and S0 2 ) are short-lived. In the H 2 atmosphere, sulfur emission 
from the surface is readily converted to elemental sulfur aerosols (S 8 ) in the atmosphere. 

Interestingly, C0 2 is fairly long-lived and well mixed in the H 2 atmosphere, meaning 
that the reduction of C0 2 in the H 2 dominated atmosphere is not efficient (Figure H]). Only 
1/3 of emitted C0 2 is reduced in the atmosphere, and the rest is deposited to the surface by 
dry and wet deposition. It is the balance between surface emission and surface deposition 
that sets the steady-state mixing ratio of C0 2 . 




(C3) 




(C4) 



4-2.2. N 2 Dominated Atmospheres 



Both reducing radicals (i.e., H) and oxidizing radicals (i.e., O and OH) are relatively 
abundant in the N 2 atmosphere compared with the H 2 atmosphere (Figure H]). Like in H 2 
dominated atmospheres, H abundance is orders of magnitude larger than OH abundance, 
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because most of OH molecules from water photolysis react with H 2 to reform H 2 and 
produce H. The molecular hydrogen that consumes OH and boosts H is emitted volcanically, 
so a general N 2 atmosphere can be more oxidizing (i.e., having lower H and higher OH) if 
there is a lower volcanic H 2 emission and mixing ratio than the benchmark model. 

An important feature of the N 2 dominated atmosphere is that both H and OH are 
relatively abundant near the surface. Comparing the N 2 dominated atmosphere with the H 2 
dominated atmosphere (Figure SJ), we find that both the OH and the H number densities are 
higher near the surface due to the lower H 2 number density. The relatively high OH number 
density leads to relatively fast removal of CO by 

CO + OH — > C0 2 + H, (C5) 

as shown in Figure HI With a low CO abundance the recombination of H via reactions AO 34 - 
IC4j) is inefficient. Counter- intuitively, a high OH number density helps preserve H in this 
specific case. This example shows the complexity and the nonlinearity of an atmospheric 
chemical network. The feature of simultaneous high OH and H abundances near the surface 
is sensitive to the specification of surface hydrogen emission and eddy diffusion coefficients 
(see section B~B1 for relevant rationale). 

The chemical lifetimes of CH 4 and CO mainly depend on the amount of OH. In the N 2 
atmosphere, CH 4 is well mixed because its chemical lifetime is long. CH 4 is photodissociated 
and oxidized slowly in the atmosphere into methanol (CH 4 0) with a chemical lifetime of 
~ 6 x 10 4 years. The photolysis of methane is a secondary source of atomic hydrogen, 
which concentrates at the pressure level of 10 Pa. Interestingly, methane photolysis causes 
the apparent trough of 2 mixing ratio profile at ~ 10 Pa. It is the shielding of the UV 
radiation that dissociates methane by methane itself that determines this pressure level (see 
Appendix D for an analytical formula for assessing the pressure level at which the photolysis 
of a certain gas is important). 

C0 2 is actively photodissociated into CO and O in the upper atmosphere, but most of 
the CO produced is efficiently converted back to C0 2 by reacting with OH (reaction IC 5p . 
In equilibrium the net chemical removal of C0 2 is minimal, so the C0 2 mixing ratio is set 
by the balance between emission and deposition. The steady-state amount of C0 2 can be 
directly estimated by the mass balance between emission and deposition, viz., 

f=/i = rw' (23) 

where the overall mixing ratio (F) is the same as the near-surface mixing ratio (fi), and 
$1/2 is the surface emission rate. For C0 2 in the N 2 dominated atmosphere, using equation 
f )23|) . with Vdep(C0 2 ) = 1 x 10~ 4 cm s -1 , we find the steady-state mixing ratio is 1.3 x 10~ 4 , 
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consistent with the full photochemical model (see Table H] and Figure H]). The steady-state 
mixing ratio of a long-lived volcanic gas that is primarily removed by surface deposition is 
inversely proportional to its dry deposition velocity. 

Sulfur-bearing gases emitted from the surface are effectively converted into elemental 
sulfur and sulfuric acids. Elemental sulfur, mainly in the condensed phase (i.e., aerosols), 
is the major sulfur-bearing species in the steady state because of relatively high H 2 mixing 
ratio in the benchmark model. In separate numerical simulations we find that sulfuric acid 
aerosols may outnumber elemental sulfur aerosols when the H 2 emission is reduced by more 
than 1 order of magnitude and the atmosphere is more oxidizing than the benchmark model. 

4-2.3. CO 2 Dominated Atmospheres 

CO 2 photodissociation produces CO, O, and 2 . Atomic oxygen is the most abundant 
reactive radicals in the C0 2 dominated atmosphere (Figure HJ), because H is readily removed 
by 2 and OH is readily removed by CO. As a result of low H and OH, CO is long-lived in 
the atmosphere and can build up to very high mixing ratios in the atmosphere depending 
on its dry deposition velocity; and CH 4 emitted from the surface is also long-lived in the 
atmosphere with a chemical lifetime of ~ 6 x 10 4 years. The steady-state abundance of CO 
and CH 4 is therefore controlled by their dry deposition velocities. Most of the emitted S0 2 
is deposited to the surface, because there are very few OH or O radicals near the surface. A 
small fraction of the S0 2 is transported upwards and converted into sulfuric acid aerosols in 
the radiative layer. 

Based on our benchmark model of C0 2 dominated atmospheres, we now revisit whether 
photochemically produced 2 can cause false positive for detecting oxygenic photosynthesis. 
Oxygen and ozone are the most studied biosignature gases for terrestrial exoplanet charac- 
terization (e.g., Owen 1980; Angel 1986; Leger et al. 1993, 1996; Beichman et al. 1999). 
One of the main concerns of using 2 /0 3 as biosignature gases is that 2 may be produced 
abiotically from photodissociation of C0 2 . A number of authors have studied the abiotic 
production of oxygen in terrestrial atmospheres, either for understanding prebiotic Earth's 
atmosphere (e.g., Walker 1977; Kasting et al. 1979; Kasting & Catling 2003), or for assess- 
ing whether abiotic oxygen can be a false positive for detecting photosynthesis on habitable 
exoplanets (Selsis et al. 2002; Segura et al. 2007). It has been proposed that abiotic oxygen 
from C0 2 photodissociation is not likely to build up in the atmosphere on an planet having 
active hydrological cycle (Segura et al. 2007). 

We find that the steady-state number density of 2 and 3 in the C0 2 dominated 
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atmosphere is mainly controlled by the surface emission of reducing gases such as H 2 and 
CH 4 , and without surface emission of reducing gas photochemically produced 2 can build 
up in a 1-bar C0 2 dominated atmosphere. In addition to the benchmark model, we have 
simulated C0 2 dominated atmospheres with relatively low and zero emission rates of H 2 and 
CH 4 (Tableland Figure [6]). We see that the 2 mixing ratio near the surface increases 
dramatically in 1-bar C0 2 atmospheres when the emission of reducing gases decreases. 2 
is virtually nonexistent at the surface for the Earth-like emission rates of H 2 and CH 4 ; but 
2 mixing ratio can be as high as 10 -3 if no H 2 or CH 4 is emitted (Figure [6]). In particular, 
if no H 2 or CH 4 is emitted, the 3 column integrated number density can reach one third 
of the present-day Earth's atmospheric levels (Table EJ), which constitutes a potential false 
positive. 

In the 1-bar C0 2 dominated atmosphere 3 can potentially build up to a false-positive 
level even on a planet with active hydrological cycle. Segura et al. (2007) have based their 
conclusion on simulations of 20% C0 2 1-bar atmospheres with and without emission of H 2 
and CH 4 and simulations of 2-bar C0 2 atmospheres with emission of H 2 and CH 4 . We 
have been able to reproduce all results of Segura et al. (2007) quantitatively to within a 
factor of two. Where we differ from Segura et al. (2007) is that we successfully simulated 
high C0 2 1-bar atmospheres with minimal volcanic reducing gas emission (Figure [6]). This 
is a parameter space that Segura et al. (2007) did not cover, but we find that this is the 
parameter space for high abiotic 2 . 



4.3. Rationale of Model Parameters 

We here provide rationale for our specification of the atmospheric temperature profile, 
the eddy diffusion coefficients, and the dry deposition velocities. 

First, the surface temperatures for the three scenarios are assumed to be 288 K, and 
the semi-major axis of the planet is adjusted according to appropriate amounts of green- 
house effect. The semi-major axis around a Sun-like star for H 2 , N 2 , and C0 2 dominated 
atmospheres is found to be 1.6 AU, 1.0 AU, and 1.3 AU. We have compared the planetary 
thermal emission flux and the incidence stellar flux to determine the semi-major axis, a sim- 
ilar procedure as Kasting et al. (1993) taking into account C0 2 , H 2 0, and CH 4 absorption, 
H 2 collision-induced absorption, and 50% cloud coverage and a Bond albedo of 30%. The 
temperature profiles are assumed to follow appropriate dry adiabatic lapse rate (i.e., the 
convective layer) until 160 K (H 2 atmosphere), 200 K (N 2 atmosphere), and 175 K (C0 2 
atmosphere) and to be constant above (i.e., the radiative layer). We simulate the atmosphere 
up to the altitude of about 10 scale heights, or the pressure level of 0.1 Pa. The adopted 
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Table 5: Mixing ratios of 2 and 3 and redox budget for C0 2 -dominated atmospheres on 
rocky exoplanets having different surface emission of reducing gases. The redox number for 
each species is defined according to equation ( f22|) . In the redox balance, all values have unit 
of molecule cm" 2 s -1 , and defined as positive for hydrogen flux into the atmosphere. The 
redox budget for our atmosphere models is balanced, meaning that the atmosphere is not 
becoming more oxidized or reduced. 



Chemical C0 9 -dominated atmospheres 

Z 1 






species $i/ 2 (H 2 ) = 3 x 10 10 cm s -1 


$ 1/2 (H 2 ) = 3 x 10 9 cm s" 1 


No H 2 Emission 


(InliiTYiTi-niiPvn npd ttiititi n vnti n 






0„ 6 4E-7 


3.8E-6 


1.3E-3 


3 7.0E-11 


3.7E-10 


1.3E-7 


Redox balance 






Atmospheric escape 






H -7.0E+8 


-2.4E+8 


-1.2E+6 


H 2 -5.9E+10 


-5.8E+9 


-2.0E+6 


Surface emission 






H 2 6.0E+10 


6.0E+9 





CH 4 2.4E+9 


2.4E+9 





H 2 S 1.8E+9 


1.8E+9 


1.8E+9 


Dry and wet deposition 






O3 


6.3E+4 


1.1E+10 


H0 2 3.1E+4 


7.8E+7 


1.2E+8 


H 2 2 1.7E+5 


1.5E+9 


2.7E+9 


CO -3.7E+9 


-5.6E+9 


-1.5E+10 


CH 2 -1.3E+6 


-7.2E+4 





OrganicHaze -2.0E+6 


-8.1E+1 





H 2 S -4.9E+8 


-2.5E+8 


-3.2E+8 


H 2 S0 4 1.2E+8 


6.3E+7 


2.9E+7 


Balance 4.9E+3 


-9.5E+2 


-1.3E+3 
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temperature profiles are consistent with significant greenhouse effects in the convective layer 
and no additional heating above the convective layer for habitable exoplanets. The results 
discussed above do not change significantly if these temperature profiles are changed by 
several tens of K. 

Second, we have used eddy diffusion coefficients empirically determined on Earth and 
scaled the values to account for different mean molecular masses. The eddy diffusion co- 
efficient at a certain pressure level is assumed to be that of current Earth's atmosphere at 
the same pressure level (see Figured]), scaled by 6.3, 1.0, and 0.68 for the H 2 , N 2 , and C0 2 
dominated atmospheres, respectively to account for different dominant molecules. We have 
roughly scaled the eddy diffusion coefficient assuming K oc Ho where Ho is the atmospheric 
scale height. The reasoning of the scaling is as follows. According to the mixing length the- 
ory, K oc Iw, where I is the typical mixing length and w is the mean vertical velocity. The 
mixing length is a fraction of the pressure scale height (e.g., Smith 1998). The mean vertical 
velocity is related to the vertical convective energy flux, as F oc pw, where p is the pressure 
(e.g., Lindzen 1990). For a certain planet F should have the same order of magnitude for 
different atmosphere compositions. As we apply the scaling from pressure surface to pressure 
surface, we have roughly K oc H . The scaling is an approximation and we only intend to 
provide a consistent description of eddy diffusion for atmospheres with very different mean 
molecular mass. The eddy diffusion coefficients for the three scenarios are also consistent 
with their temperature profiles, featuring minima near the tropopause. The general results 
discussed above are not sensitive to the variation of eddy diffusion coefficients by an order 
of magnitude. 

Third, the typical deposition velocity on Earth is sometimes not directly applicable for 
terrestrial exoplanets, because the major surface sink of a number of gases (notably, H 2 , CH 4 , 
CO) on Earth is actually microorganisms. In this paper we focus on the scenarios assuming 
no biotic contribution in neither surface emission nor surface deposition. For H 2 and CH 4 , a 
sensible dry deposition velocity without biotic surface sink is zero (J. Kasting, 2012, private 
communication). For C0 2 , the dry deposition velocity should match with the timescale of 
weathering and carbonate formation, which is 10,000 years (Archer et al. 2010). For CO, 
the rate limiting step for converting CO into bicarbonate has been proposed to the hydration 
of CO in the ocean, which corresponds to a deposition velocity of 1.0 x 10~ 9 ~ 10~ 8 cm s~ x 
(Kharecha et al. 2005). The adopted or assumed values of dry deposition velocities are 
tabulated in Table [31 

For completeness, we now comment on our assumptions for nitrogen chemistry and 
organic haze. We do not track the bulk nitrogen cycle, instead assuming 10% N 2 for H 2 and 
C0 2 dominated atmospheres. Our model does not treat abiotic nitrogen fixation by lightning 
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(see Kasting & Walker 1981; Zahnle 1986; Tian et al. 2011 for specific analysis), so N 2 is 
considered as inert in our models. We include the formation of elemental sulfur aerosols and 
sulfuric acid aerosols in our models; but we do not include the less-understood hydrocarbon 
chemical network for hydrocarbon molecules that have more than 2 carbon atoms or the 
formation of hydrocarbon haze. Organic haze may be formed in anoxic atmospheres based on 
methane photolysis and hydrocarbon polymerization. The number density of hydrocarbons 
that have more than 4 carbon atoms (usually considered to be condensable) is ~ 4 orders 
of magnitude smaller than the number density of C 2 H 6 at the steady state (e.g. Allen 
et al. 1980; Yung et al. 1984; Pavlov et al. 2001). To account for the loss of carbon 
due to possible organic haze formation and deposition, we apply an ad hoc dry deposition 
velocity of 1.0 x 10~ 5 cm s" 1 for C 2 H 6 . This small velocity results from a scaling based on 
typical sub-micron particle deposition velocity (~ 0.2 cm s _1 , Sehmel 1980) and the number 
density ratio of C 4 hydrocarbons and C 2 H 6 (e.g., Yung et al. 1984). An accurate treatment 
for organic haze is not the purpose of this paper, because the efficiency of haze formation 
depends on oxidation states of atmospheres, stellar UV radiation, and a number of less- 
understood reaction rates. Nonetheless, for all simulated scenarios, the carbon loss due to 
haze formation is less than 1% of the methane emission flux; therefore we do not expect our 
simplification of organic haze formation would impact our results. We will fully explore the 
formation of organic haze and its implication on atmospheric chemistry in a separate paper 
(Hu et al., in prep.). 



5. Summary 

We have developed a comprehensive photochemistry model for the study of terrestrial 
exoplanet atmospheres. The photochemistry model solves the one-dimensional chemical- 
transport equation for 111 O, H, C, N, S species including S 8 and H 2 S0 4 aerosols. The 
output is the steady state of molecular mixing ratios in which concentrations of all species 
at all altitudes do not vary. In order to find the steady-state solution from arbitrary initial 
conditions for a wide variety of atmospheric compositions, required for the study of exoplan- 
ets, we have designed a numerical scheme that allows the selection of chemical species to be 
treated in fast or slow reactions automatically. The steady-state solution depends on a pool 
of input parameters, among which the major chemical species, surface emission, deposition 
velocities of long-lived species, and ultraviolet radiation are found to be critical. We validate 
the photochemistry model by simulating the atmospheric composition of current Earth and 
Mars. 

Based on the photochemistry model, we have investigated the main chemistry processes 
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and lifetimes of key spectrally active species for rocky exoplanet atmospheres by simulating 
benchmark cases of atmospheres having redox states ranging from reducing to oxidizing. 
We find that atomic hydrogen is a more abundant reactive radical than hydroxyl radical 
in anoxic atmospheres, and therefore reactions with atomic hydrogen are likely to be an 
important removal pathway for spectrally important trace gases. The source of H and OH 
is water vapor photolysis in anoxic atmospheres, and the abundance of H in the atmosphere 
is always larger than the amount of OH because OH can react with H 2 or CO to produce 
H. In addition to atomic hydrogen, in weakly oxidizing N 2 atmospheres, OH, despite its 
lower abundance than H, is important in removing CH 4 and CO. In highly oxidizing C0 2 
atmospheres, atomic oxygen is the most abundant reactive species. 

As a general observation we find that volcanic carbon compounds are long-lived and 
volcanic sulfur compounds are short-lived. In particular, due to the scarcity of OH in anoxic 
atmospheres, methane is always long-lived having chemical lifetime longer than 10,000 years. 
We also find that the reduction of C0 2 to CO and formaldehyde is minimal in N 2 atmospheres 
and limited in H 2 dominated atmospheres. In contrast to carbon species, volcanic sulfur 
compounds (i.e., H 2 S and S0 2 ) are readily converted into either elemental sulfur or sulfuric 
acid aerosols in atmospheres from reducing to oxidizing. We will discuss in detail the sulfur 
chemistry in an accompanying paper (Hu et al. 2012). 

The photochemistry is critical for prospecting the possible atmospheric composition 
that will eventually be characterized by a TPF-like mission. We have shown that volcanic 
carbon compounds including CH 4 and C0 2 are likely to be abundant in terrestrial exoplanet 
atmospheres; and in the accompany paper we will show that an enhanced volcanic activity 
leads to formation of optically thick sulfur or sulfate aerosols. As for biosignatures, we 
here have shown that photochemically produced 2 and 3 can be a potential false positive 
biosignature in thick C0 2 atmospheres. We also find that oxygen and ozone can only build 
up without H 2 and CH 4 emission; so we confirm that simultaneous detection of ozone and 
methane remains a rigorous biosignature. More generally, the three benchmark models 
presented in this paper can serve as the standard atmospheres for reducing, weakly oxidizing, 
and highly oxidizing atmospheres on habitable exoplanets for assessing chemical lifetime of 
potential biosignature gases. 

We thank James Kasting for helpful suggestions about the photochemical model. We 
thank Kerry Emanuel for enlightening discussion about modeling the Earth's hydrological 
cycle. We thank Linda Elkins-Tanton for helpful suggestions on the mantle degassing. We 
thank Susan Solomon for discussions about C0 2 cycle on terrestrial exoplanets. We thank 
the anonymous referee for the improvement of the manuscript. RH is supported by the 
NASA Earth and Space Science Fellowship (NESSF/NNX11AP47H). 
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A. Formulation of Vertical Diffusion Flux 



The vertical diffusion flux can be derived rigorously from the general diffusion equation 
for a minor constituent in a heterogenous atmosphere that 
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in which the first term is the eddy diffusion flux ($>k) and the second term is the molecular 
diffusion flux ($_d)- Note that n = Nf, and then the eddy diffusion flux can be simplified as 
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in which we have used the definition of atmospheric scale height and the ideal gas law. 

The fundamental difference between the molecular diffusion and the eddy diffusion is in 
the scale height term. The molecular diffusion depends on the specific scale height for the 
molecule, whereas the total atmospheric pressure falls according to the mean scale height. 
The molecular diffusion term can be similarly simplified as 
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which yields equation 



B. Mean Stellar Zenith Angle for 1-Dimensional Photochemistry Model 

All 1-dimensional photochemistry models need to assume a zenith angle (#o) for the 
incoming stellar radiation as a global average. Various values have been adopted in previous 
photochemistry models, for example #o = 50° (Zahnle et al. 2006), #o = 57.3° (Zahnle et al. 
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2008), and 8q = 48° (Moses et al. 2011). We find that all these assumptions are plausible 
and provide justification as follows. 

At different location of star-facing hemisphere of the planet, the local stellar zenith angle 
is /j, = cos 6 = cos if) cos where ip an d 4> are the local latitude and longitude, respectively. A 
dayside disk average should be weighted by the radiation intensity at certain optical depth 
r, viz. 

exp(— r//i ) = — / / exp(— rj cos ip cos <fi) cos 2 ?/' cos (j)dipd(f). (Bl) 

n J-tt/2 J-tt/2 

Equation ( 1B1I) is well-defined for any particular level of optical depth, and can be solved 
numerically. The relationship between the average zenith angle and the optical depth of 
concern is illustrated in Figure [7J 

We see in Figure UJ that the appropriate mean zenith angle depends on the optical depth 
of concern. We find that optical depth between 0.1 and 1.0 corresponds to a mean zenith 
angle between 57° and 48°. In the extreme of zero optical depth, the appropriate mean zenith 
angle is 60°. In general, it is appropriate to assume the mean zenith angle to be 48° ~ 60° 
for the application of 1-dimensional photochemistry model. 



C. Deposition Velocities 

In the photochemical model we calculate the number density in the bottom layer of the 
atmosphere. The interaction between the bottom layer and the surface consists of two steps: 
first, the molecular transport across a thin stagnant layer of air adjacent to the surface, 
called the quasi-laminar sublayer, second, the uptake at the surface. Each step provides a 
resistance to the overall dry deposition, viz. 

*dep = r b + r c (CI) 

in which is the quasi-laminar resistance and r c is the surface resistance. The quasi-laminar 
resistance is 

5 Sc 2 / 3 

n = , C2 

in which Sc is the dimensionless Schmidt number defined as the ratio between the kinetic 
viscosity of air and the molecular diffusivity of the molecule considered, and is the friction 
velocity (Seinfeld & Pandis 2006). The friction velocity depends on the wind speed adjacent 
to the surface and the roughness of the surface. In the current model, the friction velocity 
should be treated as a free parameter varying from 0.1 to 1 m s _1 . 
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The surface resistance is more complicated and largely depends on the property of the 
surface and the solubility of the molecule. For example, we consider two general types of 
surfaces, ocean and land. For the land we envisage desert-like surface and do not consider 
the complications of foliage. 

For the deposition to the ocean, the surface resistance is 



in which ko is the gas phase mass transfer coefficient, ki is the liquid phase mass transfer 
coefficient and H is the dimensionless Henry's law constant (Seinfeld & Pandis 2006). These 
parameters depend on the wind speed (or the friction velocity), and the Schmidt number of 
the molecule in sea water. 

The surface resistance to the land is even more complicated and depends on the prop- 
erties of land, for example, surface morphology, roughness, vegetation, canopy, etc. There 
have been tremendous efforts to measure and parameterize the surface resistance of differ- 
ent molecules on different types of land on Earth. For example, the surface resistance to a 
featureless desert can be expressed as 



in which H is in the unit of M atm 1 , fo is a normalized (0 to 1) reactivity factor, r$ = 1000 
s m -1 and tq = 400 s m _1 (Seinfeld & Pandis, 2006). 

We present the deposition velocities of slightly soluble gases (such as H 2 S), highly soluble 
gases (such as S0 2 ), not soluble but reactive gases (such as 3 ) to the ocean and the land 
in Figure [HJ The deposition velocity to the ocean critically depends on the solubility of the 
gas. For poorly soluble gas, the deposition velocity to the ocean is negligible. For the highly 
soluble gas, the deposition velocity is limited by the friction velocity, or the wind speed. For 
the reactive gas, the deposition velocity can be also very high. In particular, the deposition 
velocity of H 2 S to the land the negligible, and in most cases, the deposition velocity of H 2 S 
to the ocean is smaller than 0.1 cm s _1 . The estimate of dry deposition velocity here only 
considers the surface uptake, but not the eventual loss of a gas at or beneath the surface. 
The effective dry deposition velocity might be smaller than what we estimate if no effective 
loss mechanism is available at the surface and the surface is saturated (e.g., Kharecha et al. 




(C4) 



2005). 
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D. Photochemical Stability of Atmospheric Trace Gases 

We here present an analytical treatment of the photochemical stability of radiative trace 
gases. The goal of the analytical treatment is to obtain a useful formula that computes the 
critical pressure level of any radiative trace gas above which the photodissociation of the gas 
is important. The key idea is that the photolysis rate at any certain altitude (denoted as z*) 
is proportional to the radiation flux in the dissociating wavelengths at this altitude, which 
is attenuated by absorptions above z*. The problem can be significantly simplified, if the 
absorption is mainly due to the gas itself. This is an analog to the concept of the Chapman 
layer in the ionosphere (e.g., Banks & Kockarts 1973). 

In the following we first derive the analytical formula for the case of self-shielding and 
then extend the formula to the generic situation. 

We intend to relate the ultraviolet optical depth to the pressure level. The ultraviolet 
optical depth in the dissociating wavelength range of a specified gas is 



CO 



r uv = — — / n(z) dz, (Dl) 
Mo Jz* 

where <7uv is the characteristic cross section in the wavelengths that lead to photolysis. 
Under hydrostatic conditions, 

poo 

Patm = M atm g / N(z) dz, (D2) 

J z* 

where P atm is the atmospheric pressure at the altitude of z* and M atm is the mean molecular 
mass. We can define the column averaged mixing ratio as 

f°? n(z) dz 

and then the UV optical depth at z* can be related to the atmospheric pressure as 

Tuv = — n — • ( D4 

The molecule is subject to rapid photodissociation if the UV optical depth is smaller than 
unity, or ruv < 1, so that the critical pressure level above which photolysis is important is 

p * = /ioffM atm 

The equation ( 1D5I) can be readily extended to include other gases in the atmosphere may 
provide additional shielding, viz., 

P* = ^ Matm , (D6) 
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where a{jy and F' is the cross section and the mixing ratio of the interfering species. 

Equations (ID5[ ID6[) define a critical pressure level (or altitude) for each radiative trace 
gas subject to photolysis. Above the critical altitude, the gas is readily photodissociated and 
its mixing ratio decreases with altitude rapidly if no efficient reformation pathway exists; 
below the critical altitude, photolysis is not important and the gas is likely to be well-mixed. 
The critical pressure (altitude) is also defined for a certain mixing ratio. If the mixing 
ratio is larger, the critical pressure is smaller, and the gas is photochemically stable up to 
a higher altitude. Mean cross sections in dissociating wavelengths and the critical pressure 
level corresponding to 1 ppm mixing ratio of common spectrally active gases are tabulated 
in Table EJ 

Table 6: Mean cross sections in dissociating wavelengths of common atmospheric trace gases 
and the critical pressures of photolysis. The mean cross sections are weighted by the So- 
lar spectrum and quantum yields in the dissociation wavelengths. The critical pressure 
is evaluated according to equation (1D5I) with the conditions of /i = 0.5, g = 9.8 m s _1 , 
Mitm = 4.65 x 10~ 26 kg (N 2 -dominated atmosphere), and F = 10~ 6 (1 ppm). The additional 
shielding from other gases and reforming reactions are omitted in the estimate of critical 
pressures. 

Name Dissociation Wavelength Mean Cross Section 1-ppm Critical Pressure 





(nm) 


(cm 2 ) 


(Pa) 


H 2 


18 - 85 


1.10E-17 


2.07E+02 


2 


4 - 245 


7.18E-21 


3.17E+05 


H 2 


6 -198 


1.72E-18 


1.33E+03 


NH 3 


8 - 230 


1.31E-18 


1.73E+03 


N 2 


130 - 240 


1.38E-20 


1.65E+05 


CH 4 


52 - 152 


1.81E-17 


1.26E+02 


C0 2 


35 - 205 


7.44E-22 


3.06E+06 


H 2 S 


5 - 259 


1.11E-18 


2.05E+03 


OCS 


186 - 296 


2.74E-20 


8.33E+04 


S0 2 


5 - 395 


3.74E-18 


6.09E+02 



Despite simplicity, equations (1D51 ID6I) predict the behavior of trace gases in the upper 
part of terrestrial exoplanet atmospheres, in agreement with results from the full photochem- 
istry code. For example, 10-ppm methane in the N 2 atmosphere is photodissociated above 
an altitude of ~ 10 Pa pressure level, in agreement with Table |6j Using equation flD6j) we 
estimate the 10-ppm critical pressure for methane in the C0 2 atmosphere to be about 0.1 
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Pa, consistent with results from the full photochemistry code as shown in Figure HI For 
another example, for H 2 S in the H 2 atmosphere, the near-surface mixing ratio is about 10~ 9 , 
and according to Table |H1 we find the critical pressure is in the order of 10 5 Pa, i.e., H 2 S is 
photodissociated at all altitudes of the H 2 atmosphere. This is again consistent with the full 
photochemistry model (see the upper panel of Figure 0J. 

Equations (1D51 ID6[) provide an order-of-magnitude assessment on whether or not a 
vertically well mixed distribution of gas is a good assumption when investigating spectra of 
super Earths. When the critical pressure defined by equations (1D51 ID6[) is smaller than the 
pressure level that generates the spectral feature, it is plausible to assume the gas is vertically 
well-mixed. However, if the critical pressure is larger than the pressure level that generates 
the spectral feature, cautions should be taken because the gas could be photochemically 
depleted. If data suggests an abundance of a certain gas above its critical altitude, an 
additional mechanism, such as ultraviolet shielding by other gases, or efficient reformation, 
must be at play. In summary, equations (1D51 ID6j) provide a simplified approach to take 
into account the effects of photochemistry for the interpretation of spectral features and for 
observation planning. 
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Fig. 1. — Validation of the transport and condensation schemes of our photochemistry code. 
Only vertical transport via eddy diffusion, molecular diffusion (for H 2 only), dry deposition, 
wet deposition, condensation of water vapor, diffusion-limited escape of H 2 , and photolysis 
of C0 2 , H 2 0, H 2 , S0 2 and H 2 S are considered for an Earth-like atmosphere. We omit all 
chemical kinetics to exclusively test the transport-related schemes. H 2 S and S0 2 are the only 
sulfur compounds in this test and they are considered as removed once photo-dissociated or 
deposited. The atmosphere has the temperature-pressure profile of COSPAR International 
Reference Atmosphere (CIRA) 1986 at the equator in January, as shown by the dashed line 
on the middle panel, and major constituents of 20% 2 and 80% N 2 . Also on the middle 
panel we plot the temperature profile of US Standard Atmosphere 1976 for comparison. We 
adopt the eddy diffusion coefficient empirically derived from vertical profiles of trace gases 
(Massie & Hunten, 1981), shown in the lower panel. The boundary conditions are set as 
follows: C0 2 mixing ratio at the surface 350 ppm; H 2 mixing ratio at the surface 0.01; H 2 
mixing ratio at the surface 1 ppm; S0 2 surface emission flux 9.0 x 10 9 cm -2 s 



" 1 ; H 2 S surface 



emission flux 2.0 x 10 8 cm 2 s 1 



Water vapor in the atmosphere is limited by condensation. 
Steady-state mixing ratios are shown on the left panel, and the solid and dashed lines for H 2 




Mixing Ratio Mixing Ratio Mixing Ratio 

Fig. 2. — Profiles of key molecules in Earth's atmosphere predicted by our photochemistry 
models using the US Standard Atmosphere 1976 for mid-latitudes (Model A; dashed lines)and 
the reference tropical temperature profile (Model B; solid lines), compared with selected 
observations (horizontal lines). We adopt the eddy diffusion coefficient empirically derived 
from vertical profiles of trace gases (Massie & Hunten, 1981), as shown in Figure [U The 
error bars of observations are uniformly set to 1 order of magnitude to account for diurnal 
and spatial variations. The terrestrial data are: (1) globally averaged mixing ratios of 3 , 
N 2 and CH 4 compiled by Massie & Hunten (1981); (2) number density of H 2 from the 
US Standard Atmosphere 1976 for mid-latitudes; (3) number densities of OH and H0 2 
from balloon observations in Fairbanks, AK in 1997 (blue lines; Jucks et al. 1998); (4) 
number densities of OH and H0 2 from AURA satellite observations using the Microwave 
Limb Sounder that are zonally averaged over a latitude interval of 20° for a period of 15 
days (Pickett et al. 2006); (5) mixing ratios of NO, N0 2 and HN0 3 from balloon observations 
at Lat~35°N in 1993 (Sen et al. 1998). 
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Fig. 3. — Profiles of key molecules in Mars' atmosphere predicted by our photochemistry 
model. The chemical-transport model is computed with a 2-km grid from to 120 km. For 
the UV flux computation, unattenuated solar flux at 1.524 AU and surface albedo of 0.1 are 
used. UV cross sections are computed at 200 K when their temperature dependencies are 
available from laboratory measurements. All C, H, O, N species and relevant reactions are 
computed. Ter-molecular reaction rates are multiplied by 2.5 to account for C0 2 being the 
third body (e.g., Zahnle et al. 2008). The eddy diffusion coefficient is assumed to have a 
profile as equation ([6]) with Kt = 10 6 cm 2 s _1 , Kh = 10 8 cm 2 s -1 and zt = 20 km. The 
atmospheric temperature profile is assumed as Zahnle et al. (2008): surface temperature 210 
K, lapse rate of 1.4 K km -1 untill 50 km and then isothermal of 140 km up to 120 km. The 
relative humidity at the surface is fixed at 0.19, and the water vapor is transported up into 
atmosphere by eddy mixing, limited by condensation. The maximum saturation of water 
vapor of 60% is imposed, which gives a column number density of 9.8 precipitable microns. 
Wet deposition is reduced compared to Earth, with the rainout factor assumed to be 0.01 in 
the model. 
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Fig. 4. — Compositions of the three scenarios of rocky exoplanet atmospheres tabulated in 
Table HI The left column shows mixing ratios of H and O species, and the right column shows 
mixing ratios of N, C, and S species. From top to bottom, the three panels correspond to the 
reducing (H 2 -dominated), oxidized (N 2 -dominated), and highly oxidized (C0 2 -dominated) 
atmospheres. The vertical scales are expressed in pressure, which allows comparison between 
different scenarios that have very different mean molecular masses, and hence the altitude 
difference for a given pressure change. We highlight the profiles of three reactive species, H, 
OH, and O by thick lines. 
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Fig. 5. — Schematic illustrations of key non-equilibrium processes in the 3 scenarios of rocky 
exoplanet atmospheres considered in this paper (Table HJ), in comparison with the current 
Earth. From top to bottom, the four panels correspond to the H 2 , N 2 , C0 2 atmospheres, and 
the atmosphere of Earth. The red color highlights the reactive radicals in each atmospheric 
scenario, and the blue color highlights the major photochemical products in the atmosphere. 



-60 - 




Mixing 



Fig. 6. — Effects of the surface volcanic emission for C0 2 -dominated atmospheres on rocky 
exoplanets. The upper panel shows mixing ratios of H 2 0, CO, H 2 , and CH 4 , and the lower 
panel shows mixing ratios of 2 , O, and 3 . The solid lines show the chemical composition 
of the benchmark scenario whose parameters are tabulated in Table El In particular the 
emission rate of H 2 is 3 x 10 10 cm -2 s" 1 . The dashed lines show the chemical composition 
of the same scenario, but with an H 2 emission rate of 3 x 10 9 cm -2 s -1 ; and the dotted lines 
show the chemical composition for zero emission of H 2 and CH 4 . We see a dramatic increase 
of 2 and 3 mixing ratios as a result of a decrease of the surface emission of reduced gases. 
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Fig. 7. — Mean zenith angle of stellar radiation for photochemistry models as a function of 
optical depth. 
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Fig. 8. — Deposition velocities to the ocean (solid lines) and the land (dotted lines) as 
a function of the friction velocity u*. The calculation is for the terrestrial atmosphere at 
T = 273.15 K, P — 1 atm. The ocean is assumed to have a pH of about 6.5. The land is 
assumed to be featureless desert. 
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Table 7: Photolysis reactions in our photochemistry model. The rates (s x ) are computed 



with the unattenuated Solar spectrum at 1 AU with cross sections and quantum yields at 
295 K and 200 K, and divided by 2 for diurnal average. 



No. 


Reactants 


Products 


Rate 
295 K 


Rate 
200 K 


Cross Sections 


T° 


Quantum Yields 


T a 


Note 


1 


02 


+ 


3.52E-7 


3.52E-7 


4-108nm: Brion et al. (1979) 
108-129nm: Ogawa & Ogawa 
(1975) 


N 
N 
N 


<139nm: 0.5 
139-175nm: 
175-242nm: 1.0 


N 
N 
N 




2 


02 


+ 0(lD) 


1.40E-6 


1.40E-6 


129-173nm: Yoshino et al. (2005) 
176-203nm: Kockarts (1976) 
205-245nm:JPLRec(2011) 


N 
N 


<139nm: 0.5 
139-175nm: 1.0 
175-242nm: 


N 
N 
N 




3 


03 


02 + 


6.33E-4 


6.28E-4 


110-172nm: Mason et al. (1996) 
186-230nm:JPLRec (2011) 
230-1070nm: Bogumil et al. 
(2003) 


N 
N 
Y 


<186nm: 0.71 
186-220nm: JPL Rec (2011) 
220-340nm: Matsumi et al. (2002) 
340-411nm: 0.92 
>411nm: 1.0 


N 
N 
Y 
N 
N 




4 


03 


02 + 0(lD) 


3.66E-3 


3.63E-3 






<186nm: 0.29 
186-220nm:JPL Rec (2011) 
220-340nm: Matsumi et al. (2002) 
340-411nm: 0.08 
>411nm: 


N 
N 
Y 
N 
N 




5 


H02 


O + OH 


2.87E-4 


2.87E-4 


190-260nm:JPLRec (2011) 


N 


190-260nm: 1.0 


N 




6 


H20 


H + OH 


5.78E-6 


5.78E-6 


6-120nm: Chan etal. (1993) 
121-198nm: JPL Rec (2011) 


N 
N 


6-198nm: 1.0 


N 




7 


H202 


OH + OH 


4.77E-5 


4.57E-5 


190-350nm:JPLRec (2011) 


Y 


190-350nm: 1.0 


N 




8 


N02 


IMO + O 


4.49E-3 


4.13E-3 


16-200nm: Au & Brion (1997) 
205-235nm: Atkinson et al. (2004) 
238-667nm: Vandaele et al. 
(1998) 

667-910nm: Bogumil et al. (2003) 


N 
Y 
Y 
Y 


<398nm: 1.0 

398-428nm: JPL Rec (2011) 
>428nm: 


N 
Y 
N 




9 


N03 


NO + 02 


1.05E-2 


1.43E-2 


403-691nm: JPL Rec (2011) 


Y 


403-691nm: JPL Rec (2011) 


Y 




10 


N03 


N02 + 


8.32E-2 


1.02E-1 






403-691nm:JPL Rec (2011) 


Y 




11 


N20 


N2 + 0(1D) 


6.94E-7 


5.53E-7 


130-160nm: Nicolet & 
Peetermans (1972) 
160-240nm:JPLRec (2011) 


N 
Y 


130-140nm: 
140-230nm: 1.0 
230-240nm: 


N 
N 
N 




12 


N205 


N03 + N02 


8.55E-5 


4.82E-5 


200-410nm: JPL Rec (2011) 


Y 


<248nm: 1.0 

248-300nm: Linear interpolation 
>300nm: 


N 
N 
N 




13 


N205 


N03 + NO + 



1.51E-4 


1.45E-4 






<248nm: 

248-300nm: Linear interpolation 
>300nm: 0.7 


N 
N 
N 




14 


HN02 


NO + OH 


8.40E-4 


8.40E-4 


184-396nm:JPLRec (2011) 


N 


184-396: 1.0 


N 




15 


HN03 


N02 + 0H 


6.16E-5 


5.43E-5 


105-191nm: Suto & Lee (1984) 
192-350nm:JPLRec (2011) 


N 

Y 


<193nm: 0.33 

193-200nm: Linear interpolation 
200-350nm: 1.0 


N 
N 
N 




16 


HN03 


HN02 + 


1.94E-5 


1.92E-5 






<193nm: 0.39 

193-200nm: Linear interpolation 
200-350nm: 0.0 


N 
N 
N 




17 


HN03 


HN02 + 
0(1D) 


1.40E-5 


1.38E-5 






<193nm: 0.28 

193-200nm: Linear interpolation 
200-350nm: 0.0 


N 
N 
N 




IS 


HN04 


N02 + H02 


1.83E-4 


1.83E-4 


190-350nm:JPLRec (2011) 


N 


190-350nm:JPL Rec (2011) 


N 




19 


HN04 


N03 + OH 


4.89E-5 


4.89E-5 






190-350nm:JPL Rec (2011) 


N 




20 


CH20 


CHO + H 


4.16E-5 


4.16E-5 


226-375nm:JPLRec (2011) 


Y 


226-250nm: 0.31 
250-375nm:JPL Rec (2011) 


N 
N 




21 


CH20 


H2 + CO 


3.55E-5 


3.31E-5 






226-250nm: 0.49 
250-375nm:JPL Rec (2011) 


N 
Y 


P 

dependence 
omitted 


22 


C2H40 


CH3 + CHO 


2.59E-5 


2.59E-5 


202-360nm:JPLRec (2011) 


N 


<256nm: 0.29 
256-332nm:JPL Rec (2011) 
>332nm: 


N 
N 
N 




23 


C2H40 


CH4+ CO 


3.14E-6 


3.14E-6 






<256nm: 0.48 
256-332nm:JPL Rec (2011) 
>332nm: 


N 
N 
N 
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24 


CH402 


CH30 + 0H 


2.53E-5 


2.53E-5 


210-405nm:JPL Rec (2011) 


N 


210-405nm: 1.0 


N 




25 


CH3N02 


CH30 + NO 


1.58E-3 


1.58E-3 


190-440nm:JPL Rec (2011) 


N 


190-440nm: 1.0 


N 




26 


CH3N03 


CH30 + 
N02 


1.40E-4 


1.36E-4 


190-344nm: JPL Rec (2011) 


Y 


190-344nm: 1.0 


N 




27 


CH202 


CHO + OH 


6.20E-6 


6.20E-6 


195-250nm:JPL Rec (2011) 


N 


195-250nm: 0.7 


N 




28 


CH202 


CH02 + H 


1.33E-6 


1.33E-6 


195-250nm: 0.15 


N 




29 


CH202 


H20 + CO 


1.33E-6 


1.33E-6 


195-250nm: 0.15 


N 




30 


S03 


S02 + 


2.33E-5 


2.33E-5 


140-180nm: Hintze et al. (2003) 
180-330nm:JPL Rec (2011) 


N 
N 


140-330nm: 1.0 


N 




31 


OCS 


CO + S 


1.99E-5 


1.69E-5 


186-296nm:JPL Rec (2001) 


Y 


186-296nm: 1.0 


N 




32 


CH4 


CH3 + H 


1.95E-6 


1.95E-6 


52-125nm: Kameta (2002) 
125-142nm: Chen & Wu (2004) 
142-152nm: Lee etal. (2001) 


N 
N 
N 


0.41 (Smith & Raulin 1999) 


N 




33 


CH4 


'CH2 + H2 


2.52E-6 


2.52E-6 


0.53 (Smith & Raulin 1999) 


N 




34 


CH4 


CH + H2 + H 


2.85E-7 


2.85E-7 


0.06 (Smith & Raulin 1999) 


N 




35 


CO 


C + O 


3.94E-7 


3.94E-7 


6-177nm: Chan etal. (1993) 


N 


<lllnm 
>lllnm 


1.0 



N 
N 




36 


C02 


CO + O 


5.65E-9 


3.71E-9 


35-197nm: Huestis & Berkowitz 
(2011) 

197-270nm: Ityakov et al. (2008) 
Note: 

108-118nm: Stark et al. (2007, 
195K) 

118-163nm: Yoshino et al. (1996, 
195K) 

163-193nm: Parkinson et al. 
(2003, 195K) 


Y 
N 


<167nm 
167-205 
>205nm 




im: 1.0 



N 
N 
N 




37 


C02 


C0 + 0(1D) 


5.57E-7 


5.48E-7 


<167nm: 1.0 
167-205nm: 
>205nm: 


N 
N 
N 




38 


H2 


H + H 


1.68E-8 


1.68E-8 


18-70nm:Lee etal. (1976) 
70-87nm: Cook & Metzger (1964) 
87-124nm: Backx et al. (1976) 


N 
N 
N 


<80nm: 0.1 (Mentall & Gentieu 
1970) 

80-85nm: 1.0 
>85nm: 


N 

N 
N 




39 


N2 


N + N 


5.61E-7 


5.61E-7 


6-108nm: Chan etal. (1993) 


N 


6-108nm: 1.0 


N 




40 


CH40 


CH30+ H 


5.97E-6 


5.97E-6 


16-106nm: Burton et al. (1992) 
106-165nm: Nee et al. (1985) 
165-220nm: Cheng et al. (2002) 


N 
N 
N 


16-220nm: 1.0 


N 




41 


HCN 


H + CN 


3.16E-7 


3.16E-7 


133-145nm: Macpherson & 
Simons (1978) 


N 


133-145nm: 1.0 




Incomplete 
wavelength 
coverage 


42 


NH3 


NH2 + H 


6.89E-5 


6.89E-5 


8-106nm: Samson et al. (1987) 
106-140nm: Wu et al. (2007) 
140-230nm: Cheng et al. (2006) 


N 
N 
N 


<106nm: 0.3 (Lilly etal. 1973) 
106-165nm: Linear interpolation 
>165nm: 1.0 


N 
N 
N 




43 


NH3 


NH + H2 


1.55E-6 


1.55E-6 


<106nm 
106-165 
>165nm 


0.7 

im: Linear interpolation 



N 
N 
N 




44 


NO 


N + 


4.69E-6 


4.69E-6 


7-177nm: Lida etal. (1986) 
186-227nm: Thompson et al. 
(1968) 


N 


<191nm 
>191nm 


1.0 



N 




45 


H2S 


HS + H 


1.48E-4 


1.52E-4 


5-160nm: Feng et al. (1999) 
160-260nm: Wu & Chen (1998) 


N 

Y 


5-260nm: 1.0 


N 




46 


S02 


SO + O 


1.06E-4 


8.61E-5 


5-106nm: Feng etal. (1999) 
106-171nm: Manett & Lane 
(1993) 

171-239nm:Wu etal. (2000) 
239-395nm: Bogumil et al. (2003) 


N 
N 

Y 
Y 


<220nm: 1.0 
220-340nm: 
>340nm: 


N 
N 
N 




47 


502 


S02 1 


8.26E-4 


8.17E-4 


<220nm 
220-340 
>340nm 




im: 1.0 




N 
N 
N 




48 


S02 


S02 3 


7.51E-7 


1.66E-6 


<340nm 
>340nm 




1.0 


N 
N 




49 


HSO 


HS + O 


2.87E-4 


2.87E-4 


190-260nm:JPL Rec (2011) 


N 


190-260nm: 1.0 


N 


Assumed as 
H02 


50 


H2S04 


S03+ H20 


1.41E-5 


1.41E-5 


118-180nm: Mills et al. (2005) 
510-735nm: Lane & Kjaergaard 
(2008) 


N 
N 


118-180nm: 1.0 
510-735nm: 1.0 


N 
N 
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51 


SO 


s + o 


3.70E-4 


3.70E-4 


100-190nm: Linear interpolation 
190-235nm: Phillips (1981) 


N 
N 


100-235nm: 1.0 


N 




52 


S2 


s + s 


4.00E-3 


4.00E-3 


Photolysis around 280nm, value estimated by De Almeida & Singh (1986). 


53 


S3 


S2 + S 


5.76E-1 


5.76E-1 


350-475nm: Billmers & Smith 
(1991) 


N 


350-475nm: 1.0 


N 




54 


S4 


S2 + S2 


1.10E-1 


1.10E-1 


425-575nm: Billmers & Smith 
(1991) 


N 


425-575nm: 1.0 


N 




55 


C2H2 


C2H + H 


5.86E-6 


5.42E-6 


6-137nm: Cooper et al. (1995) 
137-185nm: Smith et al. (1991) 
185-236nm: Benilan et al. (2000) 


N 

Y 
Y 


6-217nm: 1.0 (Lauter et al. 2002) 
>217nm: 


N 




56 


C2H3 


C2H2 + H 


6.25E-3 


6.25E-3 


225-238nm: Fahretal. (1999) 
400-446nm: Hunziker et al. (1983) 


N 
N 


225-446nm: 1.0 


N 


Incomplete 
Wavelength 
Coverage 


57 


C2H4 


C2H2 + H2 


6.52E-6 


6.52E-6 


50-105nm: Holland et al. (1997) 
105-175nm:Lu et al. (2004) 
175-209nm: Orkin et al. (1997) 


N 
N 
N 


50-209nm: 0.4 


N 




58 


C2H4 


C2H2 + 2H 


8.15E-6 


8.15E-6 


50-209nm: 0.5 


N 




59 


C2H4 


C2H3 + H 


1.63E-6 


1.63E-6 


50-209nm: 0.1 


N 




60 


C2H5 


CH3 + CH2 1 


4.24E-4 


4.24E-4 


200-260nm: Fagerstrom et al. 
(1993) 


N 


200-260nm: 1.0 


N 




61 


C2H6 


C2H4 + H2 


8.84E-7 


8.84E-7 


13-120nm: Kameta et al. (1996) 
120-150nm: Chen & Wu (2004) 
150-160nm: Lee etal. (2001) 


N 
N 
N 


13-160nm: 0.14 (Akimoto et al. 
1965) 


N 




62 


C2H6 


C2H4 + 2H 


1.96E-6 


1.96E-6 


13-160nm: 0.31 


N 




63 


C2H6 


C2H2 + 2H2 


1.71E-6 


1.71E-6 


13-160nm: 0.27 


N 




64 


C2H6 


CH4 + CH2 1 


1.39E-6 


1.39E-6 


13-160nm: 0.22 


N 




65 


C2H6 


CH3 + CH3 


3.79E-7 


3.79E-7 


13-160nm: 0.06 


N 




66 


N2H2 


N2 + H2 


1.13E-4 


1.13E-4 


191-291nm: Vaghjiani (1993) 


N 


191-291nm: 1.0 


N 


Assumed as 
N2H4 


67 


N2H4 


N2H3 + H 


1.13E-4 


1.13E-4 


191-291nm: Vaghjiani (1993) 


N 


191-291nm: 1.0 


N 




68 


CH5N 


HCN + 2H2 


1.21E-5 


1.21E-5 


5-83nm: Burton et al. (1994) 
83-247nm: Hubin-Franskin et al. 
(2002) 


N 
N 


5-247nm: 0.08 


N 




69 


C2H5N 


HCN + 
CH3+ H 


2.67E-5 


2.67E-5 


lll-213nm: Basch etal. (1969) 


N 


lll-213nm: 0.1 


N 




70 


HNO 


NO + H 


8.40E-4 


8.40E-4 


184-396nm:JPL Rec (2011) 


N 


184-396: 1.0 


N 


Assumed as 
HN02 



"Temperature dependence of cross sections and quantum yields: Y indicates temperature dependence is taken 
into account for the respective wavelength range. 



Tabic 10:: Reaction rates of bi-molccular reactions (R), ter-molecular reac- 
tions (M) and thermo-dissociation reactions (T) in the photochemistry model. 
In the column "Temperature" we note the temperature range for which each 
reaction rate is valid. We annotate "Low Temperature" and "High Tempera- 
ture" to the reaction rates only valid at temperatures lower than 300 K and 
higher than 1000 K, respectively. For low-temperature applications, including 
explorations of Solar-System terrestrial planets and habitable cxoplancts, we 
exclude those reactions marked as "High Temperature". 



No. 



Reactants 



Products 



Reaction Rate 



Reference 


Temperature 


Note 


NIST 


1000-4000 




NIST 


5000-8000 




NIST 


1520-2540 




NIST 


298 




NIST 


2000-5000 




NIST 


15-295 


Low Temperature 


NIST 


298 




NIST 


298 




NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 


1 


NIST 


300-2500 


1 


NIST 


300-2500 


NIST 


300-2500 




NIST 


290-310 


Low Temperature 


NIST 


294 




NIST 


480-520 


Low Temperature 


NIST 


300-2500 




NIST 


294 




NIST 


298 




NIST 


298 




NIST 


357-676 


Low Temperature 


NIST 


293 




NIST 


293 




JPL 


200-300 


Low Temperature 


NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


298 




NIST 


298 




NIST 


297-670 


Low Temperature 


NIST 


298 




NIST 


300-2500 





Rl 

R2 

R3 

R4 

R5 

R6 

R7 

R8 

R9 

R10 

Rll 

R12 

R13 

R14 

R15 

R16 

R17 

R18 

R19 

R20 

R21 

R22 

R23 

R24 

R25 

R26 

R27 

R28 

R29 

R30 

R31 

R32 

R33 

R34 



C + CH 2 
C + CN 
C + H 2 
C + H 2 S 
C + N 2 

c + o 2 
c + ocs 
c + so 2 

C 2 H + C 2 H 3 
H 2 + C 2 H 4 
C 2 H + C 2 H 5 



+ C 2 H 5 



C 2 H 2 



C 2 H 4 



■ CN 



- CN 



C 2 H 5 + 
C 2 H 5 + 
C 2 H 5 Q 



■ C 2 H 4 
-CNO 
-CNO 
O, 



CH 
CH 
CH 
CH 
CH, 



+ C 2 H s O 
+ CN 

C 2 H 6 

CH 

CH 2 

CH. 



CH- 
■C 2 4 
CH- 
CH- 
■ CN- 

co- 
co- 
co- 



C 2 H 5 



-CH 

N 

•H 

-HS 

-N 

-O 

-CS 

-so 

H 



HCN 



HCN - 



HC 2 HO 
CH 3 Q 



-C 2 H 
r C 2 H 3 



C 2 H 6 + C 2 H 2 
C 2 H 6 + C 2 H 3 



HO, 



C 2 H 4 + C 2 H 5 
HCN + C 2 H 5 
C 2 H 4 + CH 3 
C 2 H 2 
C 2 H 3 
C 2 H 4 + H 



2.69 x 10 
4.98 x 10 
6.64 x 10 
2.1 x 10- 
8.7 x 10" 
5.1 x 10" 

1.01 x 10 
6.91 x 10 
1.6 x 10" 
1.69 x 10 
3.0 x 10- 
6.0 x 10" 
4.0 x 10- 

2.0 x 10" 

4.5 x 10" 

1.6 x 10" 

2.2 x 10" 

1.35 x 10 
5.83 x 10 

2.1 x 10- 
4.4 x 10" 
2.4 x 10" 
2.1 x 10- 

8.7 x 10- 
1.1 x 10- 

6.3 x 10" 
1.46 x 10 
1.91 x 10 

1.36 x 10 
1.3 x 10" 
2.0 x 10" 
1.57 x 10 
1.7 x 10- 
3.0 x 10" 



- 12 cxp(-23573.0/T) 

- 10 exp (-18041. 0/T) 

- 10 exp(-11700.0/T) 
10 

11 cxp (-22611. 0/T) 

n (T/298.0)-°- 3 

-10 

-11 

12 
-11 
12 
12 
11 
12 



cxp (-34277.6/T) 



13 cxp(-11800.0/T) 

12 cxp (-2300.0/T) 
10 

- 13 cxp(-1850.0/T) 

-14( T / 298 . )3.13 ex p (-9060.0/T) 

10 

11 

11 

12 cxp (-4280.0/T) 
13 



14 exp (-550.0/T) 

-i3( T /298.0) 3 - 3 exp (-5280.0/T) 

- 13 (T/298.0) 2 ' 75 cxp (-8820. 0/T) 

- 11 (T/298.0) 1 ' 26 cxp (208. 0/T) 
10 

10 

-10 

11 

11 



cxp (260. 0/T) 



Continued on next page 



No. 



Rcactants 



Products 



Reaction Rate 



TTl 

11 
11 
11 

- n exp(-3980.0/T) 
- 12 (T/298.0) 1 cxp (5800.0/T) 
12 exp(-32714.0/T) 

" 9 , 
-10 , 

n 
n 

12 
11 

- 12 cxp(-5050.0/T) 
-1B( T / 298 )3.1 cxp (-3490.0/T) 
n 

- 14 (T/298.0) 2 ' 81 cxp (-2950. 0/T) 
13 exp (-6500.0/T) 
- 14 (T/298.0) 2 ' 81 cxp (-2950.0/T) 
12 exp (-5870.0/T) 
li 



Reference 


Temperature 


Note 


NIST 


298 




NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


296-707 


Low 


NIST 


298-3000 




NIST 


300-2500 




NIST 


1100-2700 


High 


NIST 


1100-2700 


High 


NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


296-707 


Low 


NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


500-2500 




NIST 


500-2500 




NIST 


300-2500 




NIST 


1640-1950 


High 


NIST 


300-900 




NIST 


298 




NIST 


300-2500 




NIST 


300-2500 




NIST 


300-1250 




NIST 


300-2500 




NIST 


300-1500 




NIST 


300-2000 




NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


293 




NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 





R35 
R36 
R37 
R38 
R39 
R40 
R41 
R42 
R43 
R44 
R45 
R46 
R47 
R48 
R49 
R50 
R51 
R52 
R53 
R54 
R55 
R56 
R57 
R58 
R59 
R60 
R61 
R62 
R63 
R64 
R65 
R66 
R67 
R68 
R69 
R70 
R71 
R72 
R73 
R74 
R75 
R76 
R77 



CH 2 
CH 2 
CH 2 
CH 



+ C 2 H 3 



+ C 2 H 5 



CH 2 + C 2 H 6 
NH + OH 
H 2 + C 2 H 2 
CH 2 
CH 2 
CH^ + CHo 



CH 2 
CH, 



-CH3O 
-CH3O 



CH 2 
CH 2 
CH 2 + CH 3 2 
CH 2 + CH 4 
CH 2 + CH 4 
CH 2 + CHO 
CH 2 + C 2 H a 
CH 2 + C 2 H 3 



CH 2 + CH 3 2 
CH 2 + CN 
CH 2 + CNO 

CHg + C2H2 

CH 3 + C 2 H 2 
CHg + C H q 
CH 3 
CH 3 
CH 3 
CH 3 
CH 3 
CH, 



C 2 H a O 
C 2 H 3 Q 



+ C 2 H 5 



CH 3 
CH 3 
CH 3 
CH 3 
CH 3 
CH 3 
CH, 



CH 2 
CH 3 



-CHO 
CHO 
-CNO 



CH s O + C 2 H 2 
CH 3 + C 2 H 3 
CH s O + C 2 H 5 
CH3O + C 2 H 5 



C 2 H 4 + CO 
C2H2 + CHg 
C 2 H 2 + CH 3 
C 2 H 4 + CH 3 

C 2 H 5 + CH 3 



O 
+- H 
fH, 



C 2 H 2 + H 



NH 2 - 
C 2 H 3 



CH 2 - 



CH, 



3 



CH 3 + CH 
CH 3 - 
CO + CH 3 " 
C 2 H 4 + CHO 
C 2 H 4 + CHO 
C 2 H 6 + CHO 
CHO + CH 4 2 
HCN + CHO 
HCNO + CHO 
CH 4 + C 2 H 
CO + C 2 H 5 
CH 4 + C 2 H 2 
C 2 H 6 + CO 
CH 4 + C 2 H 2 
CH 4 + C 2 H 3 
CH 4 + C 2 H 3 
CH 4 + C 2 H 4 
CH 4 + C 2 H 5 
CH 4 + CHO 
CH 2 + CH 4 



CH. 



CH, 



CH 3 



CO 

+ HCN 



C 2 H 3 + CH a O 
C 2 H 4 + CH 2 Q 



2.1 x 10 - 
3.0 x 10" 
3.0 x 10" 
3.0 x 10" 
1.07 x 10 
2.94 x 10 
4.0 x 10" 
2.62 x 10 
3.32 x 10 
7.0 x 10 
4.0 x 10 
2.0 x 10 
3.0 x 10 
7.12 x 10 
1.12 x 10 
3.0 x 10- 
8.07 x 10 
3.0 x 10" 
8.19 x 10 
3.3 x 10" 
7.0 x 10" 
1.0 x 10- 
3.0 x 10- 

8.3 x 10- 

1.5 x 10- 

5.4 x 10- 
3.52 x 10 
1.57 x 10 
2.97 x 10 
1.88 x 10 
1.74 x 10 

1.6 x 10- 
4.0 x 10- 
4.0 x 10- 
1.12 x 10 
3.0 x 10" 
2.0 x 10- 

4.2 x 10- 
6.0 x 10- 
1.2 x 10- 
4.0 x 10- 
4.0 x 10- 
4.0 x 10- 



1 exp (-6010.0/T) 
1 exp (-5530.0/T) 



cxp (-8700.0/T) 
exp (385.0/T) 



11 

13 
12 
11 
11 

- 12 (T/298.0)-°' s 

- 14 (T/298.0) 3 - 7 cxp (-4780.0/T) 

- 16 (T/298.0) 5 ' 64 cxp(-1240.0/T) 

- 12 (T/298.0)-°- 5 

- 16 (T/298.0) 6 cxp(-3040.0/T) 

16 (T/298.0) 6 ' 1 exp (-990.0/T) 
11 

11 

-i5( T/ / 298 )3.i cxp (_3490.0/T) 

11 

10 

13 
11 
12 
11 
12 
11 



cxp (-4529.0/T) 



OS 
-J 

I 

High Temperature 



Continued on next page . 



No. 


Rcactants 


R78 


CH 3 


+ C 2 H 6 


R79 


CH..O 


+ CH 2 


R80 


CH3O + CH 3 2 


R81 


CH..O + 2 


R82 


CH3O 


2 + C 2 H 3 


R83 


CH3O 


2 + C 2 H 5 


R84 


CH3O 


2 + C 2 H 6 


R85 


CH3O 


2 + CH 3 : 


R86 


CH3O 


2 + CH3O. 


R87 


CH3O 


2 + CH 4 


R88 


CH 3 S 


+ CO 


R89 


CH 3 S 


+ o 2 


R90 


CH3S 


+ o 3 


R91 


CH 4 + 


C 2 H 


R92 


CH 4 4 


C 2 H 2 


R93 


CH 4 + 


C 2 H 3 


R94 


CH 4 + 


C 2 H 5 


R95 


CH 4 + 


CH 3 


R96 


CH 4 + 


CH30 


R97 


CH 4 4 


CH 3 2 


R98 


CH 4 4 


CHO 


R99 


CH 4 + 


- CN 


R100 


CH 4 + 


CNO 


R101 


CH 4 + CN 


R102 


CH 4 S + CN 


R103 


CH 4 S 


+ H 


R104 


CH 4 S 


+ H 


R105 


CH 4 S 


+ NO3 


R106 


CH 4 S 


+ 


R107 


CH 4 S 


+ OH 


R108 


CHO- 


hC 2 H 


R109 


CHO- 


|-C 2 H 3 


R110 


CHO- 


hC 2 H s O 


Rill 


CHO- 


h C 2 H 4 


R112 


CHO- 


HC 2 H 5 


R113 


CHO + CH3O 


R114 


CHO + CH3O 


R115 


CHO- 


hCH 4 


R116 


CHO- 


hCN 


R117 


CHO- 


hCNO 


R118 


CHO + 2 


R119 


CN + HCNO 


R120 


CO + C 2 H 2 



Products 



Reaction Rate 



Reference 



13 cxp (-3570.0/T) 

13 cxp(-1500.0/T) 
13 

12 

11 

11 

13 

14 cxp (390.0/T) 
14 cxp (390.0/T) 
12 cxp (-6900.0/T) 
- 13 (T/298.0) 1 ' 57 exp(-3359.0/T) 
" 8 (T/298.0)" 3 - 8 cxp (-6189.2/T) 
12 cxp (360.0/T) 

12 cxp (-250.0/T) 

13 cxp (-2479.0/T) 

- 14 (T/298.0) 4 - 02 cxp (- 2479. 0/T) 

- 15 (T/298.0) 4 ' 14 cxp (-6320. 0/T) 

- n cxp(-11600.0/T) 

- 13 exp(-4450.0/T) 

13 exp (-9299.5/T) 

" 13 (T/298.0) 2 ' 85 cxp(-11300.0/T) 

- 13 (T/298.0) 2 ' 64 cxp(150.3/T) 

- n cxp(-4090.5/T) 
10 

-10 

- n exp(-840.0/T) 
~ n exp (-1309. 8/T) 
13 cxp (210.0/T) 
- 14 (T/298.0) 4 ' 7 cxp (1210.0/T) 
12 exp (360.0/T) 
10 

10 
11 

-11 
10 
10 
10 

-13 
10 
11 
12 
10 
10 



Temperature Note 



CH 4 + CHO 
CH 2 + CH 4 O a 



-C 2 H 3 



CH3O 



o 2 

0, 



-CH30 



OCS + CH 3 
CH 3 + so 2 
CH3O + so 2 
C 2 H 2 + CH 3 
C 2 H 3 + CH 3 
C 2 H 4 + CH 3 
C H ft + CHo 



CH, 



-CH, 



CH, 



HCN + CH 3 
CH 3 + HCNO 
HCN + CH3O 
HCN + CH,S 



H 2 + CH 3 S 
CH 3 S + HNO3 
CH 3 + HSO 
H 2 + CH 3 S 



CH 
CH 4 



)H,; 



HCO 
) + CO 
CO + CHO 
hCO 
+ CO 



HCN + 
HCNO 



-CH3O 
CO 
f CO 



CO + H0 2 
■ HCN + CNO 
C 2 H + CHO 



4.0 x 10" 
1.7 x 10" 

5.0 x 10" 

9.1 x 10" 
4.0 x 10" 
4.0 x 10" 
4.9 x 10" 
5.9 x 10" 
3.6 x 10" 
3.0 x 10" 
8.43 x 10 
6.23 x 10 
1.5 x 10" 
3.0 x 10" 
5.0 x 10" 
2.13 x 10 
2.51 x 10 
1.66 x 10 

2.61 x 10 
3.0 x 10" 
1.36 x 10 
5.11 x 10 

1.62 x 10 

1.2 x 10" 
2.19 x 10 
1.15 x 10 
4.82 x 10 

4.4 x 10" 
3.45 x 10 
9.9 x 10" 
1.0 x 10" 

1.5 x 10" 
1.5 x 10" 
1.33 x 10 
2.0 x 10" 
1.5 x 10" 
3.0 x 10" 
2.41 x 10 
1.0 x 10" 

6.0 x 10" 
5.2 x 10" 

1.1 x 10" 
8.0 x 10" 



exp (-523.0/T) 



(T/298.0) 2 - 9 exp (-6600.0/T) 



cxp(-53641.4/T) 



NIST 

NIST 

NIST 

JPL 

NIST 

NIST 

NIST 

JPL 

JPL 

NIST 

NIST 

NIST 

JPL 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

JPL 

NIST 

JPL 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

JPL 

NIST 

NIST 



300-2500 

300-2500 

300-2500 

200-300 

300-2500 

300-2500 

300-2500 

200-300 

200-300 

300-2500 

200-3000 

200-2000 

200-300 

300-2500 

290-310 

300-2500 

300-2500 

1950-2770 

300-2500 

300-2500 

300-2500 

290-1500 

512-1100 

294 

293 

313-454 

313-454 

200-300 

254-495 

200-300 

300-2500 

300-2500 

300-2500 

415-670 

300-2500 

300-2500 

300-2500 

300-2500 

500-2500 

500-2500 

200-400 

298 

300-2500 



Low Temperature 



Low Temperature 
Low Temperature 



Low Temperature 
Low Temperature 



00 



Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 



Low Temperature 



Low Temperature 
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No. 



Rcactants 



Products 



Reaction Rate 



10 cxp (-45583. 2/T) 

11 exp(- 30428. 9/T) 

14 

- 12 (T/298.0) 2 - 16 cxp(-13470.5/T) 

10 cxp (-760.0/T) 

16 
10 

11 cxp (-14192. 0/T) 
10 cxp(-11200.0/T) 

12 (T/298.0) 1 - 45 exp (-1398.8/T) 

- 10 (T/298.0) - 2 
n 

n 

n 

12 (T/298.0) 2 - 53 exp (-6160.0/T) 

" n exp (-3221. 0/T) 

- n cxp(-2120.4/T) 

" 10 exp(-2120.4/T) 

- 10 exp(-7600.0/T) 
n 

12 

-11 

- 11 (T/298.0) 1 ' 5 exp (-3729.6/T) 
-11 

~ 10 exp(-85.6/T) 
n cxp(900.0/T) 
n cxp(-2068.7/T) 
~ 12 (T/298.0) 1 - 62 cxp(- 
10 cxp (-7600.0/T) 

13 cxp (-956. 2/T) 
" 9 cxp (-7920.0/T) 
- 10 (T/298.0) ' 04 

10 

n 

10 

- 11 (T/298.0) 1 ' 02 cxp (- 
- 13 (T/298.0) 3 - cxp ( 
- 10 cxp(-2670.0/T) 
- 12 (T/298.0) 2 - cxp (-2270. 0/T) 
- 13 exp(-936.0/T) 
- 13 exp(-936.0/T) 

-10 

ii 



Reference 


Temperature 


Note 


NIST 


300-2500 




NIST 


1500-1900 




NIST 


300-2500 




NIST 


300-3000 




JPL 


200-300 


Low 


NIST 


298 




NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


200-2500 




NIST 


200-2000 




NIST 


300-2500 




NIST 


298 




NIST 


298 




NIST 


300-2500 




NIST 


298-2000 




NIST 


300-2000 




NIST 


300-2000 




NIST 


700-2500 


High 


NIST 


300-2000 




NIST 


300-2500 




NIST 


295 




NIST 


300-2000 




NIST 


1000-2500 


High 


NIST 


300-2000 




NIST 


300-3000 




NIST 


290-2000 




NIST 


300-1700 




NIST 


300-2500 




NIST 


223-398 


Low 


NIST 


1000-1100 


High 


NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


273-2000 




NIST 


300-2500 




NIST 


1370-1840 


High 


NIST 


300-2500 




NIST 


250-358 


Low 


NIST 


250-358 


Low ' 


NIST 


300-2500 




NIST 


500-2500 





R121 
R122 
R123 
R124 
R125 
R126 
R127 
R128 
R129 
R130 
R131 
R132 
R133 
R134 
R135 
R136 
R137 
R138 
R139 
R140 
R141 
R142 
R143 
R144 
R145 
R146 
R147 
R148 
R149 
R150 
R151 
R152 
R153 
R154 
R155 
R156 
R157 
R158 
R159 
R160 
R161 
R162 
R163 



CO + C 2 H 4 
CO + CH 3 
C0 2 + CH 2 
C0 2 + CN 
CS + O 
CS + 3 
H + C 2 H 
H + C 2 H 
H + C 2 H 2 
H + C 2 H 2 
H + C2H3 
H + CoH, 



H + C 2 H 3 
H + C 2 H 4 
H + C 2 H 4 



H + C 2 H 4 
H + N 2 
H + C 2 H 5 
H + C 2 H 5 
H + C 2 H s O 
H + C 2 H 6 
H + C 2 HO 
H + CH 
H + CH 2 
H + CH^O 
H 

H + CH 3 



CH 2 



H + CH 3 N0 2 
H + CH 3 N0 2 
H + CH 3 
H + CH 3 
H + CH 3 
H + CH 3 O a 
H + CH 3 O a 
H + CH 4 
H + CH 4 
H + CH 4 
H + CH 4 2 
H + CH 4 2 
H + CHO 
H + CNO 



C 2 H 3 + CHO 
C 2 H 2 + OH 
CH 2 + CO 
CNO + CO 
CO + S 
O, + OCS 



"2 T ^ 2 

C 2 H + H 2 
CO + CH, 



CHO 



"H 2 
CH 3 
' + H 



C 2 H s O 



CH 4 
■N 2 -t 
CH, 



-I- CHO 
OH 
+ CH 3 
^ + H 2 
O + H, 



C 2 H 5 + H 2 

CO + CH 2 
C + H 2 
CH + H 2 
CH 3 
H 2 + CHO 



CH, 



CH 3 + HN0 2 
CH 4 



CH, 



CH, 



OH 
+ H 2 
+ OH 
O, 



CH3 + H2 

CH 3 + H 2 
CH 3 + H 2 
CH 3 + H 2 
H 2 + CH 3 2 
CO + H 2 
CO + NH 



2.5 x 10" 
6.3 x 10" 
3.9 x 10- 
1.35 x 10 

2.7 x 10" 
3.0 x 10" 
3.0 x 10" 
6.0 x 10" 
1.0 x 10" 
5.0 x 10" 
2.02 x 10 
2.0 x 10" 
2.2 x 10" 
1.2 x 10" 
4.0 x 10" 
1.33 x 10 
6.64 x 10 

1.08 x 10 
1.60 x 10 
6.0 x 10" 
3.0 x 10- 
3.32 x 10 
1.23 x 10 
4.98 x 10 

1.31 x 10 
1.0 x 10- 
4.0 x 10- 
2.14 x 10 
1.0 x 10- 
2.0 x 10- 
5.28 x 10 
2.89 x 10 

1.6 x 10- 
1.0 x 10- 
1.6 x 10- 
1.17 x 10 
5.83 x 10 

3.32 x 10 
2.42 x 10 
1.21 x 10 
1.46 x 10 
1.50 x 10 

8.9 x 10- 



1090.0/T) 



-8352.9/T) 
4040.0/T) 



I 

High Temperature 
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Reactants 



Products 
> HCN 



No. 

R164 

R165 

R166 

R167 

R168 

R169 

R170 

R171 

R172 

R173 

R174 

R175 

R176 

R177 

R178 

R179 

R180 

R181 

R182 

R183 

R184 

R185 

R186 

R187 

R188 

R189 

R190 

R191 

R192 

R193 

R194 

R195 

R196 

R197 

R198 

R199 

R200 

R201 

R202 

R203 

R204 

R205 

R206 



Reaction Rate 



Reference 



- 11 (T/298.0) ' 9 cxp (-2920.2/T) 
- 10 exp(-13350.0/T) 
- 12 (T/298.0) 1 6 exp (-9720.0/T) 
- n cxp(-1800.0/r) 
- 12 cxp(-1890.0/T) 
- 12 (T/298.0) 1 66 exp (-7000.0/T) 

-12 

n cxp(-2300.0/T) 
n cxp(-500.0/T) 
- 9 (T/298.0)-°' 5 cxp(- 
- 9 (T/298.0)-°- 3 cxp(- 
" n (T/298.0) 86 exp (- 
n cxp(-3700.0/T) 
" 13 (T/298.0) 1 89 exp (• 
- 14 (T/298.0) 3 - 29 cxp(- 
12 (T/298.0) 1 - 53 cxp (- 



Tempcraturc Note 



H 


+ 


CNO 


H 


+ 


co 2 


H 


+ 


H 2 


H 


+ 


H 2 2 


H 


+ 


H 2 2 


H 


+ 


HCNO 


H 


+ 


HCNO 


H 


+ 


HCNO 


H 


+ 


HNO 


H 


+ 


HNO 


H 


+ 


HNO 


H 


+ 


HN0 2 


H 


+ 


HN0 2 


H 


+ 


HN0 2 


H 


+ 


HN0 3 


H 


+ 


HNOg 


H 


+ 


H0 2 


H 


+ 


H0 2 


H 


+ 


HO a 


H 


+ 


N 2 


H 


+ 


NH 


H 


+ 


NH 2 


H 


+ 


NH 3 


H 


+ 


NO 


H 


+ 


NO 


H 


+ 


N0 2 


H 


+ 


N0 3 


H 


+ 


o 2 


H 


+ 


O3 



>co- 



2 + c 

2 + C 2 
2 



H 
H 
H 
H 
H 
H 
H 
H 

NH 
H 2 
H 2 

H 2 + CH3O 
H 2 + CH a 2 
H 2 + CHO 



-C 2 H 

" C 2 H 2 

■ C 2 H 3 

■ C 2 H 5 

-C 2 HO 

-C 2 HO 

f NO 
-CH 
-CH 3 



+ 
OH 
H 2 + OH 
OH + H 2 
H 2 + H0 2 
H 2 + CNO 



. CN- 

> CO- 
' H 2 4 
. OH 

> NH 2 
. OH 

> H„ 4 



hH 2 
f NH 2 
- NO 
f NH 
2 + 
-HNO 
N0 2 
. H 2 + NO 
. H 2 + N0 2 
. H 2 + N0 3 
. OH + OH 
. + H 2 

Ur, + 2 

-NH 
N 

H 2 + NH 
H 2 + NH 2 
NH + O 
OH + N 
OH + NO 
OH + N0 2 
O + OH 
OH + 2 
CH- 



^2 
■ NO 
-H, 



H 

hH 
+ H 



-H 
-H 



Z 

Continued on next page 



«C 2 H- 

* C 2 H 2 

. C 2 H 4 

> C 2 H 4 

' C 2 H 6 

' C 2 H 3 
■ CH 3 + CO 
. OH + N 2 
. CH 2 + H 

> CH 4 + H 

• CH 4 + H 

■ CH 4 2 +H 
; CH 2 Q + H 



1.86 x 10 
2.51 x 10 
6.82 x 10 
1.69 x 10 
2.81 x 10 
1.91 x 10 
1.66 x 10 
5.0 x 10~ 
3.0 x 10- 
2.41 x 10 

1.05 x 10 
1.26 x 10 
2.0 x 10 
6.39 x 10 
1.39 x 10 

5.65 x 10 

7.2 x 10- 

1.6 x 10- 
6.9 x 10- 

5.03 x 10 
1.69 x 10 

1.05 x 10 
7.80 x 10 

9.3 x 10- 

3.6 x 10" 

4.0 x 10- 

9.4 x 10- 
6.73 x 10 
1.4 x 10- 
6.64 x 10 

1.1 x 10- 

8.95 x 10 
5.0 x 10- 
1.61 x 10 
4.12 x 10 
1.0 x 10- 

1.2 x 10- 
5.86 x 10 
3.75 x 10 
6.86 x 10 

9.96 x 10 
5.0 x 10- 

2.66 x 10 



9010.0/T) 

-14673.2/T) 

-2500.0/T) 

-1940.0/T) 
-3160. 0/T) 
-8250.0/T) 



(T/298.0)- 



i cxp(-18642.2/T) 



- 10 exp (-4450. 1/T) 
-i3( T /298.0) 2 - 4 exp (-4990.1/T) 
10 (T/298.0)- 



01 exp(- 
24896. 0/T) 



10 cxp ( 
10 cxp (-340. 0/T) 



35240.0/T) 



11 

- 10 (T/298.0)-°- 59 exp (-8152. 0/T) 
10 cxp (-470.0/T) 
- 10 exp(-11700.0/T) 

10 cxp (-4000.0/T) 
- 13 (T/298.0) 2 ' 57 cxp(-130.0/T) 
13 exp(-19604.4/T) 
- 13 (T/298.0) 2 ' 63 cxp (-4300. 0/T) 
-l5( T /298.0) 3 ' 6 exp (-4250.0/T) 
n cxp(-2000.0/T) 
n cxp(-2000.0/T) 
- 12 (T/298.0)-°- 5 cxp (-60. 0/T) 

10 exp(-1660.0/T) 
- 14 (T/298.0) 2 - 74 cxp(-4740.0/T) 
- 14 (T/298.0) 2 cxp(-6720.0/T) 

11 cxp (-13110. 0/T) 
- 13 (T/298.0) 2 ' cxp (-8970.0/T) 



NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

JPL 

JPL 

JPL 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

JPL 

NIST 

NIST 

JPL 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 



500-2500 

300-2500 

300-2500 

300-1000 

300-1000 

500-2500 

1800-2600 

500-2500 

298-2000 

350-3000 

550-3000 

300-2500 

300-2500 

300-2500 

300-3000 

300-3000 

298 

298 

298 

300-4000 

1500-2500 

1100-1800 

490-1780 

1300-3000 

1500-4500 

200-300 

298 

800-3500 

200-300 

1520-2540 

2580-4650 

200-2000 

300-2500 

300-2500 

300-2500 

438-761 

438-761 

300-4000 

327-397 

300-2500 

300-2500 

300-2500 

300-2500 



High Temperature 



High Temperature 
High Temperature 



Low Temperature 

High Temperature 
Low Temperature 

High Temperature 



Low Temperature 



No. 



Rcactants 



Products 



Reaction Rate 



" 10 exp (-8059. 4/T) 
11 cxp (-13110. 0/T) 
" 12 (T/298.0)°- 5 
n cxp(-7758.0/T) 
-i4( T /298.0) 2 6 exp (-3006.8/T) 

11 cxp (-14553. 0/T) 

12 cxp (-19845. 0/T) 

-14 



Reference 


Temperature 


Note 


NIST 


300-2500 




NIST 


1300-1900 


High 


NIST 


300-2500 




NIST 


1700-3000 


High 


NIST 


833-1432 


High 


NIST 


400-2200 




NIST 


760-1000 




NIST 


1000-4000 




NIST 


300-2000 




JPL 


298 




NIST 


300-2500 




NIST 


298-669 


Low ' 


NIST 


500-2500 




NIST 


500-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


330 




NIST 


295 




NIST 


190-2237 




JPL 


200-300 


Low ' 


JPL 


200-300 


Low ' 


NIST 


500-2500 




NIST 


883-953 


High 


NIST 


298-2500 




NIST 


298 




NIST 


300-2500 




NIST 


500-2500 




NIST 


500-2500 




NIST 


1200-1900 


High 


NIST 


1000-1100 


High 


NIST 


500-2500 




NIST 


500-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


900-1200 


High 



R207 
R208 
R209 
R210 
R211 
R212 
R213 
R214 
R215 
R216 
R217 
R218 
R219 
R220 
R221 
R222 
R223 
R224 
R225 
R226 
R227 
R228 
R229 
R230 
R231 
R232 
R233 
R234 
R235 
R236 
R237 
R238 
R239 
R240 
R241 
R242 
R243 
R244 
R245 
R246 
R247 
R248 
R249 



H 



■ CN 



H 2 + H0 2 



H„ 



- NH„ 



H 9 + NO, 



HSO + O, 



H a O 



- C 2 H 5 



CH 
H 2 + CN 
H 2 + CNO 

H2O2 H - (-^Hg 

H 2 2 + C 2 H 3 
H 2 2 + C 2 H 5 

H 2 2 + CHg 

H 2 2 + CH3O 
H 2 O a + CH 3 2 
H 2 2 + CHO 
H 2 S + CH 

+ 

H 2 S + H 
H 2 S + 
H 2 S + OH 
HCN + CNO 
HNO + C 2 H 5 
HNO + CH 3 
HNO + CH3O 
HNO + CHO 



HNO - 
HNO - 
HNO - 
HNO, 



CN 
CNO 
CO 
f CH, 



HN0 2 + CNO 
H0 2 + C 2 H 2 
H0 2 + C 2 H 2 
H0 2 + C„H, 
H0 2 
H0 2 + C 2 H 4 
HQ 2 + C 2 H 4 Q 



HCN - 
NH + 

'2 

H 2 + N 2 
NH, + H 



NH 3 
H + 
CH- 



HSO 



+ H 

HN0 2 
f OH 
2 + OH 
2 + 2 
+ OH 



HCN + OH 
HCNO + OH 

C2H4 + HO2 

C 2 H 4 + H0 2 

^2^6 ~^ 

CH 4 + H0 2 
CH 4 + H0 2 
CH 4 2 +H0 2 
CH 2 + H0 2 
HS + CH 2 
HS + CH 3 
H 2 + HS 
OH + HS 
H 2 + HS 
HCNO + CN 
C 2 H 6 +NO 
CH 4 + NO 
CH 4 + NO 
CH 2 + NO 
HCN + NO 
HCNO + NO 



co 2 

CH, 



-NH 
-N0 2 
HCN + N0 2 
HCNO + N0 2 
C 2 H 3 + 2 
C 2 H 2 + OH 
OH + CH 3 + CO 
OH + CH, + CO, 



OH 



C 2 H 3 Q - 



5.65 x 10 

1.66 x 10 
5.0 x 10" 

5.73 x 10 

3.5 x 10" 
6.75 x 10 
4.0 x 10" 
1.3 x 10" 

7.74 x 10 
3.0 x 10" 

2.06 x 10 
9.48 x 10 
1.3 x 10 
9.13 x 10 
2.0 x 10" 
3.0 x 10- 
1.45 x 10 
2.0 x 10" 
5.0 x 10- 
4.0 x 10" 
1.69 x 10 
2.8 x 10" 

2.0 x 10" 
3.66 x 10 

9.2 x 10" 

6.1 x 10" 
2.0 x 10" 
1.66 x 10 
1.85 x 10 
5.0 x 10" 
1.0 x 10" 
3.0 x 10" 
3.0 x 10" 

3.3 x 10" 
3.55 x 10 
2.0 x 10" 
6.0 x 10" 
5.0 x 10" 
1.0 x 10" 
5.0 x 10" 
5.0 x 10" 
1.0 x 10" 
5.0 x 10" 



(T/298.0) 2 ' 45 exp (-1131. 0/T) 



14 



(T/298.0) 3 05 exp (-376.0/T) 

-10200.0/T) 



14 (T/298.0) 1 - 44 cxp ( 
" 12 exp(380.0/T) 
n cxp(-3760.0/T) 
" 14 (T/298.0) 2,17 cxp (—3050. 0/T) 
14 cxp (300.0/T) 

13 cxp (-4140.0/T) 
" 14 cxp (-490.0/T) 

14 exp (300.0/T) 

15 cxp(-1300.0/T) 

12 cxp (-5000.0/T) 

" 13 exp(-3490.0/T) 
10 

12 

~ 12 (T/298.0) 1 - 94 cxp (-454.6/T) 
12 cxp(-1800.0/T) 
12 , 

n cxp(-4460.0/T) 

-12 

11 (T/298.0) - 76 ' 
11 

12 , 
11 

11 

12 cxp(-6190.0/T) 

" 7 cxp (-10100. 0/T) 
11 

12 

14 (T/298.0) 1 61 exp (-7129.7/T) 

14 cxp (-4000.3/T) 
11 



1 cxp (-75.0/T) 

cxp (-4460 

2 

L (T/298.0) - 76 cxp (-175.6/T) 
! cxp(-1000.0/T) 



14 cxp (-4000.3/T) 
12 cxp (-6000.4/T) 
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No. Reactants 



Products 



Reaction Rate 



T3" 

13 

12 exp (-5870.5/T) 
- 12 cxp(-9620.6/T) 

13 , 

13 exp (-6330.0/T) 
-11 

13, 
11 

12 . 
12 

14 exp (-490.0/T) 

-11 

11 
n 
n 
10 

11 

19 
12 
12 
14 
11 
13 
11 
12 
11 
11 
11 
14 



Reference 


Temperature 


Note 


NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


1200-1800 


High Temperature 


JPL 


200-300 


Low Temperature 


NIST 


300-2500 




NIST 


623-773 


High Temperature 


JPL 


200-300 


Low Temperature 


JPL 


298 




JPL 


200-300 


Low Temperature 


JPL 


298 




JPL 


200-300 


Low Temperature 


NIST 


681-723 


High Temperature 


NIST 


298 




NIST 


298 




JPL 


220-420 


Low Temperature 


NIST 


298 




Kasting (1990) 


298 




NIST 


298 




JPL 


200-300 


Low Temperature 


JPL 


298 


-J 

to 


JPL 


298 


JPL 


200-300 


Low Temperature 


JPL 


298 




Kasting (1990) 


298 




JPL 


200-300 


Low Temperature 


NIST 


298 




NIST 


298 




NIST 


298 




NIST 


298 




NIST 


990-1010 


High Temperature 


NIST 


216-584 


Low Temperature 


NIST 


200-423 


Low Temperature 


NIST 


309-409 


Low Temperature 


NIST 


300-2500 




NIST 


1100-1400 


High Temperature 


NIST 


1400-2000 


High Temperature 


NIST 


800-3000 




NIST 


1000-3000 




NIST 


300-2500 




NIST 


1000-4000 


High Temperature 


JPL 


200-300 


Low Temperature 


JPL 


200-300 


Low Temperature 



R250 H0 2 + C 2 H 5 
R251 H0 2 + C 2 H 5 
R252 H0 2 + CH 2 



► C 2 H 4 



H 9 



R254 H0 2 + CH 3 2 
R255 H0 2 + CH 4 



R258 H0 2 + NH 2 

R259 H0 2 + NO 

R260 H0 2 + N0 3 

R261 H0 2 + 3 

R262 HS + CH 3 

R263 HS + H0 2 

R264 HS + HS 

R265 HS + N0 2 

R266 HS + O 

R267 HS + CHO 

R268 HS + 2 

R269 HS + 3 

R270 HSO + N0 2 

R271 HSO + 3 

R272 HS + CH 2 

R273 HSO a + O a 

R274 HS0 3 + O 

R275 HS0 3 + 2 

R276 N + C 2 

R277 N + C 2 H 3 

R278 N + C 2 H 3 

R279 N + C 2 H 4 

R280 N + CH 

R281 N + CH 

R282 N + CH 3 

R283 N + CH 4 

R284 N + CN 

R285 N + CNO 

R286 N + CNO 

R287 N + H 2 

R288 N + HCNO 

R289 N + NH 

R290 N + NH 2 

R291 N + NO 

R292 N + NQ 2 



• C 2 H fi 



'2^2 



► CHO 



h0 2 
-H 2 
> CH 3 + 2 
- CH 4 2 +0 
> CH..O 



2 

3 W ^2^2 

HCNO + 2 



f o 2 
+ o 2 

-N0 2 

-N0 2 + 2 

-o, + o, 



>H 2 2 

1 NH 3 ' 

> OH + 
. OH 
. OH- 

> CH 4 S 

- OH + HSO 

> S + H 2 S 

> NO + HSO 
■ H + SO 

hCO 
OH 
HSO 
HSO a + NO 
HS 



>H 2 S- 

> SO + 

> o 



. H,S- 



o 2 + o 2 



CHO 
H0 2 + S0 2 
OH + SO s 
H0 2 + S0 3 
CN + C 
C 2 H 2 N + H 
C 2 H 2 + NH 
HCN + CH 3 
C + NH 
-H 
HCN 
4-HNO 



. CN 
>H 2 - 

> CH. 

> C + 

> CO 
• CN 
. OH 

> NH 



N 



' iN 2 

NO 
-NH 
-CNO 
H 

NH + NH 

+ N 2 
N 2 Q + Q 



- N. 



5.0 x 10" 

5.0 x 10~ 

3.3 x 10~ 
5.63 x 10 

4.1 x 10~ 
1.6 x 10- 
3.32 x 10 
3.5 x 10- 

3.4 x 10- 

3.5 x 10- 
3.5 x 10- 
1.0 x 10- 
1.66 x 10 
1.0 x 10- 

1.5 x 10- 
2.9 x 10- 

1.6 x 10- 
5.0 x 10- 
4.0 x 10- 
9.0 x 10- 

9.6 x 10- 
7.0 x 10- 

1.7 x 10- 
3.0 x 10- 
1.0 x 10- 
1.3 x 10- 

2.8 x 10- 

6.2 x 10- 
1.2 x 10- 

4.2 x 10- 
3.0 x 10- 
1.66 x 10 

4.3 x 10- 
4.0 x 10- 
3.0 x 10- 
3.3 x 10- 
1.66 x 10 
6.03 x 10 
3.85 x 10 
1.95 x 10 

3.0 x 10 

2.1 x 10 
5.8 x 10 



1 exp (750.0/T) 
1 exp (-6330.0/' 

Ll 

1 exp (430.0/T) 
1 exp (250.0/T) 



exp (240.0/T) 



exp (-280.0/T) 



exp (-800.0/T) 



exp (-330.0/T) 



"(T^gS.O) ' 65 exp (-1203.0/T) 
- 10 (T/298.0)-°- 09 
1() cxp(-420.0/T) 
10 exp (-4330.0/T) 

10 

n 

-12 

- 11 (T/298.0) 1 - 2 exp (-19243.6/T) 

- 5 cxp(-26460.0/T) 

- 11 (T/298.0) ' 51 cxp(-9.6/T) 

13 exp (-7600.0/T) 

n cxp(100.0/T) 

12 exp (220.0/T) 
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No. Reactants 



Products 



Reaction Rate 



Reference 



TT 

16 
■11 
■12 
■11 
■12 
■11 

- 10 exp(-13470.0/T) 
- 10 (T/298.0)- 1 03 cxp (-420.0/T) 
13 cxp(1140.2/T) 
13 cxp (1140.2/T) 

-10 

-11 

- 14 (T/298.0) - 79 exp(-601.4/T) 
-i4(2y 2 98.o)2-0 cxp (-3270.2/T) 

-11 

12 (T/298.0) 1 2 
" 13 exp(-10500.0/T) 
13 cxp (-2780.0/T) 
- 14 cxp(-1320.6/T) 
13 cxp(-1250.0/T) 
12 cxp (-4678.6/T) 
- 14 (T/298.0) 3 ' 46 cxp (-2820.0/T) 
10 cxp (-4834.9/T) 
- 15 (T/298.0) 3 exp(-2130.0/T) 
- n (T/298.0) 1 91 cxp (1039.0/T) 
13 (T/298.0) 1 9 exp (-5725.0/T) 
12 cxp (450.0/T) 

12 cxp (650.0/T) 

-11 

-11 

12 

12 exp (-930.0/T) 

- 13 (T/298.0) 1 5 cxp (-230.0/T) 

-11 . 



Temperature Note 



R293 N + O a 

R294 N + 3 

R295 N + OH 

R296 N a O + C 

R297 N 2 + CH 

R298 N a O + CN 

R299 N 2 + HCf 

R300 NH + NH 3 

R301 NH + NO 

R302 NH + NO a 

R303 NH + N0 2 

R304 NH + O 

R305 NH + O 

R306 NH + 2 

R307 NH + 2 

R308 NH + OH 

R309 NH + OH 

R310 NH 2 + C 



R313 NH 2 + C 2 H 4 

R314 NH 2 + C 2 H 5 

R315 NH 2 + C 2 H 6 

R316 NH 2 + CH 3 

R317 NH 2 + CH 4 

R318 NH 2 + CNO 



R321 NH, + NO, 



R324 NH 2 + O 

R325 NH 2 + 3 

R326 NH 2 + OH 

R327 NH 3 + CH 

R328 NH a + CH 3 

R329 NH 3 + CN 

R330 NH 3 + CNO 

R331 NO + C 

R332 NO + C 

R333 NO + C 2 H 

R334 NO + C 2 H 3 

R335 NO + C 2 H 5 Q 



NO + O 
NO + 2 
NO + H 
CO + N 2 
HCN + NO 
N 2 + CNO 
■ HCNO + N 2 
NH 2 + NH 2 
N„0 + H 
- HNO 
N 2 
H 
-N 
•NO 



NO - 
OH 
NO - 
OH 
OH 
O + HNO 
H + HNO 
H 2 + N 
CH + NH 
C„H + NH 



C 2 H 3 



3 

NH, 



3 

NH 



C 2 H 6 +NH 
C 2 H S +NH 3 
CH 4 + NH 
CH 3 + NH 
HCNO 
OH + NH 3 
H 2 + N 2 
H 2 + N 2 
H + HNO 
OH + NH 
H 2 + NO 
H 2 + N0 2 
H 2 + NH 
HCN + H 2 



H 



HCN + NH 2 
CO + N 2 H 3 
CO + N 
CN + O 
HCN + CO 
CH 2 + HCN 



1.5 x 10 
2.0 x 10 
4.7 x 10 

8.0 x 10 

3.1 x 10 
2.0 x 10 

2.0 x 10 
5.25 x 10 
1.17 x 10 

2.1 x 10" 

1.4 x 10" 
1.16 x 10 
1.16 x 10 
6.74 x 10 
6.79 x 10 
3.32 x 10 

3.1 x 10" 
9.61 x 10 

8.2 x 10" 
3.42 x 10 

3.5 x 10" 
5.1 x 10" 
2.74 x 10 
8.4 x 1G- 
8.77 x 10 
1.71 x 10 
2.1 x 10- 

4.0 x 10- 

2.1 x 10- 
7.47 x 10 
1.16 x 10 

8.3 x 10- 
4.3 x 10- 
7.69 x 10 
7.24 x 10 
9.55 x 10 
1.66 x 10 
6.29 x 10 
3.49 x 10 
5.57 x 10 

3.6 x 10 
5.0 x 10 
1.3 x 10 



exp (-3600.0/T) 



exp (257.4/T) 



exp (317.0/T) 
cxp (-4895.0/T) 



-14 
-11 

14 (T/298.0) 2 - 48 exp (-493.0/T) 
11 (T/298.0)-°- 02 
- n (T/298.0)-°- 32 



n (T/298.0)- 3 ' 38 exp (-516.0/T) 
li 



JPL 

JPL 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

JPL 

JPL 

NIST 

NIST 

NIST 

JPL 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 



200-300 
298 

300-2500 
298 

200-1250 
297 

1400-1800 

1300-1700 

300-4000 

200-300 

200-300 

250-3000 

250-3000 

298-2200 

300-3300 

298-3000 

298-3000 

1000-4000 

340-510 

250-465 

297- 453 
300-2000 
300-2000 
300-2000 
300-2000 
1000-3000 
300-3000 
200-300 
200-300 

298- 3000 
298-3000 
298-3000 
200-300 
250-3000 
300-1300 
350-600 
300-700 
295-882 
200-5000 
200-5000 
298 

200-2000 
298 



Low Temperature 



High Temperature 



Low Temperature 



Low Temperature 

Low Temperature 

Low Temperature 
-4 
GO 

I 



High Temperature 

Low Temperature 
Low Temperature 



Low Temperature 
Low Temperature 
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No. 



Rcactants 



Products 



Reaction Rate 



Reference 



11 cxp (320.0/T) 

" 12 (T/298.0)"°- 72 exp (200.0/T) 
-11 

11 

10 

-10 

10 

" 12 exp(-3007.0/T) 
" 13 cxp(-1443.0/T) 

12 cxp (-7900.0/T) 
12 (T/298.0)-° 7 



Temperature Note 



R336 NO + C 2 HO 

R337 NO + C 2 HO 

R338 NO + CH 

R339 NO + CH 

R340 NO + CH 

R341 NO + CH 

R342 NO + CH 

R343 NO + CH 2 

R344 NO + CH 2 

R345 NO + CH a 

R346 NO + CHgO 

R347 NO + CH 3 2 

R348 NO + CHO 

R349 NO + CN 

R350 NO + CNO 

R351 NO + CNO 

R352 NO + CNO 

R353 NO + N0 3 

R354 NO + 3 



R357 N0 2 + C 2 H 3 

R358 N0 2 + C 2 H 3 

R359 N0 2 + C 2 H s O 

R360 N0 2 + C 2 HO 

R361 NO, + C,HO 



R363 NO, + CH, 



R365 NO, + CH 3 



R367 NO, + CH 4 



CO + HCNO 
HCN + C0 2 
N + CHO 
H + CNO 
CO + NH 
O + HCN 
CN + OH 
HCNO + H 
HCN + OH 
HCN + H 2 
CH a O + HNO 
CH 3 + N0 2 
CO + HNO 
CO + N 2 
CO + N 2 
CO, 



CO + N, + O 



CO + NO 



C0 2 + CH 3 + NO 
C 2 H 4 + HNO a 
C0 2 + HCNO 
CO + CHO + NO 
CHO + NO 
CH 2 + NO 
CH 3 + NO 
CH 2 + HNO 
CH 2 + HN0 2 
HN0 2 + CH 3 



NO 
NO - 
NO - 
NO - 
CO- 

co 2 

CO, 



CH0 2 
CO + 

co 2 ^ 

CNO 
N 2 



OH 
-H 



CO- 
NO - 



R377 N0 2 + O s 
R378 NQ 3 + C 2 H 3 Q 



NO + NO 
-N0 2 + 2 
NO s + 2 
C 2 H 2 Q + HNQ 3 



1.4 x 10" 
6.12 x 10 
1.33 x 10 
4.0 x 10" 

2.5 x 10- 
1.37 x 10 

1.4 x 10" 
4.17 x 10 
8.32 x 10 
4.0 x 10 
4.0 x 10 
2.8 x 10 

1.2 x 10 
1.79 x 10 
5.15 x 10 
1.29 x 10 
2.78 x 10 

1.5 x 10" 

3.0 x 10- 

1.1 x 10- 

7.6 x 10" 
4.21 x 10 

7.3 x 10" 

1.3 x 10" 
9.72 x 10 
1.46 x 10 
1.45 x 10 
6.91 x 10 
2.26 x 10 

5.4 x 10" 

1.1 x 10" 
1.58 x 10 
2.71 x 10 
4.0 x 10" 
1.94 x 10 
4.0 x 10" 
7.11 x 10 

5.2 x 10" 
1.64 x 10 

1.3 x 10" 

4.5 x 10" 
1.2 x 10" 
2.5 x 10" 



12 exp (300.0/T) 
11 



- 10 exp(-4040.0/T) 
" n (T/298.0)" 1 - 34 exp (-360.0/T) 
- 10 (T/298.0)- 1 - 97 exp (-560.0/T) 
- 11 (T/298.0)- 1 - 73 exp (-380.0/T) 
11 exp (170.0/T) 
12 cxp(-1500.0/T) 
10 

n cxp(130.0/T) 

" n (T/298.0)" ' 6 

11 cxp (80.6/T) 

11 (T/298.0)" 1 - 

" 12 cxp(171.0/T) 

- n exp(m.0/T) 
-10 

-11 

-11 

12 

n cxp(-1200.0/T) 

- 9 exp(-19364.0/T) 
-11 

11 

" 10 (T/298.0)" - 75 exp (-970.0/T) 
n cxp(186.0/T) 

-12 

12 

-11 
12 

14 cxp(-1260.0/T) 

13 cxp (-2450.0/T) 
11 



NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

JPL 

NIST 

NIST 

NIST 

NIST 

NIST 

JPL 

JPL 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

JPL 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

JPL 

JPL 

NIST 



296-2000 

296-2000 

298 

298 

298 

298 

289 

1100-2600 

1100-2600 

1300-2080 

200-400 

200-300 

295-2500 

687-1500 

300-3000 

300-3000 

300-3000 

200-300 

200-300 

300 

288-800 

220-336 

298-328 

200-300 

298-423 

298-423 

298 

295 

233-1700 
298 

200-300 
600-1100 
296 
298 

294-1650 

500-2500 

298 

298 

296 

296 

200-300 
200-300 
298 



High Temperature 
High Temperature 
High Temperature 
Low Temperature 
Low Temperature 

High Temperature 



Low Temperature 
Low Temperature 



Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 



Low Temperature 



Low Temperature 
Low Temperature 
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No. 



Rcactants 



Products 



Reaction Rate 



11 
11 

12 
12 
12 

13 cxp (-2646.0/T) 

- 12 cxp(-5027.0/T) 

13 cxp (-2450.0/T) 

10 

n 

11 cxp(-2280.0/T) 

~ 12 (T/298.0) 1 ' 5 cxp (-850.0/T) 

- 11 (T/298.0) ' 64 

11 (T/298.0) ' 2 exp (215.0/T) 

11 (T/298.0) ' 2 exp (215.0/T) 

10 

- 13 (T/298.0) 1 - 88 exp (-91.4/T) 
- 12 (T/298.0) 1 ' 91 cxp(-1880.0/T) 
12 (T/298.0) 1 - 55 exp (-215.0/T) 
14 (T/298.0) 1 ' 88 exp (-91.4/T) 

12 cxp (-754.0/T) 
n cxp(-1100.0/T) 

11 (T/298.0) ' 03 exp (198.45/T) 
-10 

12 (T/298.0) ' 03 exp (198.45/T) 
-11 

- 12 (T/298.0) 1 ' 5 cxp (-2920.0/T) 

- n exp(-2380.0/T) 
n 

12 
11 
10 
11 

" n («[i]/298.0, 0.57) cxp (-1390.0/T) 
- 14 (T/298.0)- 2 ' 12 cxp (-314.0/T) 

10 

-li 

11 cxp (-2620. 0/T) 

- n exp(-2650.0/T) 
n 

n 

14 



Reference 


Temperature 


Note 


JPL 


200-300 


Low Temperature 


NIST 


298 




NIST 


298 




NIST 


298 




NIST 


298 




NIST 


298 




NIST 


258-367 


Low Temperature 


NIST 


623-773 


High Temperature 


JPL 


200-300 


Low Temperature 


NIST 


8000 


High Temperature 


NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2500 




NIST 


253-353 


Low Temperature 


NIST 


200-2500 




NIST 


200-2500 




NIST 


300-2500 




NIST 


300-2000 




NIST 


290-1510 




NIST 


300-2500 




NIST 


300-2000 


Low Temperature 


NIST 


226-380 


JPL 


200-300 


Low Temperature 


NIST 


200-2500 




NIST 


300-2500 




NIST 


200-2500 




NIST 


300-2500 




NIST 


300-1200 




NIST 


10-6000 




NIST 


300-2000 




NIST 


300-2500 




NIST 


1200-1800 


High Temperature 


NIST 


300-2500 




NIST 


300-2500 




NIST 


250-2200 




NIST 


300-2500 




NIST 


300-2500 




NIST 


290-900 




NIST 


300-410 


Low Temperature 


NIST 


294-473 


Low Temperature 


NIST 


300-2500 




NIST 


298 




JPL 


200-300 


Low Temperature 



R379 N0 3 + C 2 H 4 
R380 NO s + C 2 H 5 



R381 NO, 



CH, 



R382 NO s + CH 3 

R383 N0 3 + CH3O 

R384 N0 3 + CH 3 2 

R385 N0 3 + CH 4 



R388 O + C 2 

R389 O + C 2 H 

R390 O + C 2 H 2 

R391 O + C 2 H 2 

R392 O + C 2 H 2 r 

R393 O + C 2 H 3 

R394 O + C 2 H 3 

R395 O + C 2 H 3 

R396 O + C 2 H 4 



R399 O + C 2 H 4 

R400 O + C 2 H 4 

R401 O + C 2 H 4 

R402 O + C 2 H 5 

R403 O + C 2 H 5 

R404 O + C 2 H B 

R405 O + C 2 H 5 

R406 O + C 2 H 6 

R407 O + CH 

R408 O + CH 

R409 O + CH 2 

R410 O + CH 2 

R411 O + CH 2 

R412 O + CH 2 

R413 O + CH 2 

R414 O + CH 3 

R415 O + CH 3 

R416 O + CH 3 

R417 O + CH 3 N0 2 

R418 O + CH3NO3 

R419 O + CH3O 

R420 O + CH3O 

R421 Q 3 + C 2 H 2 



C 2 H 3 4-HN0 3 
C 2 H 5 + N0 2 
CH3O + NO a 
CH 3 2 +N0 2 



CH3O + NO a 
CH3O + HNO 



C0 2 
N0 2 
CO- 
CO- 



3 

HNO 



-I- NO 
4- N0 2 + 2 
C 

CH 



CO + CH, 



CH, 



) + H 
) + H 
f OH 
CHO 



CH 2 + CH 2 
C,H,0 + OH 



OH 4 
CO 4 
CH 4 
CHO 
CO 4 
CO 4 
CHO 



O + H 
4- OH 
) + CH 3 
4- OH 
C 
H 

OH 
+ H 
H + H 
H 2 

4- OH 



CH 2 + H 
CO + H 2 4-H 
OH4-CH 2 4 
OH4-CH 2 4 
CH 2 + OH 
CH 3 4- 2 
CO + CH 



NO 
NO, 



1.4 x 10 
4.0 x 10 

3.5 x 10 
1.8 x 10 
1.5 x 10 

1.3 x 10 

9.4 x 10 

1.66 x 10 

8.5 x 10" 
6.0 x 10- 

1.7 x 10" 
1.5 x 10" 
3.49 x 10 
8.49 x 10 
5.5 x 10" 

5.5 x 10" 

1.6 x 10" 
3.85 x 10 
1.33 x 10 
1.5 x 10- 

5.5 x 10- 
8.3 x 10- 

1.8 x 10- 
6.3 x 10- 
1.33 x 10 
6.3 x 10- 

2.67 x 10 
8.54 x 10 
2.52 x 10 

6.6 x 10- 
7.2 x 10- 
5.0 x 10- 

1.2 x 10- 

7.3 x 10- 
1.78 x 10 
7.51 x 10 

1.4 x 10- 
5.72 x 10 
2.32 x 10 
2.49 x 10 
7.0 x 10 

2.5 x 10 
1.0 x 10 



exp(-1900.0/T) 



exp (-4100/T) 



Continued on next page . 



Products 
► CH 3 



No. 



Reactants 



Reaction Rate 



Reference 



TT 

- 12 (T/298.0) 2 ' 2 cxp (-3820.0/T) 
- n cxp(-2360.0/T) 
11 « 
11 
11 
-11 

- n exp (-13711. 0/T) 

11 cxp (-210. 0/T) 
11 

- 10 (T/298.0)- 143 exp (-3500.0/T) 
- 13 (T/298.0) 2 - 67 exp(-3160.0/T) 

12 exp (-2000.0/T) 

12 cxp (-4000.0/T) 

12 exp (-4000.0/T) 

-i3( T /298.0) 2 11 cxp (-5750.0/T) 

- 13 (T/298.0) 141 cxp (-4290. 0/T) 
11 

n cxp(-3400.0/T) 

n cxp(200.0/T) 

12 « 
11 

12 cxp (-2060.0/T) 
n cxp(120.0/T) 

-11 

-10 
-10 
11 
-10 
12 
11 
10 

11 « 

10 

-10 , 
11 

" n cxp (110. 0/T) 

" n exp(20.0/T) 

- n cxp(20.0/T) 
10 

11 

10 

11 , 



Temperature Note 



R422 
R423 
R424 
R425 
R426 
R427 
R428 
R429 
R430 
R431 
R432 
R433 
R434 
R435 
R436 
R437 
R438 
R439 
R440 
R441 
R442 
R443 
R444 
R445 
R446 
R447 
R448 
R449 
R450 
R451 
R452 
R453 
R454 
R455 
R456 
R457 
R458 
R459 
R460 
R461 
R462 
R463 
R464 



O + CH 3 2 
O + CH, 



O + CHO 
O + CHO 
O + CH0 2 
O + CN 
O + CN 
O + CNO 
O + CNO 
+ H 2 
+ H 2 2 
O + HCN 
O + HCN 
O + HCNO 
O + HCNO 
O + HNO 
O + HN0 4 
+ H0 2 
+ N0 2 
O + NO3 
O + O3 
O + OH 

o + s 2 






45)- 


hC 2 U 4 





45)- 


HC 2 H 6 





45)- 


hCH 4 





45)- 


hCH 4 





45)- 


hCH 4 





45)- 


hCH 4 





45)- 


HCH 4 





45)- 


hC0 2 





45)- 


KH 2 





4)) + H 2 





45)- 


h HCNO 





45)- 


h N 2 





45)- 


hN 2 





45)- 


HN 2 





45)- 


I" NHg 





45)- 


I- NO 





45)- 


hN0 2 





45)- 


h0 2 



+ o 2 

> CH 3 + OH 

■ CH3O + OH 

• CH 3 2 +OH 
. CO + OH 

> C0 2 + H 

> C0 2 + OH 

> C + NO 
. CO + N 

- CO + NO 
. CN + O a 

. H + OH 

> OH + H0 2 
. H + CNO 

. CO + NH 
. CNO + OH 
. C0 2 + NH 

- OH + NO 

> HNO3 + o 2 
. OH + O a 

. NO + 2 

• N0 2 + 2 

- o 2 + o 2 

>0 2 + H 

■ s + so 

OH 
OH 
CH 3 6 + H 
CH 3 + OH 
CH 2 + H 2 
CH 3 2 +H 
CH3O + OH 

co 2 + o 

H + OH 
OH + OH 
C0 2 + NH 
+ N 2 
O a + N 2 
NO + NO 
OH + NH 2 
O a + N 
O a + NO 

o + o 2 



'C 2 H 3 
. C 9 H, 



8.3 x 10" 
2.26 x 10 
1.66 x 10 

3.3 x 10" 
5.0 x 10" 
5.0 x 10- 
1.44 x 10 
5.37 x 10 

3.4 x 10" 

7.5 x 10" 
4.05 x 10 
3.44 x 10 

1.4 x 10 - 
7.0 x 10- 
3.0 x 10- 
6.16 x 10 
4.93 x 10 
6.0 x 10- 
7.8 x 10- 

3.0 x 10" 

5.1 x 10- 
1.0 x 10- 
8.0 x 10- 

2.2 x 10- 
1.66 x 10 
2.19 x 10 
6.29 x 10 

3.5 x 10" 
1.31 x 10 
9.0 x 10- 
9.0 x 10- 
4.2 x 10- 

7.5 x 10- 

1.2 x 10- 
1.63 x 10 

4.6 x 10- 
2.15 x 10 
4.63 x 10 
7.25 x 10 
2.5 x 10 
8.5 x 10 
3.0 x 10 

3.3 x 10 



exp(-2390.0/T) 



1 cxp (210.0/T) 



- cxp (115.0/T) 
1 exp (60.0/T) 



cxp (55.0/T) 



NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

JPL 

NIST 

NIST 

NIST 

NIST 

NIST 

JPL 

JPL 

JPL 

NIST 

NIST 

NIST 

JPL 

JPL 

JPL 

JPL 

JPL 

JPL 

NIST 

NIST 

NIST 

JPL 

JPL 

JPL 

NIST 

NIST 

JPL 

JPL 

JPL 

NIST 

JPL 

JPL 

JPL 

JPL 

NIST 

NIST 

JPL 



300-2000 
420-1520 
300-1000 
300-1000 
300-2500 
300-2500 
298 

200-5000 

500-2500 

500-2500 

500-2500 

300-2500 

200-300 

470-900 

470-900 

500-2500 

500-2500 

300-2500 

200-300 

200-300 

200-300 

200-300 

200-300 

200-300 

298 

298 

300 

200-300 

200-300 

200-300 

300 

300 

200-300 
200-300 
200-300 
298 

200-300 

200-300 

200-300 

200-300 

300 

300 

200-300 



Low Temperature 



Low Temperature 
Low Temperature 
Low Temperature 
Low Tempcr^Jjirc 
Low Temperature 
Low Temperature 



Low Temperature 
Low Temperature 
Low Temperature 



Low Temperature 
Low Temperature 
Low Temperature 

Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 



Low Temperature 



Continued on next page 



No. 


Rcactants 


R465 


0( i D) + 3 


R466 


0( 1 D) + 3 


R467 


O^DJ + OCS 


R468 


3 + C 2 H 4 


R469 


3 + 


C 2 H 5 


R470 


O3 + 


CH 3 


R471 


3 + 


CH3O 


R472 


3 + 


CH 3 2 


R473 


ocs 


4- OH 


R474 


OCS + CH 


R475 


OCS + O 


R476 


OH 4 


-c 2 


R477 


OH 4 


-C 2 H 


R478 


OH 4 


-C 2 H 


R479 


OH 4 


- C 2 H 2 


R480 


OH4-C 2 H 2 


R481 


OH 4 


" C 2 H 2 


R482 


OH 4 


- C 2 H 2 


R483 


OH 4 


-C 2 H 2 


R484 


OH 4 


-C 2 H 2 


R485 


OH 4 


-C 2 H 3 


R486 


OH 4 


- C2H2 


R487 


OH 4 


-C 2 H 3 


R488 


OH 4 


- C 2 H 4 


R489 


OH 4 


-C 2 H 4 


R490 


OH 4 


-C 2 H 5 


R491 


OH 4 


" C 2 H 6 


R492 


OH 4 


-CH 2 


R493 


OH 4 


-CH 2 


R494 


OH 4 


-CH 2 2 


R495 


OH 4 


-CH 3 


R496 


OH 4 


-CH 3 


R497 


OH 4 


-CH 3 


R498 


OH 4 


- CH3NO3 


R499 


OH 4 


-CH..O 


R500 


OH + CH,0, 


R501 


OH 4 


-CH 4 


R502 


OH 4 


-CH 4 


R503 


OH 4 


-CH 4 2 


R504 


OH 4 


-CHO 


R505 


OH 4 


-CH0 2 


R506 


OH4-CN 


R507 


OH4-CN 



Products 



Reaction Rate 



Reference 



Temperature Note 



2 4- 2 
2 4- O + O 
CO 4- SO 
CH 2 4- CH 2 2 
C 2 H s O- 



■o 2 

O, 



o 9 



CHO 4- CS 
CO 4- SO 
CO 4- CH 
CO 4- CH 2 
C 2 H 2 4-0 
C 2 H 2 + H 
C 2 H4-H 2 
C 2 HO 4- H 2 



CO- 
CO- 
CO, 



CH 3 
CH3O 
f CH, 



C 2 H 2 4-H 2 

C2H2 + 

C 2 H 3 4-H 2 



CH 2 4- H 
CH0 4-H 2 
CH0 2 4- H 2 
H 



CH3O - 
CH 4- 



H„ 



CH,0 4- NO, 



CH 4 4- 2 



CH3O - 



C0 4-H 2 
C0 2 + H 2 ( 

H 4- CNO 
O 4- HCN 



1.2 x 10" 
1.2 x 10" 

3.0 x 10- 

1.2 x 10" 
3.32 x 10 
5.4 x 10" 
2.53 x 10 
2.9 x 10 - 

1.1 x 10- 

2.0 x 10" 

2.1 x 10- 

8.3 x 10" 
3.0 x 10- 
3.0 x 10- 
4.12 x 10 

1.03 x 10 
1.91 x 10 
6.34 x 10 
1.0 x 10- 
4.98 x 10 
5.0 x 10- 
5.0 x 10- 
2.0 x 10- 
1.66 x 10 

4.4 x 10- 
4.0 x 10- 

8.7 x 10- 
3.0 x 10" 

5.5 x 10- 
4.0 x 10- 
6.45 x 10 

1.2 x 10- 
2.59 x 10 
8.0 x 10- 
3.0 x 10- 
1.0 x 10- 
2.45 x 10 
2.9 x 10- 

3.8 x 10- 
1.69 x 10 
1.0 x 10- 
7.0 x 10- 
1.0 x 10- 



10 
10 
10 

14 cxp (-2630.0/T) 

-14 

12 cxp(-220.0/T) 

-14 

16 cxp(-1000.0/T) 

13 cxp(-1200.0/T) 

10 cxp (190.0/T) 

n cxp(-2200.0/T) 
12 

11 

11 

- 13 (T/298.0) 2 ' 3 exp (-6790.0/T) 
- 13 (T/298.0) 2 ' 68 exp (-6060.0/T) 

-13 
-18,- 
11 
-12 
11 
11 
11 



(T/298.0) 4 ' exp (1010.0/T) 



- 13 (T/298.0) 2 ' 75 cxp (-2100.0/T) 

12 cxp (364.0/T) 

11 



cxp(-1070.0/T) 
cxp (125.0/T) 



- 13 (T/298.0) 1 exp(-6012.0/T) 

1() cxp(-1400.0/T) 

- 13 (T/298.0)" - 53 exp (-5440.0/T) 

13 cxp(-1000.0/T) 
11 

10 

" 12 exp (-1775. 0/T) 
12 cxp (-345. 0/T) 

12 ( 

-10 

10 
11 
11 . 



' exp (200.0/T) 



exp(-1000.0/T) 



JPL 

JPL 

NIST 

JPL 

NIST 

JPL 

NIST 

JPL 

JPL 

NIST 

JPL 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

JPL 

NIST 

JPL 

JPL 

NIST 

NIST 

NIST 

JPL 

NIST 

NIST 

JPL 

JPL 

JPL 

NIST 

NIST 

NIST 

NIST 



200-300 
200-300 
300 

178-360 
298 

200-300 
298 

200-300 

200-300 

297-667 

200-300 

2200 

300-2500 

300-2500 

500-2500 

300-2500 

295 

500-2500 
298 

800-1200 

300-2500 

300-2500 

300-2500 

300-2500 

200-350 

300-2500 

200-300 

300-2500 

200-300 

200-300 

300-3000 

300-1000 

600-2500 

200-300 

300-2500 

300-2500 

200-300 

200-300 

200-300 

300-2500 

300 

500-2500 
500-2500 



Low Temperature 
Low Temperature 

Low Temperature 

Low Temperature 

Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
High Temperature 



High Temperature 

-J 

I 

Low Temperature 

Low Temperature 

Low Temperature 
Low Temperature 



High Temperature 
Low Temperature 



Low Temperature 
Low Temperature 
Low Temperature 



Continued on next page 



No. 



Reactants 



Products 



Reaction Rate 



— 

-14 
15 

14 (T/298.0) 16 exp (250.0/T) 
12 exp(-1800.0/T) 

12 exp (-160. 0/T) 
- n exp(-1860.0/T) 

13 exp (-5890.0/T) 

- 13 (T/298.0) 1 5 exp (-3890.0/T) 
- 14 (T/298.0) 2 exp (-1290. 0/T) 
- 14 (T/298.0) 1 ' 5 exp (-1810. 0/T) 
11 exp (-500. 0/T) 

11 exp (-390. 0/T) 

12 exp (380.0/T) 

11 exp (250.0/T) 

12 exp (-710. 0/T) 
li 



Reference 


Temperature 


Note 


NIST 


1500 


High Temperature 


NIST 


500-2500 




NIST 


1000 


High Temperature 


NIST 


300-2000 




JPL 


200-300 


Low Temperature 


JPL 


200-300 


Low Temperature 


NIST 


1500-2400 


High Temperature 


NIST 


500-2500 




NIST 


298-2840 




NIST 


500-2500 




NIST 


1250-1860 


High Temperature 


NIST 


300-2000 




JPL 


200-300 


Low Temperature 


JPL 


200-300 


Low Temperature 


JPL 


200-300 


Low Temperature 


JPL 


200-300 


Low Temperature 


NIST 


298 




JPL 


200-300 


Low Temperature 


JPL 


200-300 


Low Temperature 


JPL 


298 




JPL 


200-300 


Low Temperature 


NIST 


880-1150 


High Tempeig^urc 


NIST 


1140-1480 


High Temperature 


NIST 


910-1660 


High Temperature 


NIST 


1050-1540 


High Temperature 


NIST 


295 




NIST 


2420-3870 




NIST 


296-410 


Low Temperature 


JPL 


200-300 


Low Temperature 


JPL 


298 




NIST 


300-2570 




NIST 


1120-1540 


High Temperature 


JPL 


298 




JPL 


200-300 


Low Temperature 


JPL 


298 




Kasting (1990) 




Low Temperature 


JPL 


200-300 


Low Temperature 


JPL 


200-300 


Low Temperature 


NIST 


298 




NIST 


295-665 


Low Temperature 


NIST 


295 




NIST 


298 




NIST 


900-1800 


High Temperature 



R508 OH + CNO 

R509 OH + CNO 

R510 OH + CNO 

R511 OH + CO 

R512 OH + H 2 

R513 OH + H 2 2 

R514 OH + HCN 

R515 OH + HCN 

R516 OH + HCN 

R517 OH + HCNO 

R518 OH + HCNO 

R519 OH + HNO 

R520 OH + HNO a 

R521 OH + HN0 4 

R522 OH + H0 2 

R523 OH + NH 3 

R524 OH + N0 3 

R525 OH + O a 

R526 OH + OH 

R527 OH + S 

R528 OH + SO 

R529 S + C 2 H 6 

R530 S + CH 4 

R531 S + H 2 

R532 S + H 2 S 

R533 S + HS 

R534 S + NO 

R535 S + N0 2 

R536 S + 2 

R537 S + 3 

R538 S + OCS 

R539 S + S0 2 

R540 CS + 2 

R541 SO + N0 2 

R542 CS + NO a 

R543 SO + CHO 

R544 SO + 2 

R545 SO + 3 

R546 S0 2 + C 2 H 

R547 S0 2 + CH 

R548 S0 2 + CH 2 

R549 S0 2 + CN 

R550 SQ 2 + H 



CHO + NO 
HCNO + O 
CO + NO + H 
C0 2 + H 
H 2 + H 
H 2 + H0 2 
HCNO + H 
CO + NH 2 
CN + H 2 
CNO + H 2 
CO, + NH, 



H 2 0- 



NO 
NO, 



w 2 

-NH, 



H0 2 + N0 2 
H0 2 + 2 
H 2 + 
SO 
+ H 



H + 

so, 



C 2 H 5 + HS 
CH 3 + HS 
H + HS 
HS + HS 
H + S 2 
SO + N 
NO + SO 
SO + O 
o 2 + SO 
co + s 2 
so + so 
OCS + o 
S0 2 + NO 
OCS + NO 
HSO + CO 
so 2 + O 
so 2 + o 2 

CO + OCS + 
CHO + SO 
CH 2 + SO 
CNO + SO 
OH + SO 



4.0 x 10 
5.38 x 10 

9.3 x 10- 

5.4 x 10- 

2.8 x 10" 

2.9 x 10" 
4.65 x 10 

1.1 x 10" 
1.84 x 10 
9.43 x 10 
1.54 x 10 
8.0 x 10- 
1.8 x 10" 

1.3 x 10" 
4.8 x 10" 
1.7 x 10" 
2.0 x 10" 

1.7 x 10" 

1.8 x 10" 

6.6 x 10" 

2.7 x 10" 

2.04 x 10 

3.4 x 10" 

2.24 x 10 
1.38 x 10 

4.0 x 10" 
1.75 x 10 

4.9 x 10" 

2.3 x 10" 

1.2 x 10" 

1.5 x 10" 
9.77 x 10 
2.9 x 10" 

1.4 x 10" 

7.6 x 10" 

5.5 x 10" 

1.25 x 10 
3.4 x 10" 

1.1 x 10" 
1.32 x 10 
4.9 x 10" 
4.4 x 10" 
4.58 x 10 



(T/298.0) 2 ' 27 exp (-497.0/T) 



exp (-940.0/T) 



12 
12 
11 

n exp(355.0/T) 

" 10 cxp(-7420.0/T) 

1() cxp(-10000.0/T) 

" 10 cxp(-9710.0/T) 

" 10 exp(-3720.0/T) 
n 

" 10 cxp(-20206.0/T) 

11 exp (84.2/T) 
12 

11 

12 exp(-1260.0/T) 
" 12 exp (-4546. 0/T) 
19 

11 

17 

11 (T)" - 4 

" 13 exp(-2190.0/T) 

12 cxp(-1100.0/T) 
n 

-10 

12 

12 

2.3 



1 exp (250.0/T) 



(T/298.0)" 2 - 3 exp (-15635.0/T) 
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No. 



Rcactants 



Products 



Reaction Rate 



Reference 



Temperature Note 



R551 S0 2 + H0 2 

R552 S0 2 + NO a 

R553 S0 2 + O 

R554 S0 3 + H 

R555 S0 3 + O 

R556 S0 3 + OH 

R557 SO ;j + SO 

R558 HSO + 2 

R559 SO + SO 

R560 HS + H 

R561 CHO + CHO 

R562 S0 2 + 2 

R563 SOi + SO, 



R564 SO 3 , 



SO, 



R565 S0 3 +H 2 

R566 HSO3 + OH 

R567 HSO3 + H 

R568 N0 2 + CHO 

R569 S + CHO 

R570 S + H0 2 

R571 S + H0 2 

R572 S + C0 2 

R573 HSO + NO 

R574 HSO + OH 

R575 HSO + H 

R576 HSO + H 

R577 HSO + HS 

R578 HSO + O 

R579 HSO + S 

R580 H + C 3 H 2 

R581 H + C 3 H B 

R582 H + C 3 H 5 

R583 H + C 3 H B 

R584 H + C 3 H 6 

R585 H + C 3 H 7 

R586 H + C 3 H 7 

R587 H + C 4 H 3 

R588 H + C 4 H 3 

R589 H ■ 

R590 H ■ 

R591 H + C 6 H 3 

R592 H + C 6 H 3 

R593 CH + C 2 H 2 



-C 4 H 5 
- C 4 H Q 



HS0 2 + 2 
SO3 + NO 



SO- 
OH 



o 2 

SO, 



H0 2 + S0 2 

so 2 + so 2 

OH 

S0 2 + £ 
H, + S 



SO, 



+ 

so 
so 



so 3 

hSOg 

hSO, 



H 2 + S0 3 
H 2 + S0 3 
CO + HNO, 



HS + 
HS + 
SO + 
SO + 
HNO 



CO 

o 2 

OH 

CO 
+ SO 



H 2 + SO 
HS + OH 
H 2 + SO 
H,S + SO 
hSO 
SO 



OH 

HS- 

h + c 3 nl 

C 3 H\ + H 



-'3 ii 4 

CH 4 + C 2 H ; 

C 3 H 5 + H 2 



C 3 H 6 



+ H 2 

+ CH 3 
+ C 2 H 2 
+ H, 



H, + H, 



8.6 x 10 

1.05 x 10 

8.3 x 10" 
1.46 x 10 

1.06 x 10 
9.73 x 10 
2.0 x 10" 
1.69 x 10 
3.5 x 10" 
2.16 x 10 
3.0 x 10" 
1.0 x 10 _ 
4.0 x 10 _ 
7.0 x 10- 

1.2 x 10- 
1.0 x 10- 
1.0 x 10- 
3.49 x 10 
5.0 x 10- 
1.5 x 10" 
1.5 x 10" 
1.0 x 10- 

3.7 x 10" 

4.8 x 10- 
7.29 x 10 
6.48 x 10 
1.0 x 10" 
3.0 x 10- 
1.0 x 10- 
1.0 x 10- 
6.0 x 10" 
6.0 x 10" 
6.0 x 10- 

4.4 x 10" 
5.0 x 10- 
6.0 x 10- 

3.3 x 10" 
6.0 x 10- 
6.0 x 10- 

1.5 x 10- 
3.3 x 10- 
6.0 x 10- 
3.5 x 10- 



lo exp(-5230.0/T) 
- n exp (-13591. 0/T) 
12 cxp (-9800.0/T) 
- 11 (T/298.0) 1 22 cxp (-1672. 0/T) 
(T/298.0) 2 ' 57 exp (-14673.0/T) 
14 (T/298.0) 2 ' 46 exp (-13711. 0/T) 



-13/ 



(T/298.0) 3 - 29 exp (-1180.0/T) 



15 
-15 
15 
-11 
11 
16 
12 
14 
15 
11 
11 
-24 
11 
11 
11 
20 

12 exp (250.0/T) 

n cxp(250.0/T) 
-11 

-12 

12 

n cxp(200.0/T) 
11 

n cxp(-1000.0/T) 
12 

12 

12 

13 (T/298.0) 2 ' 5 exp (-1250. 0/T) 
11 

11 

12 

12 

12 

12 

12 

12 
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'exp (61. 0/T) 



NIST 
NIST 
NIST 
NIST 
NIST 
NIST 
NIST 
NIST 
NIST 
NIST 
NIST 

Turco ct al. (1982) 
Turco ct al. (1982) 
Turco et al. (1982) 
NIST 

Kasting (1990) 
Kasting (1990) 
NIST 

Kasting (1990) 
Kasting (1990) 
Kasting (1990) 
Yung & Demore (1999) 
Kasting (1990) 
Kasting (1990) 
Kasting (1990) 
Kasting (1990) 
Kasting (1990) 
Kasting (1990) 
Kasting (1990) 
Yung & Demore (1999) 
Yung & Demore (1999) 
Yung & Demore (1999) 
Yung & Demore (1999) 
NIST 

Yung & Demore (1999) 
Yung & Demore (1999) 
Yung & Demore (1999) 
Yung & Demore (1999) 
Yung & Demore (1999) 
NIST 

Yung & Demore (1999) 
Yung & Demore (1999) 
Yung & Demore (1999) 



400-2000 

703-1850 

440-3000 

700-2000 

1000-2000 

800-2000 

298 

298 

298 

295 

300-2500 



298 



294-1650 



300-2500 



300-2500 



Low Temperature 
Low Temperature 
Low Temperature 

Low Temperature 
Low Temperature 

Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 

Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 

Low Temperature 
Low Temperature 
Low Temperature 



No. 



Reactants 



Products 



Reaction Rate 



Reference 



Temperature 



Note 

Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 



Low Temperature 
00 

Low Tempergjjirc 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 



Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 



R594 
R595 
R596 
R597 
R598 
R599 
R600 
R601 
R602 
R603 
R604 
R605 
R606 
R607 
R608 
R609 
R610 
R611 
R612 
R613 
R614 
R615 
R616 
R617 
R618 
R619 
R620 
R621 
R622 
R623 
R624 
R625 
R626 
R627 
R628 
R629 
R630 
R631 
R632 
R633 
R634 
R635 
R636 



CH- 
CH 2 
CH„ 



2 H„ 

+ H 2 
+ H 2 



CH 2 i + CH 4 



CH 2 - 
CH 2 
CH 2 
CH 3 
CH, 



CH 3 + C 3 H 7 
CH., + C o Ho 



C 2 + CH 4 
C 2 H + C 2 H 2 
C 2 H + C 2 H 4 
NO + CN 



C 2 H- 



C 6 H 2 



C„H + C„H„ 



C 2 H 5 



C 3 H 2 


+ C 2 H. 


C3H3 


+ C 2 H. 


C3H5 


+ H 2 


C 3 H 7 


+ H 2 


C 4 H4 


H 2 


C 4 H4 


-CH 4 


C 4 H4 


C 2 H 2 


C 4 H4 


C 2 H 6 


C 4 H4 


- C 4 H 2 


C 4 H 4- C 6 H 2 


C 4 H4 


- C 8 H 2 


C 4 H 5 


+ H 2 


C 4 H 5 


+ C 2 H. 


C 6 H4 


H 2 


C 6 H4 


-CH 4 


C 6 H4 


" C 2 H 2 


C 6 H 4- C 2 H 6 


C 6 H4 


' C 4 H 2 



C 3 H^ + H 
C 3 H 6 + H 
CH 4" H 



CH 3 
CH, 



H 

CH 4 
CH 3 4- CH 
C 3 H 2 ■ 

C 3 H 3 



3 

f H 2 
+■ H 
+ H 

C 3 H 4 
C 3 H 4 



CH 4 
CH 4 
CH 4 4- C 3 H 6 
CH 4 4- C 3 H 7 
CH 4 4- C 4 H 4 
C 



CH, 



C 4 H 2 + H 
C 4 H 4 + H 
CNO + N 



C K H 



- C 3 H 7 
H 



C 8 H 2 + H 
Haze 



C 3 H 6 



+ CH 2 
+ CH 3 
+ H 
+ H 
+ H 
+ CH 3 
+ H 



C 8 H 2 + H 

Haze 

Haze 



C 4 Hg 

C 6 H 6 
C 6 H 2 
C 6 H 2 
H, 



-'8 1 - 1 2 



C 6 H 2 
Haze 



+ H 
+ H 
+ H 
+ CH 3 
+ H 
+ C 2 H B 



2.23 x 10 

1.8 x 10" 
1.26 x 10 

9.24 x 10 

1.2 x 10" 

5.9 x 10" 
1.0 x 10" 
1.0 x 10" 
6.5 x 10" 
6.5 x 10" 
6.5 x 10" 
1.85 x 10 

1.3 x 10" 
6.5 x 10" 

5.5 x 10" 

1.25 x 10 
2.0 x 10" 

1.6 x 10" 
6.0 x 10" 
1.5 x 10" 
1.0 x 10" 
3.75 x 10 
3.75 x 10 
1.8 x 10" 

1.2 x 10" 
5.0 x 10- 
2.0 x 10" 
1.39 x 10 
3.19 x 10 

6.3 x 10" 

9.0 x 10" 
1.5 x 10" 

2.1 x 10" 
3.75 x 10 
3.75 x 10 
3.75 x 10 
1.8 x 10" 
3.0 x 10" 
6.3 x 10" 

9.0 x 10" 
3.75 x 10 

2.1 x 10" 
3.75 x 10 



- 10 exp(173.0/T) 

10 cxp (132.0/T) 
-11 

-11 

11 

11 

10 cxp (-3332.0/T) 
10 cxp (-3332.0/T) 

13 
13 
13 

~ 12 exp(339.0/T) 

24 (T) 4 -°exp(-4175.0/T) 
13 

n cxp(-297.0/T) 

" 10 (T/298.0)°- 24 cxp(37.3/T) 
11 



10 exp (-21168.0/T) 
12 



10 
11 
-12 

-12 
11 
11 



exp (-540.0/T) 



1450.0/T) 
■250.0/T) 



13 
13 

" 13 (T/298.0) 2 ' 3S cxp (-9560.0/T) 
" 14 (T/298.0) 2 - 84 exp(-4600.0/T) 
12 cxp(- 
12 exp (- 
10 
11 
-12 
-12 
-12 

n cxp(-1699/T) 
12 exp(-1861/T) 
12 cxp(-1450.0/T) 

12 cxp (-250.0/T) 

-12 

11 

-12 



Yung & Dcmorc 
Yung & Dcmorc 
Yung & Demore 
Yung & Demore 
Yung & Dcmorc 
Yung & Dcmorc 
Yung & Demore 
Yung & Demore 
Yung & Demore 
Yung & Demore 
Yung & Demore 
Yung & Dcmorc 
Yung & Dcmorc 
Yung & Demore 
Yung & Demore 
NIST 
NIST 
NIST 
NIST 

Yung & Demore 
NIST 

Yung & Demore 
Yung & Demore 
Yung & Demore 
Yung & Demore 
Yung & Demore 
Yung & Demore 
NIST 
NIST 

Yung & Demore 
Yung & Dcmorc 
Yung & Dcmorc 
Yung & Demore 
Yung & Demore 
Yung & Demore 
Yung & Demore 
Yung & Demore 
Yung & Demore 
Yung & Demore 
Yung & Demore 
Yung & Demore 
Yung & Demore 
Yung & Demore 



(1999) 
(1999) 
(1999) 
(1999) 
(1999) 
(1999) 
(1999) 
(1999) 
(1999) 
(1999) 
(1999) 
(1999) 
(1999) 
(1999) 
(1999) 



(1999) 

(1999) 
(1999) 
(1999) 
(1999) 
(1999) 
(1999) 



(1999) 
(1999) 
(1999) 
(1999) 
(1999) 
(1999) 
(1999) 
(1999) 
(1999) 
(1999) 
(1999) 
(1999) 
(1999) 
(1999) 



200-2000 
300-2500 
500-2500 
300-2500 

300-2500 



300-2500 
300-2500 
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No. 



Rcactants 



Products 



Reaction Rate 



Reference 



3.75 X 10" 12 
3.75 x 1CT 12 

1.17 x 1CT 11 cxp (-1260.0/T) 
3.44 x 10- n (T/298.0)-°' 6 
2.66 x 1CT 12 

i.9 x kt 11 

5.0 x 10~ 12 

3.32 x Kr 13 (T/298.0)°' 4 cxp(-250.2/T) 
2.0 x 10- 12 



Temperature Note 



R637 
R638 
R639 
R640 
R641 
R642 
R643 
R644 
R645 



C 6 H- 



6 H 2 

CgH 2 
N 2 H 4 + H 
N0 2 + CH 3 

n 2 h 3 + h' 



NH 2 + OH 
Q 2 + CHQ 2 



Haze 
Haze 



CH 3 N0 2 
NH 9 + NH, 



NH 3 + NH 3 + N 2 



NH 3 
HO, 



O 

CO„ 



Yung , 

Yung , 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 



Demore (1999) 
Demore (1999) 



5.46 x 10 _31 (T/298.0) -1 - 6 x N 
6.89 x 10~ 32 x N 

1.68 x 10 - 24 (T/298.0)- 7 ' cxp (-1390. 0/T) x N 
k = 1.2 x 10- 12 (T/300.0) 1:1 

k x =4.0 x 10- 31 (T/300.0)- 3 ' 6 

4.1 x 10- 30 (T/298.0)- 2 :1 x N 

7.83 x 10- 29 (T/298.0)- 7 ' 56 exp (-5490.0/T) x N 

3.31 x 10" 30 cxp (-740.0/T) x N 

1.38 x 10" 30 x N 

7.69 x 10~ 30 cxp (-380.0/T) x N 

3.21 x 10" 30 (T/298.0)" 2 ' 57 cxp(-215.0/T) x JV 

9.35 x 10- 30 (T/298.0)- 2 -°cxp(-521.0/T) x N 
5.29 x 10~ 34 exp (-370.0/T) x N 

6.04 x 10- 33 (T/298.0)- 1 x N 
3.0 x 10" 30 x N 

k = 2.44 x 10- 10 (T/298.0)-°- 41 

koa = 1.34 x 10- 31 (T/298.0)- 1 - 32 cxp (-370.5/T) 

4.36 x 10- 32 (T/298.0)- 10 x N 
k = 7.5 x 10- n (T/300.0)°' 2 
k x = 4.4 x 10- 32 (T/300.0)- 1 3 
6.87 x 10- 31 (T/298.0)- 2 ' x N 
1.7 x 10~ 33 exp (1000.0/T) x N 
k = 2.9 x 10- 12 (T/300.0)- 1:1 
koo = 2.0 x 10- 31 (T/300.0)- 3 ' 4 
9.18 x 10- 34 (T/298.0)- 1 69 x N 
9.4 x 10~ 33 x N 

5.0 x 10" 32 x N 
1.0 x 10" 36 x N 
1.38 x 10" 33 cxp(502.7/T) x N 
5.46 x 10- 33 (T/298.0)- 1 29 x N 
k = 3.8 x 10- n (T/300.0)-°' 6 
koo = 2.3 x 10- 29 (T/300.0)- 2 ' 8 
k = 1.9 x 10- 11 (T/300.0)- 1 - 8 
k x = 5.3 x 10- 29 (T/300.0)- 4 ' 4 



222-657 

300-1400 

300 

298 

423 

250-3000 
220-430 



Low Temperature 
Low Temperature 
Low Temperature 



Low Temperature 



Ml 

M2 



C + C 
C + H 2 



M3 


CH 3 + CH 3 


M4 


CH 3 + 2 ' 


M5 


CO + CH 


M6 


CO + CH 3 


M7 


H + C2H2 


M8 


H + C 2 H 2 


M9 


H + C 2 H 4 


M10 


H + CHO 


Mil 


H + CN 


M12 


H + CO 


M13 


H + H 


M14 


H + NH 2 


M15 


H + NO 


M16 


H + O 


M17 


H + 2 


M18 


H + OH 


M19 


H0 2 + H0 2 


M20 


H0 2 + N0 2 


M21 


HS + 2 


M22 


N + C 


M23 


N + H 


M24 


N + H 2 


M25 


N + N 


M26 


N + O 


M27 


NO + CH 3 


M28 


N0 2 + CH 3 



CH, 



CH 3 2 



C 2 HO 
C 2 H 3 

C 2 H 3 



HCN 

CHO 

H 2 

NH 3 

HNO 



OH 

HO, 



O, 



H 2 2 
HN0 4 

HSO a 
■ CN 

NH 

NH 2 
■N 2 

NO 



CH 3 N0 3 



NIST 
NIST 
NIST 

JPL 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 
NIST 

JPL 

NIST 

JPL 

JPL 

NIST 

NIST 

NIST 

NIST 

JPL 

NIST 

JPL 

JPL 



5000-6000 
300 

300-2000 

200-300 

53-584 

300-1700 

200-400 

1000-1700 

300-800 

1560-2270 

500-2500 

300-2500 

300-2500 

298 

230-4200 
300-2500 

200-300 

300-3000 

200-300 

200-300 

200-400 

298 

298 

298 

90-6400 
196-297 

200-300 

200-300 



High Temperature 



Low Temperature 
Low Temperature 

Low Temperature 
High Temperature 



00 
1— 1 

I 



Low Temperature 

Low Temperature 

Low Temperature 
Low Temperature 



Low Temperature 
Low Temperature 
Low Temperature 
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No. 


Reactants 


M29 


JN (J 2 + IN U 3 


M30 


o + c 


M31 


o + co 


M32 


O + NO 


M33 


+ N0 2 


M34 


O + O 


M35 


o + o 2 


M36 




M37 


OH + C 2 H 2 


M38 


OH + C 2 H 4 



Products 



Reaction Rate 



Reference 



Tempcraturc Note 



M39 
M40 



M44 
M45 
M46 
M47 
M48 



M51 
M52 
M53 
M54 
M55 

M56 
M57 



NH 2 + C 2 H 3 
OH + HNO a 



M41 OH + NO 
M42 OH + N0 2 
M43 OH + OH 



s + s 
so + o 

SO + OH 
S0 2 +H 

so 2 + 



M49 S0 2 + OH 
M50 H + SO 



s + s 2 
S + S3 

^2 + ^2 

S4 + S 4 
H + CH 2 

H + CH 3 
H + C 2 H 5 



N 2 5 



CO 

co 2 

NO, 



NO., 



o 2 
•o. 



C 2 H s O 



C 2 H 5 N 
H 2 + NO, 



HNO, 



HNO, 



so 2 

HS0 2 
HSO a 
SO., 



HSO3 
HSO 

S, 



Sg(A) 
CH 3 

CH 4 

C 2 H 6 



1.4 x 10~ 12 (T/300.0) 
2.0 x 10" 3(, (T/300.0)" 



ko 
koo 

2.0 x 10" 34 x N 

1.7 x 10- 33 exp(-1509.0/T) x N 
k = 3.0 x 10- 11 

koo = 9.0 x 10- 32 (T/300)- 1 5 
k = 2.2 x 10- n (T/300.0)-°- 7 
fcoo = 2.5 x 10- 31 (T/300.0)- 1 ' 8 
5.21 x 10~ 35 exp (900.0/T) x N 
6.0 x 10- 34 (T/300.0)- 2 ' 4 x N 

2.8 x 10- 36 (T/300.0)-°- 9 x N 
k = 8.3 x 10- 13 (T/300.0) 2 ' 
A; 00 = 5.5 x 10" 30 

k = 7.5 x 10- 12 (T/300.0)-°- 85 
k x = 1.0 x 10- 2S (T/300.0)- 4 ' 5 
1.8 x 10" 27 (T/298.0)- 3 ' 85 x N 
k = 2.4 X 10~ 14 cxp (460.0/T) 
k 2 = 2.7 X 10- 17 cxp(2199.0/r) 
k 3 = 6.5 x 10~ 34 cxp (1355.0/T) 
k = 3.6 x 10- n (T/300.0)-° 1 
k^ = 7.0 x 10- 31 (T/300.0)- 2 ' 6 
fcrj = 2.8 X 10- 11 
k x = 1.8 x 10- 30 (T/300.0)- 3 
k = 2.6 x 10- 11 



: 6.9 x 10- 31 (T/300.0) 



-1.0 



5.0 x 10~ 32 cxp (900.0/T) x N 
4.82 x 10- 31 (T/298.0)- 2 ' 17 x N 

6.45 x 10" 29 (T/298.0)- 3 ' 48 cxp (-490.0/T) x N 

5.61 x 10" 30 (T/298.0)- 5 ' 19 exp (-2271.0/T) x N 

k = 4.2 x 10- 14 (T/300.0) 1 8 

fcoo = 1.8 x 10- 33 (r/300.0) 2 () 

ko = 1.6 x 10- 12 

koa = 3.3 x 10- 31 (T/300.0)- 4 ' 3 

fe = 7.5 x 10- 11 

k x = 5.7 x 10- 32 (T/300.0)- 1 - 6 

1.1 x 10" 30 (T/300.0)- 2 '° x N 
x N 



-2.0 



1.1 x 10- 30 (T/300.0) 
4.0 x 10~ 31 exp (900.0/T) x N 
4.0 x 10~ 31 cxp (900.0/T) x N 
k = 2.37 x 10- 12 cxp (523.0/T) 
k x = 5.8 x 10" 30 cxp (355.0/T) 

6.2 x 10" 29 (T/298.0) 1 8 x N 
k = 1.5 x 10- 10 



JPL 

NIST 

NIST 

JPL 

JPL 
NIST 
JPL 
JPL 

JPL 

JPL 



JPL 
JPL 
JPL 
JPL 

Moses et al. (2002) 

NIST 

NIST 

NIST 

JPL 

JPL 

Kasting (1990) 

Moses et al. (2002) 

Moses et al. (2002) 

Moses et al. (2002) 

Moses et al. (2002) 



Yung , 
NIST 



Demore (1999) 



200-300 

8000 

300-2500 

200-300 

200-300 
200-4000 
200-300 
200-300 

200-300 

200-300 
200-400 



200-300 

200-300 

200-300 

200-300 

100-500 

298-3000 

200-2000 

200-1000 

200-300 

200-300 



300-1000 



Low Temperature 
High Temperature 



Low Temperature 

Low Temperature 

Low Temperature 
Low Temperature 

Low Temperature 

Low Temperature 
Low Temperature 

Low Temperature 
Low Temperature 

00 
to 

Low Temperature 

Low Temperature 
Low Temperature 



Low Temperature 

Low Temperature 

Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 
Low Temperature 

Low Temperature 
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No. 



Reactants 



Products 



Reaction Rate 



Reference 



-,-23/ 



\0.22 



Temperature Note 



M58 


H + C a H 2 


M59 


H + C3H3 


M60 


H + C 3 H, 

H + c 3 n\ 


M61 


M62 


H + C 3 H* 


M63 


H + CgH^ 


M64 


H + c 3 h 2 


M65 


H + C 3 H 5 


M66 


H + C 3 H 6 


M67 


H + C 3 H 7 


M68 


H + C 3 H 8 


M69 


H + C 4 H 2 



M70 H + C 4 H 5 
M71 



H + C fi H, 



M72 


CH + 


H 2 


M73 


CH 3 -{ 


-C 2 H 3 


M74 


CH 3 A 


" C 2 H 5 


M75 


CH 3 -i 


" C 3 H 3 


M76 


CH 3 -i 


" C 3 H 3 


M77 




-C 3 H 5 


M78 


ch^ 


"C 3 H 6 


M79 




-C 3 H 7 


M80 


C 2 H 3 


+ C 2 H 2 


M81 


C 2 H 3 


+ C 2 H 3 


M82 


C 2 H 3 


+ C 2 H B 


M83 


C 2 H 5 


+ C 2 H 5 


M84 


NH 2 - 


-NH 2 


M85 


NH 2 - 


-CH 3 


Tl 


O3 




T2 


H0 2 




T3 


H 2 2 




T4 


N 2 




T5 


N0 2 




T6 


N0 3 




T7 


N 2 O s 




T8 


HNO 




T9 


HNO- 




T10 


HNOg 



C 3 H 3 
C 3 H 4 
C 3 H 4 
CH 3 + 

C 3 H 5 

CH, + 



C 3 H 6 

C 3 H 7 



C 3 H 7 



C 4 Hg 
C„H, 



CH., 



C 3 H 6 



C 4 Hg 



C 4 H 10 



C 4 H 6 



CH C N 



exp (43.0/T) 



J (T/298.0)° 
- 1() (T/298.0) ' 1 exp(16.0/T) 
- 1() (T/298.0)°- 21 exp(86.6/T) 



4.42 x lCT 10 ^ 
1.9 x 10" 
1.1 x 1CT 

9.49 x 10- 12 (T/298.0) 1 ' 74 cxp(3872.8/T) 

9.61 x 10- 12 cxp(-1560.0/T) 

9.49 x Kr 12 (T/298.0) 1 ' 74 cxp(3872.8/T) 

1.01 x 10" 1:L (T/298.0) - 69 exp (1509.0/T) 
2.64 x 10- 10 (T/298.0) ' 18 cxp(62.5/T) 
7.59 x Kr 12 (T/298.0) ' 51 exp (-1318.2/T) 
3.32 x KT 11 

4.23 x 10- 12 (T/298.0) 2 ' 54 cxp(-3400/T) 
k = 1.39 X 1(T 10 exp (-1184/T) 
fcoo = 1.0 x 10" 28 

6.68 x 10- n (T/298.0) ' 31 exp(93.8/T) 
k = 1.39 X 10~ 10 exp (-1184/T) 
fcoo = 1.0 x 10- 28 
fco = 2.37 x 10" 12 exp (523.0/T) 
k x = 5.8 x 10" 30 exp (355.0/T) 

1.2 x 10- 10 
5.6 x lO" 11 

8.3 x 10~ 12 
8.3 x 10~ 12 

2.72 x 10" n (T/298.0)" ' 32 exp(66.1/T) 
1.27 x 10" 14 (T/298.0) 2 ' 67 exp(-3447.0/T) 
5.16 x 10" 1:L (T/298.0)- - 32 
2.09 x 10- 14 (T/298.0) 1 ' 9 cxp(-1060/T) 
fe = 8.2 x 10- 11 

feoo = 8.76 x 10- 6 (T)- 7 03 exp (-1390/T) 
6.5 x lO" 11 
1.9 x 10" 11 

1.96 x 10" 29 (T/298.0)" 3 ' 9 x N 
1.8 x 10" 27 (T/298.0)- 3 ' 85 x N 



Yung & Demore (1999) 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

Yung & Demore (1999) 
NIST 

Yung & Demore (1999) 

Yung & Demore (1999) 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 

NIST 



200-2000 

200-2000 

200-2000 

200-3000 

195-503 

200-3000 

350-1200 

200-2000 

298-2000 

300-2000 

298-2500 



200-2000 



298 

300-2000 

1200-1500 

1200-1500 

300-2500 

298-2000 

300-2500 

300-1500 



7.16 x 10 _1 ° exp (-11200.0/T) x N 
2.41 x 10" 8 (T/298.0)" 1 ' 18 cxp(-24415.0/T) x N 
2.01 x 10~ 7 exp (-22852. 0/T) x N 
7.07 x 10~ 9 exp (-29467.0/T) x N 
1.88 x 10- 4 Or/298.0)- 3 ' 37 cxp(-37645.0/T) x N 
2.5 x 10 6 exp (-6100.0/T) 

1.0 x 10- 3 (T/298.0)- 3 - 5 exp(-11000.0/T) x N 
5.48 x 10- 7 (T/298.0)- 1 ' 24 cxp(-25312.0/T) x N 
1.98 x 10- 3 (T/298.0)- 3 - 8 cxp(-25257.0/T) x N 
1.15 x 10- 6 cxp(-23092.0/T) x N 



Low Temperature 



Low Temperature 



Low Temperature 



Low Temperature 



Low Temperature 

I 

00 

High Temperature 
High Temperature 



Yung & Demore (1999) 




Low Temperature 


NIST 


298 




NIST 


300-1200 




NIST 


200-400 


Low Temperature 


NIST 


200-400 


Low Temperature 


NIST 


300-3000 




NIST 


200-2000 




NIST 


700-1500 




NIST 


900-2000 




NIST 


1400-2500 




NIST 


298-934 




NIST 


200-400 


Low Temperature 


NIST 


300-2500 




NIST 


300-2500 




NIST 


300-2000 





o 2 + o 

O, + H 



OH 

•N 2 - 
NO 
NO 
NO, 
H + 
OH 



f OH 
-O 

to 
to 2 

+ NO s 
NO 
-I- NO 



OH + NQ 2 
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No. 



Rcactants 



Products 



Reaction Rate 



Reference 


Temperature 


Note 


NIST 


295-1200 




NIST 


245-328 


Low 


NIST 


260-300 


Low ' 


NIST 


1400-3200 




NIST 


1280-2030 




NIST 


1370-2000 




NIST 


249 


Low 


NIST 


298-343 


Low 


NIST 


1400-2500 




NIST 


1000-2000 




NIST 


300-2500 




NIST 


1000-3000 




NIST 


500-800 




NIST 


300-1000 




NIST 


500-2500 




NIST 


1500-3200 




NIST 


700-900 




NIST 


200-2500 




NIST 


1500-2500 




NIST 


1000-2000 




NIST 


627-713 




NIST 


500-2500 




NIST 


500-2000 




NIST 


1000-1700 




NIST 


1000-1700 




NIST 


298-2000 




NIST 


544-673 




NIST 


298-2000 




NIST 


1420-1750 


High 


NIST 


1450-2600 




NIST 


1273-1348 




NIST 


1100-2000 




NIST 


200-2000 




NIST 


200-2000 




NIST 


200-2000 




NIST 


298-1229 




NIST 


298-2000 




NIST 


200-2200 




NIST 


200-2200 




NIST 


680-955 




NIST 


440-473 




NIST 


700-1200 




NIST 


150-1500 





Til 
T12 
T13 
T14 CH 2 6 
T15 CH 2 2 
T16 CH 2 2 
T17 
T18 
T19 
T20 
T21 
T22 
T23 
T24 
T25 
T26 
T27 
T28 
T29 



HN0 3 
HN0 4 
HNO, 



CH 3 2 



CH 4 



C 2 H 5 



T31 
T32 
T33 
T34 
T35 
T36 

T37 C 2 H s O 
T38 C 2 H s O 
T39 HCNO 
T40 CNO 
T41 SO s 
T42 HSO 
T43 HSO a 
T44 
T45 
T46 
T47 
T48 
T49 

T50 CH 3 N0 2 
T51 CH,NO, 
T52 

T53 HSO 



HS0 2 
HSO a 
CHO 



CHO a 
CHO, 



CH 3 N0 3 



H0 2 

0,4 



-I- NO 
HN0 2 
H0 2 + N0 2 
CO + H 2 
C0 2 + H 2 
CO + H 2 
CH 2 + OH 
CH 2 2 +H 
CH,0 + H 
OH 



CH, 



CH, 



H 2 
+ OH 



CO- 
CO- 



+ H 
+ H 2 
+ H 
+ H 
-CH 
CH 2 
C 2 H 2 + OH 
CO + CH 3 
CHO + CH 3 
C 2 H 2 + H 2 
CH 4 + CO 
CH 2 + CH 3 
C 2 H 4 + OH 
C 2 H 4 + H 
CN + OH 
CO + N 

so 2 + O 

SO + H 
S0 2 +H 
O + HSO 
2 + HS 
CO + H 
CH 2 + H 
CO + OH 



CO, 



H 



CH 3 + NO 
CH 2 + HNO 
CH 3 + N0 2 
OH + SO, 



8.0 x 10~ 2 (T/298.0)~ b " ,J cxp (-26099.0/T) x N 

6.31 x 10 17 cxp(-13110.0/T) 

4.10 x 10" 5 cxp(-10600.0/T) x N 

7.7 x lO" 9 cxp (-33075. 0/T) x N 
2.49 x 10~ 8 exp (-28745.0/T) x N 
3.47 x 10" 10 cxp(-20326.0/T) x N 
5.0 x 10 4 

1.0 x 10 14 cxp(-7500.0/T) 

2.16 x 10~ 8 cxp (-33556. 0/T) x N 

1.1 x 10- 7 cxp(-33075.0/T) x N 

2.12 x 10 13 (T/298.0) 1 ' 22 cxp (-43539.0/T) 
9.51 x 10 15 (T/298.0)- 1 02 cxp (-46185.0/T) 
6.0 x 10 14 cxp(-21288.0/T) 

4.15 x 10 13 exp (-22731. 0/T) 
0.19(T/298.0)- 7 ' 5 cxp(-22852.0/T) X N 

5.8 x 10~ 8 cxp (-35961. 0/T) x N 
1.69 x 10- 6 cxp (-16838. 0/T) x N 

8.11 x 10 17 (T/298.0)- 1 ' 23 cxp (-51356.0/T) 
1.08 x 10~ 8 cxp (-29828. 0/T) x N 

5.98 x 10~ 9 cxp (-29828. 0/T) x N 

2.2 x 10 12 cxp(-15154.0/T) 

6.82 x 10- 3 (T/298.0)- 8 ' 62 cxp(-11300.0/T) x N 

2.0 x 10 15 cxp (-39810. 0/T) 

3.0 X 10 14 cxp (-42216. 0/T) 

1.0 x 10 15 cxp (-42817.0/T) 

1.33 x 10 15 (T/298.0)" 2 ' 02 cxp (-12749.0/T) 

6.19 x 10 11 cxp (-11900.0/T) 

1.07 x 10 14 (T/298.0)-°- 69 cxp (-11187.0/T) 

1.66 x 10- 9 cxp(-64.0/T) x N 
1.69 x 10" 9 cxp (-23453. 0/T) x N 
9.47 x 10~ 7 cxp (-40051. 0/T) x N 
1.4 x 10" 8 cxp(-29467.0/T) x N 

4.16 x 10~ 7 (T/298.0)- 3 - 29 cxp (-9610. 0/T) x N 
4.61 x 10- 4 (T/298.0)- 5 - 87 cxp(-15635.4/T) x N 
2.73 x 10~ 8 (T/298.0)- 2 - 82 exp (3750.0/T) x N 
6.0 x 10- 11 cxp (-7721. 0/T) x N 

1.69 x 10 14 (T/298.0)-°' 39 cxp (-13230.0/T) 

7.67 x 10- fi (T/298.0)- 1 - 89 cxp(-17800.0/T) x N 
1.25 x 10" 5 (T/298.0)" 3 - 02 cxp(-17560.0/T) x N 

1.32 x 10" 6 cxp (-17199. 0/T) x N 
3.98 x 10 13 cxp(-19364.0/T) 

1.0 x 10 13 cxp(-16838.0/T) 

3.16 x 10- 4 (T/298.0)- 4 - 72 cxp (-13591. 0/T) x N 



oo 



High Temperature 
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No. 



Rcactants 

SO^ 

SO2 

SO^ 

HCN 

H 2 S 



ocs 

C 3 H 5 
C 3 H 5 



Products 



Reaction Rate 



Reference 


Tempcrati 


Turco ct al. (1982) 




Turco et al. (1982) 




Turco ct al. (1982) 




NIST 


1800-5000 


NIST 


1890-2890 


NIST 


873-1470 


NIST 


1900-3230 


NIST 


300-2500 


NIST 


300-2500 


NIST 


300-2500 


NIST 


300-2500 


NIST 


300-2500 


NIST 


1200-1800 


NIST 


298-2000 


NIST 


298-1000 


NIST 


300-2500 


NIST 


600-1653 


NIST 


3392 


NIST 


1900-2500 


NIST 


1170-1690 


NIST 


1700-2400 


NIST 


900-1100 


NIST 


300-1500 


NIST 


1100-1600 


NIST 


300-2500 


NIST 


300-2500 


NIST 


1100-1600 


NIST 


300-2500 


NIST 


1560-2210 


NIST 


1100-1600 


NIST 


300-2500 


NIST 


298-2000 


NIST 


298-2000 


NIST 


300-2500 


NIST 


300-2500 


NIST 


300-2500 


NIST 


300-2500 


NIST 


1280-1580 


NIST 


1750-2450 


NIST 


1530-1980 



Note 



T54 
T55 
T56 
T57 
T58 
T59 
T60 
T61 
T62 
T63 
T64 
T65 
T66 
T67 
T68 
T69 
T70 
T71 
T72 
T73 
T74 
T75 
T76 
T77 
T78 
T79 
T80 
T81 
T82 
T83 
T84 
T85 
T86 
T87 
T88 
T89 
T90 
T91 
T92 
T93 



C 3 H 6 



C 3 H 6 



C 3 H 7 

C 3 H S 



C4H5 
C 4 H 6 



C 4 H 9 



C 4 H 10 



CH 5 N 



SO^ 

so 2 
so 2 

H + CN 
■ H 2 + S 
HS + H 
CO + S 
C 3 H 4 + H 



C 2 H 3 + CH 3 
C 3 H, +H 
C 3 H 4 ■ 



■H„ 



C 2 H 4 + CH 3 
C 3 H 7 + H 

C 2 H 5 + CH 3 

C 3 H 6 + H 2 
C 4 H 2 + H 
C 4 H 3 + H 



C 2 H 2 + C 2 H 3 



C 4 H 5 
C 3 H 3 



H 

-CH, 



C 3 H 3 + CH 3 



CH, 



H + C 
CH 3 + 
H + C 
CH, 4 



NH 2 
CH 4 
CH, 



C 3 H 6 



C 3 H 7 



- NH 2 
NH 
NH, 



1.0 x 10~ 12 x N + 1.5 x 10 3 
1.0 x lO" 11 xff + 2.2x 10 4 

1.5 x 10~ 13 x N + 1.13 x 10 3 

1.93 x 10- 4 (T/298.0)- 2 - 44 cxp(-62782.1/T) x N 
3.16 x 10- 10 cxp(-32954.6/T) x N 

2.92 x 10~ 8 cxp (-33315. 4/T) x N 

4.07 x 10~ 10 cxp (-30910.0/T) x N 

3.98 x 10 12 cxp(-18402.0/T) 

1.40 x 10 13 cxp (-30188.3/T) 

1.26 x 10 13 cxp(-16838.0/T) 

1.18 x 10 18 (T/298.0)" 1 ' 2 cxp (-49191. 0/T) 

2.50 x 10 15 cxp (-43659.0/T) 

4.00 x 10 13 cxp (-40291.0/T) 

1.09 x 10 13 (T/298.0)°' 17 exp (-17921. 6/T) 

1.31 x 10 13 (T/298.0)°' 87 cxp(-15274.6/T) 

1.58 x 10 16 cxp (-49071. 0/T) 

1.07 x 10 12 (T/298.0)- 15 ' 74 cxp (-49672. 0/T) x N 
3.38 x 10 10 

1.66 x 10~ 8 cxp(-29707.0/T) x N 
6.3 x 10 13 cxp (-43779. 0/T) 
1.15 x 10 15 cxp (-43539.0/T) 

1.6 x 10 14 cxp(-20807.0/T) 

3.57 x 10- 8 (T/298.0) ' 7 cxp (-21288. 0/T) x N 

7.7 x 10 14 cxp(-44260.0/T) 
5.0 x 10 15 cxp (-38126. 0/T) 
2.0 x 10 15 cxp (-43178. 0/T) 
2.0 x 10 15 cxp (-37645. 0/T) 

1.58 x 10 16 cxp (-55445.0/T) 

1.8 x 10 13 cxp(-42817.0/T) 
1.0 x 10 14 cxp (-37645.0/T) 

1.0 x 10 16 cxp(-36683.0/T) 

1.1 x 10 13 (T/298.0)°' 25 cxp(-17921.6/T) 
7.7 x 10 13 (T/298.0) ' 77 cxp (-15394.8/T) 
1.58 x 10 16 exp (-49312.0/T) 

1.0 x 10 17 cxp(-42576.0/T) 

7.94 x 10 16 cxp (-40411. 0/T) 
6.31 x 10 14 cxp (-31151. 0/T) 
6.61 x 10~ 9 cxp(-20567.0/T) x N 
2.66 x 10" 11 cxp (-24295. 0/T) x N 
1.36 x 10- 7 exp(-30669.0/T) x N 



Low Temperature 
Low Temperature 
Low Temperature 



High Temperature 



00 
01 



